1116         PetscScalar Ke[8 * NODE_SIZE][8 * NODE_SIZE], Kcpl[8 * NODE_SIZE][4 * PRMNODE_SIZE];  in THIJacobianLocal_Momentum()  local
1121         PetscCall(PetscMemzero(Ke, sizeof(Ke)));  in THIJacobianLocal_Momentum()
1143 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + dp[2] * …  in THIJacobianLocal_Momentum()
1144 …           Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * eta * 2. * dpl[1] + dp[1] * jw * eta * dpl[0];  in THIJacobianLocal_Momentum()
1145 …           Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * eta * 2. * dpl[0] + dp[0] * jw * eta * dpl[1];  in THIJacobianLocal_Momentum()
1146 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + dp[2] * …  in THIJacobianLocal_Momentum()
1148 …Ke[l * 2 + 0][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d…  in THIJacobianLocal_Momentum()
1149 …Ke[l * 2 + 0][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * d…  in THIJacobianLocal_Momentum()
1150 …Ke[l * 2 + 1][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d…  in THIJacobianLocal_Momentum()
1151 …Ke[l * 2 + 1][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * d…  in THIJacobianLocal_Momentum()
1153               Ke[l * 2 + 0][ll * 2 + 0] += pp * jw * thi->inertia * pp;  in THIJacobianLocal_Momentum()
1154               Ke[l * 2 + 1][ll * 2 + 1] += pp * jw * thi->inertia * pp;  in THIJacobianLocal_Momentum()
1169             Ke[0][0] = thi->dirichlet_scale * diagu;  in THIJacobianLocal_Momentum()
1170             Ke[0][1] = 0;  in THIJacobianLocal_Momentum()
1171             Ke[1][0] = 0;  in THIJacobianLocal_Momentum()
1172             Ke[1][1] = thi->dirichlet_scale * diagv;  in THIJacobianLocal_Momentum()
1188 …               Ke[l * 2 + 0][ll * 2 + 0] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * u * u * ppl;  in THIJacobianLocal_Momentum()
1189                   Ke[l * 2 + 0][ll * 2 + 1] += pp * jw * dbeta2 * u * v * ppl;  in THIJacobianLocal_Momentum()
1190                   Ke[l * 2 + 1][ll * 2 + 0] += pp * jw * dbeta2 * v * u * ppl;  in THIJacobianLocal_Momentum()
1191 …               Ke[l * 2 + 1][ll * 2 + 1] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * v * v * ppl;  in THIJacobianLocal_Momentum()
1204               Ke[ll * 2 + 0][l * 2 + 0] = Ke[l * 2 + 0][ll * 2 + 0];  in THIJacobianLocal_Momentum()
1205               Ke[ll * 2 + 1][l * 2 + 0] = Ke[l * 2 + 0][ll * 2 + 1];  in THIJacobianLocal_Momentum()
1206               Ke[ll * 2 + 0][l * 2 + 1] = Ke[l * 2 + 1][ll * 2 + 0];  in THIJacobianLocal_Momentum()
1207               Ke[ll * 2 + 1][l * 2 + 1] = Ke[l * 2 + 1][ll * 2 + 1];  in THIJacobianLocal_Momentum()
1211 …PetscCall(MatSetValuesBlockedLocal(B, 8, rc3blocked, 8, rc3blocked, &Ke[0][0], ADD_VALUES)); /* ve…  in THIJacobianLocal_Momentum()