121 |
out->u[i]=0.; |
out->u[i]=0.; |
122 |
} |
} |
123 |
/* identify the main diagonals */ |
/* identify the main diagonals */ |
124 |
#pragma omp parallel for schedule(static) private(i,iptr,iptr_main,k) |
#pragma omp parallel for schedule(static) private(i,iptr,iptr_main) |
125 |
for (i = 0; i < n; ++i) { |
for (i = 0; i < n; ++i) { |
126 |
iptr_main=pattern->mainPattern->ptr[0]-1; |
iptr_main=pattern->mainPattern->ptr[0]-1; |
127 |
for (iptr=pattern->mainPattern->ptr[i];iptr<pattern->mainPattern->ptr[i+1]; iptr++) { |
for (iptr=pattern->mainPattern->ptr[i];iptr<pattern->mainPattern->ptr[i+1]; iptr++) { |