285 …d4(const PetscReal dphi[][4][2], PetscReal hx, PetscReal hy, const PrmNode pn[4], PrmNode dp[4][2])  in QuadComputeGrad4()
289   PetscScalar (*restrict dpg)[2][PRMNODE_SIZE]   = (PetscScalar (*)[2][PRMNODE_SIZE])dp;  in QuadComputeGrad4()
803             const PetscReal   pp = phi[l], *dp = dphi[l];  in THIFunctionLocal_3D()  local
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()
1136             const PetscReal pp = phi[l], *restrict dp = dphi[l];  in THIJacobianLocal_Momentum()  local
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 … * 2 + 0][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta …  in THIJacobianLocal_Momentum()
1149 … * 2 + 0][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta …  in THIJacobianLocal_Momentum()
1150 … * 2 + 1][ll * 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta …  in THIJacobianLocal_Momentum()
1151 … * 2 + 1][ll * 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta …  in THIJacobianLocal_Momentum()