710 … PetscScalar *PETSC_RESTRICT u, PetscScalar *PETSC_RESTRICT v, PetscScalar du[PETSC_RESTRICT], Pet…  in PointwiseNonlinearity()
715   du[0] = du[1] = du[2] = 0;  in PointwiseNonlinearity()
723       du[ll] += dphi[l][ll] * n[l].u;  in PointwiseNonlinearity()
727 …gam = PetscSqr(du[0]) + PetscSqr(dv[1]) + du[0] * dv[1] + 0.25 * PetscSqr(du[1] + dv[0]) + 0.25 * …  in PointwiseNonlinearity()
731 …l phi[], PetscReal dphi[4][2], PetscScalar *u, PetscScalar *v, PetscScalar du[], PetscScalar dv[],…  in PointwiseNonlinearity2D()  argument
736   du[0] = du[1] = 0;  in PointwiseNonlinearity2D()
744       du[ll] += dphi[l][ll] * n[l].u;  in PointwiseNonlinearity2D()
748   gam = PetscSqr(du[0]) + PetscSqr(dv[1]) + du[0] * dv[1] + 0.25 * PetscSqr(du[1] + dv[0]);  in PointwiseNonlinearity2D()
792           PetscScalar du[3], dv[3], u, v;  in THIFunctionLocal()  local
795           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIFunctionLocal()
802 …n[l]->u += dp[0] * jw * eta * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp[…  in THIFunctionLocal()
803 …fn[l]->v += dp[1] * jw * eta * (2. * du[0] + 4. * dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp…  in THIFunctionLocal()
992         PetscScalar u, v, du[2], dv[2], h = 0, rbeta2 = 0;  in THIJacobianLocal_2D()  local
1001         PointwiseNonlinearity2D(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIJacobianLocal_2D()
1008             dgdu = 2. * du[0] * dpl[0] + dv[1] * dpl[0] + 0.5 * (du[1] + dv[0]) * dpl[1];  in THIJacobianLocal_2D()
1009             dgdv = 2. * dv[1] * dpl[1] + du[0] * dpl[1] + 0.5 * (du[1] + dv[0]) * dpl[0];  in THIJacobianLocal_2D()
1016 …][ll * 2 + 0] += dp[0] * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu *…  in THIJacobianLocal_2D()
1017 …][ll * 2 + 1] += dp[0] * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv *…  in THIJacobianLocal_2D()
1018 … 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu * (du[1…  in THIJacobianLocal_2D()
1019 … 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv * (du[1…  in THIJacobianLocal_2D()
1081           PetscScalar du[3], dv[3], u, v;  local
1084           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);
1091 …p42, p42, _MM_SHUFFLE2(0, 1)), du0 = _mm_set1_pd(du[0]), du1 = _mm_set1_pd(du[1]), du2 = _mm_set1_…
1110 …dgdu = 2. * du[0] * dpl[0] + dv[1] * dpl[0] + 0.5 * (du[1] + dv[0]) * dpl[1] + 0.5 * du[2] * dpl[2…
1111 …dgdv = 2. * dv[1] * dpl[1] + du[0] * dpl[1] + 0.5 * (du[1] + dv[0]) * dpl[0] + 0.5 * dv[2] * dpl[2…
1118 … * jw * deta * dgdu * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu * (du[1] + dv[0]) + dp[…
1119 … * jw * deta * dgdv * (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdv * (du[1] + dv[0]) + dp[…
1120 … 2 + 0] += dp[1] * jw * deta * dgdu * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdu * (du[1…
1121 … 2 + 1] += dp[1] * jw * deta * dgdv * (4. * dv[1] + 2. * du[0]) + dp[0] * jw * deta * dgdv * (du[1…