ikronpki.c

java-fft-demo,用java开发的实现fft的平台

/********************************************************/
/* This function performe IkronPkI                      */
/********************************************************/
COMPLEX **IkronPkI(int L,int N,int S,int M)
{
   COMPLEX **PkronI_VAR,**IkronPkI_VAR;
   int temp,l=0,m,k,i,j;

   IkronPkI_VAR = getComplexMatrix(L*N*M);
   PkronI_VAR = getComplexMatrix(N*M);
   PkronI_VAR = PermkI(N,S,M);

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

	temp = L;

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