ssantisi
/
simulacion-permeabilidad
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add more prints..

Browse Source
milestone_5
chortas 4 years ago
parent 0a7f2fe1ba
commit 90f76bb61d
1 changed files with 7 additions and 2 deletions
Show Stats Download Patch File Download Diff File

9
fftma_module/gen/lib_src/fourt.c
Unescape Escape View File

@ -527,8 +527,13 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
j -= ifp1; j -= ifp1;
if (j > jmin) if (j > jmin)
goto L620; goto L620;
workr[i] = wr * sr - wi * si - oldsr + valuerj; workr[i] = wr * sr - wi * si - oldsr + valuerj;
worki[i] = wi * sr + wr * si - oldsi + valueij; 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; jmin += ifp2;
i++; 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) { for (j2 = i3; j2 <= j2max; j2 += ifp1) {
j3max = j2 + np2 - ifp2; j3max = j2 + np2 - ifp2;
for (j3 = j2; j3 <= j3max; j3 += ifp2) { for (j3 = j2; j3 <= j3max; j3 += ifp2) {
printf("[68] Saving in datar the value %f in pos %d\n", workr[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\n", worki[i], j3); 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(datar, j3, workr[i]);
chunk_array_save(datai, j3, worki[i]); chunk_array_save(datai, j3, worki[i]);

Write Preview
Loading…
Cancel
Save