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

Add indexes to print

Browse Source
milestone_5
chortas 3 years ago
parent 4a90597609
commit 30a14002ae
1 changed files with 45 additions and 45 deletions
Show Stats Download Patch File Download Diff File

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

@ -190,9 +190,9 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
i = 1; i = 1;
for (j = 1; j <= ntot; j++) { for (j = 1; j <= ntot; j++) {
chunk_array_get(datar, i, &valueri); chunk_array_get(datar, i, &valueri);
printf("[1] datar[i] = %f\n", valueri); printf("[1] datar[%d] = %f\n", i, valueri);
chunk_array_get(datar, i+1, &valueri1); chunk_array_get(datar, i+1, &valueri1);
printf("[2] datar[i1] = %f\n", valueri1); printf("[2] datar[%d] = %f\n", i+1, valueri1);
chunk_array_save(datar, j, valueri); chunk_array_save(datar, j, valueri);
chunk_array_save(datai, j, valueri1); chunk_array_save(datai, j, valueri1);
i += 2; i += 2;
@ -217,10 +217,10 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datar, j3, &valuerj3); chunk_array_get(datar, j3, &valuerj3);
chunk_array_get(datai, j3, &valueij3); chunk_array_get(datai, j3, &valueij3);
printf("[3] datar[i3] = %f\n", valueri3); printf("[3] datar[%d] = %f\n", i3, valueri3);
printf("[4] datai[i3] = %f\n", valueii3); printf("[4] datai[%d] = %f\n", i3, valueii3);
printf("[5] datar[j3] = %f\n", valuerj3); printf("[5] datar[%d] = %f\n", j3, valuerj3);
printf("[6] datai[j3] = %f\n", valueij3); printf("[6] datai[%d] = %f\n", j3, valueij3);
tempr = valueri3; tempr = valueri3;
tempi = valueii3; tempi = valueii3;
@ -256,8 +256,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(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[7] datar[j] = %f\n", valuerj); printf("[7] datar[%d] = %f\n", j, valuerj);
printf("[8] datai[j] = %f\n", valueij); printf("[8] datai[%d] = %f\n", j, valueij);
if (icase != 3) { if (icase != 3) {
workr[i] = valuerj; workr[i] = valuerj;
@ -309,10 +309,10 @@ 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(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[9] datar[i] = %f\n", tempr); printf("[9] datar[%d] = %f\n", i, tempr);
printf("[10] datai[i] = %f\n", tempi); printf("[10] datai[%d] = %f\n", i, tempi);
printf("[11] datar[j] = %f\n", valuerj); printf("[11] datar[%d] = %f\n", j, valuerj);
printf("[12] datai[j] = %f\n", valueij); printf("[12] datai[%d] = %f\n", j, valueij);
chunk_array_save(datar, i, valuerj - tempr); chunk_array_save(datar, i, valuerj - tempr);
chunk_array_save(datai, i, valueij - tempi); chunk_array_save(datai, i, valueij - tempi);
@ -339,22 +339,22 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datai, i, &tempr); chunk_array_get(datai, i, &tempr);
chunk_array_get(datar, i, &tempi); chunk_array_get(datar, i, &tempi);
printf("[13] datai[i] = %f\n", tempr); printf("[13] datai[%d] = %f\n", i, tempr);
printf("[14] datar[i] = %f\n", tempi); printf("[14] datar[%d] = %f\n", i, tempi);
tempi = -tempi; tempi = -tempi;
} else { } else {
chunk_array_get(datai, i, &tempr); chunk_array_get(datai, i, &tempr);
printf("[15] datai[i] = %f\n", tempr); printf("[15] datai[%d] = %f\n", i, tempr);
tempr = -tempr; tempr = -tempr;
chunk_array_get(datar, i, &tempi); chunk_array_get(datar, i, &tempi);
printf("[16] datar[i] = %f\n", tempi); printf("[16] datar[%d] = %f\n", i, tempi);
} }
chunk_array_get(datar, j, &valuerj); chunk_array_get(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[17] datar[j] = %f\n", valuerj); printf("[17] datar[%d] = %f\n", j, valuerj);
printf("[18] datai[j] = %f\n", valueij); printf("[18] datai[%d] = %f\n", j, valueij);
chunk_array_save(datar, i, valuerj - tempr); chunk_array_save(datar, i, valuerj - tempr);
chunk_array_save(datai, i, valueij - tempi); chunk_array_save(datai, i, valueij - tempi);
@ -362,8 +362,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(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[19] datar[j] = %f\n", valuerj); printf("[19] datar[%d] = %f\n", j, valuerj);
printf("[20] datai[j] = %f\n", valueij); printf("[20] datai[%d] = %f\n", j, valueij);
chunk_array_save(datar, j, valuerj + tempr); chunk_array_save(datar, j, valuerj + tempr);
chunk_array_save(datai, j, valueij + tempi); chunk_array_save(datai, j, valueij + tempi);
@ -408,10 +408,10 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datai, i, &valueii); chunk_array_get(datai, i, &valueii);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[21] datar[i] = %f\n", valueri); printf("[21] datar[%d] = %f\n", i, valueri);
printf("[22] datar[j] = %f\n", valuerj); printf("[22] datar[%d] = %f\n", j, valuerj);
printf("[23] datai[i] = %f\n", valueii); printf("[23] datai[%d] = %f\n", i, valueii);
printf("[24] datai[j] = %f\n", valueij); printf("[24] datai[%d] = %f\n", j, valueij);
tempr = valueri * wr - valueii * wi; tempr = valueri * wr - valueii * wi;
tempi = valueri * wi + valueii * wr; tempi = valueri * wi + valueii * wr;
@ -474,8 +474,8 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datar, j, &sr); chunk_array_get(datar, j, &sr);
chunk_array_get(datai, j, &si); chunk_array_get(datai, j, &si);
printf("[25] datar[j] = %f\n", sr); printf("[25] datar[%d] = %f\n", j, sr);
printf("[26] datai[j] = %f\n", si); printf("[26] datai[%d] = %f\n", j, si);
oldsr = 0.; oldsr = 0.;
oldsi = 0.; oldsi = 0.;
@ -487,8 +487,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(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[27] datar[j] = %f\n", valuerj); printf("[27] datar[%d] = %f\n", j, valuerj);
printf("[28] datai[j] = %f\n", valueij); printf("[28] datai[%d] = %f\n", j, valueij);
sr = twowr * sr - oldsr + valuerj; sr = twowr * sr - oldsr + valuerj;
si = twowr * si - oldsi + valueij; si = twowr * si - oldsi + valueij;
@ -560,10 +560,10 @@ 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(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[29] datar[i] = %f\n", valueri); printf("[29] datar[%d] = %f\n", i, valueri);
printf("[30] datai[i] = %f\n", valueii); printf("[30] datai[%d] = %f\n", i, valueii);
printf("[31] datar[j] = %f\n", valuerj); printf("[31] datar[%d] = %f\n", j, valuerj);
printf("[32] datai[j] = %f\n", valueij); printf("[32] datai[%d] = %f\n", j, valueij);
sumr = (valueri + valuerj) / 2.; sumr = (valueri + valuerj) / 2.;
sumi = (valueii + valueij) / 2.; sumi = (valueii + valueij) / 2.;
@ -592,7 +592,7 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
goto L740; goto L740;
for (i = imin; i <= ntot; i += np2) { for (i = imin; i <= ntot; i += np2) {
chunk_array_get(datai, i, &valueii); chunk_array_get(datai, i, &valueii);
printf("[33] datai[i] = %f\n", valueii); printf("[33] datai[%d] = %f\n", i, valueii);
chunk_array_save(datai, i, -valueii); chunk_array_save(datai, i, -valueii);
} }
L740: L740:
@ -608,8 +608,8 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datar, i, &valueri); chunk_array_get(datar, i, &valueri);
chunk_array_get(datai, i, &valueii); chunk_array_get(datai, i, &valueii);
printf("[34] datar[i] = %f\n", valueri); printf("[34] datar[%d] = %f\n", i, valueri);
printf("[35] datai[i] = %f\n", valueii); printf("[35] datai[%d] = %f\n", i, valueii);
chunk_array_save(datar, j, valueri); chunk_array_save(datar, j, valueri);
chunk_array_save(datai, j, -valueii); chunk_array_save(datai, j, -valueii);
@ -621,8 +621,8 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datar, imin, &valuerimin); chunk_array_get(datar, imin, &valuerimin);
chunk_array_get(datai, imin, &valueiimin); chunk_array_get(datai, imin, &valueiimin);
printf("[36] datar[imin] = %f\n", valuerimin); printf("[36] datar[%d] = %f\n", imin, valuerimin);
printf("[37] datai[imin] = %f\n", valueiimin); printf("[37] datai[%d] = %f\n", imin, valueiimin);
chunk_array_save(datar, j, valuerimin - valueiimin); chunk_array_save(datar, j, valuerimin - valueiimin);
chunk_array_save(datai, j, 0.); chunk_array_save(datai, j, 0.);
@ -635,8 +635,8 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datar, i, &valueri); chunk_array_get(datar, i, &valueri);
chunk_array_get(datai, i, &valueii); chunk_array_get(datai, i, &valueii);
printf("[38] datar[i] = %f\n", valueri); printf("[38] datar[%d] = %f\n", i, valueri);
printf("[39] datai[i] = %f\n", valueii); printf("[39] datai[%d] = %f\n", i, valueii);
chunk_array_save(datar, j, valueri); chunk_array_save(datar, j, valueri);
chunk_array_save(datai, j, valueii); chunk_array_save(datai, j, valueii);
@ -649,8 +649,8 @@ void fourt(chunk_array_t* datar, chunk_array_t* datai, int nn[3], int ndim, int
chunk_array_get(datar, imin, &valuerimin); chunk_array_get(datar, imin, &valuerimin);
chunk_array_get(datai, imin, &valueiimin); chunk_array_get(datai, imin, &valueiimin);
printf("[40] datar[imin] = %f\n", valuerimin); printf("[40] datar[%d] = %f\n", imin, valuerimin);
printf("[41] datai[imin] = %f\n", valueiimin); printf("[41] datai[%d] = %f\n", imin, valueiimin);
chunk_array_save(datar, j, valuerimin + valueiimin); chunk_array_save(datar, j, valuerimin + valueiimin);
chunk_array_save(datai, j, 0.); chunk_array_save(datai, j, 0.);
@ -686,8 +686,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(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[44] datar[j] = %f\n", valuerj); printf("[44] datar[%d] = %f\n", j, valuerj);
printf("[45] datai[j] = %f\n", valueij); printf("[45] datai[%d] = %f\n", j, valueij);
chunk_array_save(datar, i, valuerj); chunk_array_save(datar, i, valuerj);
chunk_array_save(datai, i, -valueij); chunk_array_save(datai, i, -valueij);
@ -699,8 +699,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(datar, j, &valuerj);
chunk_array_get(datai, j, &valueij); chunk_array_get(datai, j, &valueij);
printf("[46] datar[j] = %f\n", valuerj); printf("[46] datar[%d] = %f\n", j, valuerj);
printf("[47] datai[j] = %f\n", valueij); printf("[47] datai[%d] = %f\n", j, valueij);
chunk_array_save(datar, i, valuerj); chunk_array_save(datar, i, valuerj);
chunk_array_save(datai, i, -valueij); chunk_array_save(datai, i, -valueij);

Write Preview
Loading…
Cancel
Save