ikrontki.c

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

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

   IkronTkI_VAR = getComplexMatrix(L*N*M);
   TkronI_VAR = getComplexMatrix(N);
   TkronI_VAR = TridekI(N,S,M);

   l=0;
   for(m = 0; m < (L*N*M); ++m) /* este es el valor de filas *columnas*/
   {
    k=0;
	  if (m %(N*M) ==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 < N*M; ++j)/* numero de filas o columnas de B*/
      {
        if((l % M) == i){
	       IkronTkI_VAR[m][k].real =  TkronI_VAR[m % (N*M)][j].real; /* m%filas o columnas de B*/
	       IkronTkI_VAR[m][k].imag = TkronI_VAR[m % (N*M)][j].imag;
	 	   k++;
     	}
     	else{
     	    IkronTkI_VAR[m][k].real = 0.0;
	        IkronTkI_VAR[m][k].imag = 0.0;
	 	    k++;
	 	}
      }
    }
   }
    free(TkronI_VAR);
   return IkronTkI_VAR;
}