Commit 6bc59071 authored by Michel Henry's avatar Michel Henry
Browse files

wip

parent 80974a1a
......@@ -165,7 +165,7 @@ void fluid_problem_node_force_volume(FluidProblem *problem, const double *soluti
for (int e = 0; e < D; ++e) {
double supg = drag_in_stab?sold[U+e]*problem->taup[iel]:0;
f10[(U+d)*n_fields+U+d][e] += supg*dfdu;
f11[(U+d)*n_fields+P][d][e] += supg*dfddp;
f11[(U+d)*n_fields+P][e][d] += supg*dfddp;
}
}
fe_fields_add_to_local_matrix(problem->fields, f00, f01, f10, f11, sf, dsf, 1, local_matrix);
......@@ -476,7 +476,7 @@ static void fluid_problem_f(const FluidProblem *problem, const double *sol, doub
double pspg = taup/rho;
f1[P][i] += pspg*(f0[U+i]+(problem->drag_in_stab?0:(1-c)*dp[i])) + problem->stab_param*dp[i];
f11[P*n_fields+P][i][i] += pspg*(f01[(U+i)*n_fields+P][i]+(problem->drag_in_stab?0:1-c))+problem->stab_param;
f10[P*n_fields+U][i] += pspg*f00[(U+i)*n_fields+(U+i)];
f10[P*n_fields+U+i][i] += pspg*f00[(U+i)*n_fields+(U+i)];
for (int j = 0; j < D; ++j) {
f11[P*n_fields+U+i][i][j] = pspg*f01[(U+i)*n_fields+(U+i)][j];
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment