729 …strict u, PetscScalar *restrict v, PetscScalar du[restrict 3], PetscScalar dv[restrict 3], PetscRe…  in PointwiseNonlinearity()
735   dv[0] = dv[1] = dv[2] = 0;  in PointwiseNonlinearity()
743       dv[ll] += dphi[l][ll] * n[l].v;  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 …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 …v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp[2] * j…  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 … 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu * (du[1] + …  in THIJacobianLocal_Momentum()
1149 … 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv * (du[1] + …  in THIJacobianLocal_Momentum()
1150 … * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu * (du[1] + dv[0]) + dp[…  in THIJacobianLocal_Momentum()
1151 … * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv * (du[1] + dv[0]) + dp[…  in THIJacobianLocal_Momentum()