53 |
if (in!=NULL) { |
if (in!=NULL) { |
54 |
++(in->reference_counter); |
++(in->reference_counter); |
55 |
} |
} |
56 |
|
return in; |
57 |
} |
} |
58 |
|
|
59 |
Paso_SystemMatrix* Paso_FCTransportProblem_borrowTransportMatrix(Paso_FCTransportProblem* in) { |
Paso_SystemMatrix* Paso_FCTransportProblem_borrowTransportMatrix(Paso_FCTransportProblem* in) { |
117 |
Paso_Pattern_color(pattern->mainPattern,&(out->num_colors),out->colorOf); |
Paso_Pattern_color(pattern->mainPattern,&(out->num_colors),out->colorOf); |
118 |
|
|
119 |
|
|
120 |
|
#pragma omp parallel for schedule(static) private(i) |
121 |
|
for (i = 0; i < n; ++i) { |
122 |
|
out->lumped_mass_matrix[i]=0.; |
123 |
|
out->row_sum_flux_matrix[i]=0.; |
124 |
|
out->u[i]=0.; |
125 |
|
} |
126 |
|
|
127 |
/* identify the main diagonals */ |
/* identify the main diagonals */ |
128 |
#pragma omp parallel for schedule(static) private(i,iptr,iptr_main,k) |
#pragma omp parallel for schedule(static) private(i,iptr,iptr_main,k) |
129 |
for (i = 0; i < n; ++i) { |
for (i = 0; i < n; ++i) { |