local matrix

src/fluid_problem.c

@@ -355,7 +355,7 @@ static void fluid_problem_assemble_system(FluidProblem *problem, double *rhs, co
             for (int j = 0; j < DIMENSION; ++j) {
               int U = j;
               LOCAL_MATRIX(U,P) += -jw*dphii[j]*phij;
-              LOCAL_MATRIX(P,U) +=  jw*phii*dphij[j];
+              LOCAL_MATRIX(Q,U) +=  jw*phii*dphij[j];
               double utau =0;
               double dphisigmadq = 0;
               double dutaudq = 0;
@@ -374,9 +374,9 @@ static void fluid_problem_assemble_system(FluidProblem *problem, double *rhs, co
               LOCAL_MATRIX(U,U) += jw*mu[ff]*2*0.5*dphiidtau + jw*rho[ff]*phii*(phij/dt+phij/q*divu+udtau/q);
               LOCAL_MATRIX(U,Q) += jw*rho[ff]*phii*(-u[j]*divu*phij/(q*q)-utau*phij/(q*q)+dutaudq/q)+jw*dphisigmadq;
             }
-            LOCAL_MATRIX(P,P) += jw*epsilon*dphiidphij;
-            LOCAL_MATRIX(Q,Q) += -jw*phii*phij;
-            LOCAL_MATRIX(P,Q) += jw*phii*phij/dt;
+            LOCAL_MATRIX(Q,P) += jw*epsilon*dphiidphij;
+            LOCAL_MATRIX(P,Q) += -jw*phii*phij;
+            LOCAL_MATRIX(Q,Q) += jw*phii*phij/dt;
           }
         }
       }