729 …trict 8][3], PetscScalar *restrict u, PetscScalar *restrict v, PetscScalar du[restrict 3], PetscSc…  in PointwiseNonlinearity()
734   du[0] = du[1] = du[2] = 0;  in PointwiseNonlinearity()
742       du[ll] += dphi[l][ll] * n[l].u;  in PointwiseNonlinearity()
746 …gam = Sqr(du[0]) + Sqr(dv[1]) + du[0] * dv[1] + 0.25 * Sqr(du[1] + dv[0]) + 0.25 * Sqr(du[2]) + 0.…  in PointwiseNonlinearity()
790           PetscScalar du[3], dv[3], u, v, udot = 0, vdot = 0;  in THIFunctionLocal_3D()  local
797           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIFunctionLocal_3D()
804 …n[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()
1129           PetscScalar du[3], dv[3], u, v;  in THIJacobianLocal_Momentum()  local
1132           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIJacobianLocal_Momentum()
1140 …dgdu = 2. * du[0] * dpl[0] + dv[1] * dpl[0] + 0.5 * (du[1] + dv[0]) * dpl[1] + 0.5 * du[2] * dpl[2…  in THIJacobianLocal_Momentum()
1141 …dgdv = 2. * dv[1] * dpl[1] + du[0] * dpl[1] + 0.5 * (du[1] + dv[0]) * dpl[0] + 0.5 * dv[2] * dpl[2…  in THIJacobianLocal_Momentum()
1148 … * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu * (du[1] + dv[0]) + dp[…  in THIJacobianLocal_Momentum()
1149 … * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv * (du[1] + dv[0]) + dp[…  in THIJacobianLocal_Momentum()
1150 … 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu * (du[1…  in THIJacobianLocal_Momentum()
1151 … 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv * (du[1…  in THIJacobianLocal_Momentum()