|
|
|
@ -319,8 +319,8 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
|
|
|
|
|
chunk_array_get(datar, j, &valuerj);
|
|
|
|
|
chunk_array_get(datai, j, &valueij);
|
|
|
|
|
|
|
|
|
|
printf("[9] datar[%d] = %f\n", i, tempr);
|
|
|
|
|
printf("[10] datai[%d] = %f\n", i, tempi);
|
|
|
|
|
printf("[9] tempr = %f\n", i, tempr);
|
|
|
|
|
printf("[10] tempi = %f\n", i, tempi);
|
|
|
|
|
printf("[11] datar[%d] = %f\n", j, valuerj);
|
|
|
|
|
printf("[12] datai[%d] = %f\n", j, valueij);
|
|
|
|
|
|
|
|
|
@ -528,9 +528,15 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
|
|
|
|
|
if (j > jmin)
|
|
|
|
|
goto L620;
|
|
|
|
|
|
|
|
|
|
chunk_array_get(datar, j, &valuerj);
|
|
|
|
|
chunk_array_get(datai, j, &valueij);
|
|
|
|
|
|
|
|
|
|
workr[i] = wr * sr - wi * si - oldsr + valuerj;
|
|
|
|
|
worki[i] = wi * sr + wr * si - oldsi + valueij;
|
|
|
|
|
|
|
|
|
|
printf("[85] wr = %f, sr = %f, wi = %f, si = %f, oldsr = %f, datar[j] = %f\n", wr, sr, wi, si, oldsr, valuerj);
|
|
|
|
|
printf("[86] wi = %f, sr = %f, wr = %f, si = %f, oldsi = %f, datai[j] = %f\n", 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);
|
|
|
|
|
|
|
|
|
|