// ------------------------------------------------------------------------
// File: diagint.c
// Date: April 5, 2002
// Description: Generate a DIAGONAL interleaver. Output is the interleaver
// array in a file with name specified at command line
//
// Example: N1=3 N2=5
//
//                   0  1  2  3  4
//                   5  6  7  8  9
//                  10 11 12 13 14
//
//                i   0  1  2  3  4  5  6  7  8  9 10 11 12 13 14
//    interleaver[i]  0  6 12  1  7 13  2  8 14  3  9 10  4  5 11
// ------------------------------------------------------------------------
// This program is complementary material for the book:
//
// R.H. Morelos-Zaragoza, The Art of Error Correcting Coding, Wiley, 2002.
//
// ISBN 0471 49581 6
//
// This and other programs are available at http://the-art-of-ecc.com
//
// You may use this program for academic and personal purposes only. 
// If this program is used to perform simulations whose results are 
// published in a journal or book, please refer to the book above.
//
// The use of this program in a commercial product requires explicit
// written permission from the author. The author is not responsible or 
// liable for damage or loss that may be caused by the use of this program. 
//
// Copyright (c) 2002. Robert H. Morelos-Zaragoza. All rights reserved.
// ------------------------------------------------------------------------

#include <stdio.h>
#include <math.h>
#include <float.h>
#include <limits.h>
#define NMAX 16384

main(int argc, char *argv[])
{

int interleaver[NMAX];
int i, j, position;
int N1, N2;
int N;                                 // Interleaver length
char name1[40];
FILE *fp1;

  // Command line processing
  if (argc != 4)
    {
      printf("Usage %s N1  N2  INT_file\n", argv[0]);
      exit(0);
    }

  sscanf(argv[1],"%d", &N1);
  sscanf(argv[2],"%d", &N2);
  sscanf(argv[3],"%s", name1);

  N = N1*N2;

  fp1 = fopen(name1,"w");

  for (i=0; i<N; i++)
    {
      position = (i*N2) % N + (i/N1+i%N1) % N2 ;
      interleaver[i] = position;
    }    

  for (i=0; i<N; i++)
    {
      fprintf(fp1, "%d\n", interleaver[i]);
      printf("%5d -> %5d\n",i,interleaver[i]);
    }

  fclose(fp1);

}