/********************************************************/
/*IkronP						*/
/*Maritza Rodr韌uez Mart韓ez				*/
/*Computational Signal Processing Group - march-98   	*/
/*CSPG - Dr. Domingo Rodr韌uez - Coordinator         	*/
/*							*/

/********************************************************/


COMPLEX **IkronP(int M, int N,int S)
{
  COMPLEX **per;
  COMPLEX **IkronP_VAR;
  int temp,l=0,m,k,i,j;
  int MN = M*N;

  per = getComplexMatrix(N);
  IkronP_VAR = getComplexMatrix(MN);


{
	int k, inicio=0;
	int pos=inicio;

for( k=0; k < N; ++k){
   per[k][pos].real=1;
   pos=pos+S;

   if(pos >= N){
      pos=inicio+1;
      inicio=inicio+1;
   }

}
}
l=0;
   for(m = 0; m < (M*N); ++m) /* este es el valor de filas *columnas*/
   {
	 k=0;
	 if (m % N ==0 & m!=0 ) /* el valor para comparar m debe ser el num de columnas de B*/
	   l++;

	 if(M > N)
	    temp = M;
	 else
	    temp = N;

   for (i = 0; i < temp; ++i)/* numero de filas o columnas mayor entre A y B*/
   {
	for ( j = 0; j < N; ++j)/* numero de filas o columnas de B*/
    {
       if((l % M) == i){
          IkronP_VAR[m][k].real = per[m % N][j].real; /* m%filas o columnas de B*/
          IkronP_VAR[m][k].imag = per[m % N][j].imag;
	   	  k++;
      	}
      	else{
      	  IkronP_VAR[m][k].real = 0.0;
          IkronP_VAR[m][k].imag = 0.0;
	   	  k++;
	   	}
	}
  }
   }
free(per);
return IkronP_VAR;
}