|
|
|
@ -25,9 +25,9 @@ void build_real(int n[3], int NTOT, chunk_array_t* covar, chunk_array_t* realiza
|
|
|
|
|
double temp;
|
|
|
|
|
|
|
|
|
|
/*decomposition and multiplication in the spectral domain*/
|
|
|
|
|
for (i = 1; i <= n[0]; i++) {
|
|
|
|
|
for (k = 1; k <= n[2]; k++) {
|
|
|
|
|
for (j = 1; j <= n[1]; j++) {
|
|
|
|
|
for (k = 1; k <= n[2]; k++) {
|
|
|
|
|
for (j = 1; j <= n[1]; j++) {
|
|
|
|
|
for (i = 1; i <= n[0]; i++) {
|
|
|
|
|
maille1 = i + (j - 1 + (k - 1) * n[1]) * n[0];
|
|
|
|
|
chunk_array_get(covar, maille1, &temp);
|
|
|
|
|
if (temp > 0.) {
|
|
|
|
|