204 … PetscReal hx, PetscReal hy, const PetscReal dz[PETSC_RESTRICT], PetscReal phi[PETSC_RESTRICT], Pe…  in HexComputeGeometry()
221     phi[i]                 = HexQInterp[q][i];  in HexComputeGeometry()
710 static void PointwiseNonlinearity(THI thi, const Node n[PETSC_RESTRICT], const PetscReal phi[PETSC_…  in PointwiseNonlinearity()
720     *u += phi[l] * n[l].u;  in PointwiseNonlinearity()
721     *v += phi[l] * n[l].v;  in PointwiseNonlinearity()
731 static void PointwiseNonlinearity2D(THI thi, Node n[], PetscReal phi[], PetscReal dphi[4][2], Petsc…  in PointwiseNonlinearity2D()  argument
741     *u += phi[l] * n[l].u;  in PointwiseNonlinearity2D()
742     *v += phi[l] * 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()
827               const PetscReal jw = 0.25 * hx * hy / thi->rhog, *phi = QuadQInterp[q];  in THIFunctionLocal()  local
831                 u += phi[l] * n[l].u;  in THIFunctionLocal()
832                 v += phi[l] * n[l].v;  in THIFunctionLocal()
833                 rbeta2 += phi[l] * pn[l].beta2;  in THIFunctionLocal()
838                 const PetscReal pp = phi[l];  in THIFunctionLocal()
991         PetscReal   phi[4], dphi[4][2], jw, eta, deta, beta2, dbeta2;  in THIJacobianLocal_2D()  local
994           phi[l]     = QuadQInterp[q][l];  in THIJacobianLocal_2D()
997           h += phi[l] * pn[l].h;  in THIJacobianLocal_2D()
998           rbeta2 += phi[l] * pn[l].beta2;  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);
1149               const PetscReal jw = 0.25 * hx * hy / thi->rhog, *phi = QuadQInterp[q];  local
1153                 u += phi[l] * n[l].u;
1154                 v += phi[l] * n[l].v;
1155                 rbeta2 += phi[l] * pn[l].beta2;
1159                 const PetscReal pp = phi[l];
1161                   const PetscReal ppl = phi[ll];