wip

fluid/fluid_problem.c

@@ -443,8 +443,7 @@ static void fluid_problem_f(const FluidProblem *problem, const double *sol, doub
       f10[(U+i)*n_fields+U+j][j] += -mu*dc[i]/c;
       f10[(U+i)*n_fields+U+i][i] += -mu*dc[j]/c;
       f11[(U+i)*n_fields+U+j][j][i] += mu;
-      // f11[(U+i)*n_fields+U+i][j][j] += mu;
-      f11[(U+i)*n_fields+U+j][i][j] += mu;
+      f11[(U+i)*n_fields+U+i][j][j] += mu;
 
       // SUPG
       double supg = uold[j]*taup;