192 …tscReal dz[restrict], PetscReal phi[restrict], PetscReal dphi[restrict][3], PetscReal *restrict jw)  in HexComputeGeometry()  argument
210   *jw = 1.0 * jdet;  in HexComputeGeometry()
789           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  in THIFunctionLocal_3D()  local
796           HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw);  in THIFunctionLocal_3D()
798           jw /= thi->rhog; /* scales residuals to be O(1) */  in THIFunctionLocal_3D()
804 …fn[l]->u += dp[0] * jw * eta * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp…  in THIFunctionLocal_3D()
805 …fn[l]->v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp…  in THIFunctionLocal_3D()
806             fn[l]->u += pp * jw * udot * thi->inertia * pp;  in THIFunctionLocal_3D()
807             fn[l]->v += pp * jw * vdot * thi->inertia * pp;  in THIFunctionLocal_3D()
831               const PetscReal jw = 0.25 * hx * hy, *phi = QuadQInterp[q];  in THIFunctionLocal_3D()  local
843                 fn[ls + l]->u += pp * jw * beta2 * u;  in THIFunctionLocal_3D()
844                 fn[ls + l]->v += pp * jw * beta2 * v;  in THIFunctionLocal_3D()
1128           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  in THIJacobianLocal_Momentum()  local
1131           HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw);  in THIJacobianLocal_Momentum()
1133           jw /= thi->rhog; /* residuals are scaled by this factor */  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 …][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu *…  in THIJacobianLocal_Momentum()
1149 …][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv *…  in THIJacobianLocal_Momentum()
1150 …][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu *…  in THIJacobianLocal_Momentum()
1151 …][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv *…  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()
1158 …          Kcpl[FieldIndex(Node, l, u)][FieldIndex(PrmNode, ll, h)] += pp * jw * thi->rhog * dpl[0];  in THIJacobianLocal_Momentum()
1159 …          Kcpl[FieldIndex(Node, l, u)][FieldIndex(PrmNode, ll, b)] += pp * jw * thi->rhog * dpl[0];  in THIJacobianLocal_Momentum()
1160 …          Kcpl[FieldIndex(Node, l, v)][FieldIndex(PrmNode, ll, h)] += pp * jw * thi->rhog * dpl[1];  in THIJacobianLocal_Momentum()
1161 …          Kcpl[FieldIndex(Node, l, v)][FieldIndex(PrmNode, ll, b)] += pp * jw * thi->rhog * dpl[1];  in THIJacobianLocal_Momentum()
1175               const PetscReal jw = 0.25 * hx * hy, *phi = QuadQInterp[q];  in THIJacobianLocal_Momentum()  local
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()