diff --git a/fftma_module/gen/lib_src/fourt.c b/fftma_module/gen/lib_src/fourt.c index 25faaf9..b0757cb 100644 --- a/fftma_module/gen/lib_src/fourt.c +++ b/fftma_module/gen/lib_src/fourt.c @@ -527,8 +527,13 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int j -= ifp1; if (j > jmin) goto L620; + workr[i] = wr * sr - wi * si - oldsr + valuerj; worki[i] = wi * sr + wr * si - oldsi + valueij; + + printf("[83] Saving in workr the value %f in pos %d\n", wr * sr - wi * si - oldsr + valuerj, i); + printf("[84] Saving in worki the value %f in pos %d\n", wi * sr + wr * si - oldsi + valueij, i); + jmin += ifp2; i++; } @@ -540,8 +545,8 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int for (j2 = i3; j2 <= j2max; j2 += ifp1) { j3max = j2 + np2 - ifp2; for (j3 = j2; j3 <= j3max; j3 += ifp2) { - printf("[68] Saving in datar the value %f in pos %d\n", workr[i], j3); - printf("[69] Saving in datai the value %f in pos %d\n", worki[i], j3); + printf("[68] Saving in datar the value %f in pos %d, i = %d\n", workr[i], j3, i); + printf("[69] Saving in datai the value %f in pos %d, i = %d\n", worki[i], j3, i); chunk_array_save(datar, j3, workr[i]); chunk_array_save(datai, j3, worki[i]);