801             const PetscReal pp = phi[l], *dp = dphi[l];  in THIFunctionLocal()  local
802 …* dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp[2] * jw * eta * du[2] + pp * jw * thi->rhog * d…  in THIFunctionLocal()
803 …* dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp[2] * jw * eta * dv[2] + pp * jw * thi->rhog * d…  in THIFunctionLocal()
838                 const PetscReal pp = phi[l];  in THIFunctionLocal()  local
839                 fn[ls + l]->u += pp * jw * beta2 * u;  in THIFunctionLocal()
840                 fn[ls + l]->v += pp * jw * beta2 * v;  in THIFunctionLocal()
1004           const PetscReal pp = phi[l], *dp = dphi[l];  in THIJacobianLocal_2D()  local
1011 …l * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + pp * jw * (beta2 / h) *…  in THIJacobianLocal_2D()
1014 …l * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + pp * jw * (beta2 / h) *…  in THIJacobianLocal_2D()
1016 …* (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu * (du[1] + dv[0]) + pp * jw * (dbeta2 / h) …  in THIJacobianLocal_2D()
1017 …* (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv * (du[1] + dv[0]) + pp * jw * (dbeta2 / h) …  in THIJacobianLocal_2D()
1018 …* (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu * (du[1] + dv[0]) + pp * jw * (dbeta2 / h) …  in THIJacobianLocal_2D()
1019 …* (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv * (du[1] + dv[0]) + pp * jw * (dbeta2 / h) …  in THIJacobianLocal_2D()
1159                 const PetscReal pp = phi[l];  local
1162 …               Ke[l * 2 + 0][ll * 2 + 0] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * u * u * ppl;
1163                   Ke[l * 2 + 0][ll * 2 + 1] += pp * jw * dbeta2 * u * v * ppl;
1164                   Ke[l * 2 + 1][ll * 2 + 0] += pp * jw * dbeta2 * v * u * ppl;
1165 …               Ke[l * 2 + 1][ll * 2 + 1] += pp * jw * beta2 * ppl + pp * jw * dbeta2 * v * v * ppl;