Commit ca37deae authored by Michel Henry's avatar Michel Henry
Browse files

wip

parent 67f967a6
......@@ -440,8 +440,10 @@ static void fluid_problem_f(const FluidProblem *problem, const double *sol, doub
}
for (int j = 0; j < D; ++j) {
f1[U+i][j] = mu*(du[i][j]-u[i]*dc[j]/c+ du[j][i]-u[j]*dc[i]/c);
// ca ne correspond pas aux bons indices
f10[(U+i)*n_fields+U+j][j] += -mu*dc[i]/c;
f10[(U+i)*n_fields+U+i][i] += -mu*dc[j]/c;
// equivalent a [(U+i)*n_fields+U+j]*D + j != (U+i)*D+j)*n_fields + U+j
f11[(U+i)*n_fields+U+j][j][i] += mu;
f11[(U+i)*n_fields+U+i][j][j] += mu;
......
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