193 static void HexGrad(const PetscReal dphi[][3], const PetscReal zn[], PetscReal dz[])  in HexGrad()
198     dz[0] += dphi[i][0] * zn[i];  in HexGrad()
199     dz[1] += dphi[i][1] * zn[i];  in HexGrad()
200     dz[2] += dphi[i][2] * zn[i];  in HexGrad()
204 …st PetscReal dz[PETSC_RESTRICT], PetscReal phi[PETSC_RESTRICT], PetscReal dphi[PETSC_RESTRICT][3],…  in HexComputeGeometry()
222     dphi[i][0]             = dphir[0] * ijac[0][0] + dphir[1] * ijac[1][0] + dphir[2] * ijac[2][0];  in HexComputeGeometry()
223     dphi[i][1]             = dphir[0] * ijac[0][1] + dphir[1] * ijac[1][1] + dphir[2] * ijac[2][1];  in HexComputeGeometry()
224     dphi[i][2]             = dphir[0] * ijac[0][2] + dphir[1] * ijac[1][2] + dphir[2] * ijac[2][2];  in HexComputeGeometry()
710 …st Node n[PETSC_RESTRICT], const PetscReal phi[PETSC_RESTRICT], PetscReal dphi[PETSC_RESTRICT][3],…  in PointwiseNonlinearity()
723       du[ll] += dphi[l][ll] * n[l].u;  in PointwiseNonlinearity()
724       dv[ll] += dphi[l][ll] * n[l].v;  in PointwiseNonlinearity()
731 static void PointwiseNonlinearity2D(THI thi, Node n[], PetscReal phi[], PetscReal dphi[4][2], Petsc…  in PointwiseNonlinearity2D()
744       du[ll] += dphi[l][ll] * n[l].u;  in PointwiseNonlinearity2D()
745       dv[ll] += dphi[l][ll] * n[l].v;  in PointwiseNonlinearity2D()
791           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  in THIFunctionLocal()  local
794           HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw);  in THIFunctionLocal()
795           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIFunctionLocal()
801             const PetscReal pp = phi[l], *dp = dphi[l];  in THIFunctionLocal()
991         PetscReal   phi[4], dphi[4][2], jw, eta, deta, beta2, dbeta2;  in THIJacobianLocal_2D()  local
995           dphi[l][0] = QuadQDeriv[q][l][0] * 2. / hx;  in THIJacobianLocal_2D()
996           dphi[l][1] = QuadQDeriv[q][l][1] * 2. / hy;  in THIJacobianLocal_2D()
1001         PointwiseNonlinearity2D(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);  in THIJacobianLocal_2D()
1004           const PetscReal pp = phi[l], *dp = dphi[l];  in THIJacobianLocal_2D()
1006             const PetscReal ppl = phi[ll], *dpl = dphi[ll];  in THIJacobianLocal_2D()
1080           PetscReal   dz[3], phi[8], dphi[8][3], jw, eta, deta;  local
1083           HexComputeGeometry(q, hx, hy, dz, phi, dphi, &jw);
1084           PointwiseNonlinearity(thi, n, phi, dphi, &u, &v, du, dv, &eta, &deta);
1088             const PetscReal *PETSC_RESTRICT dp = dphi[l];
1104               const PetscReal *PETSC_RESTRICT dpl = dphi[ll];