290 static void PrmHexGetZ(const PrmNode pn[], PetscInt k, PetscInt zm, PetscReal zn[])  in PrmHexGetZ()  argument
292 …pn[0].b + pn[0].h * (PetscScalar)k / zm1,       pn[1].b + pn[1].h * (PetscScalar)k / zm1,       pn…  in PrmHexGetZ()
293 …pn[0].b + pn[0].h * (PetscScalar)(k + 1) / zm1, pn[1].b + pn[1].h * (PetscScalar)(k + 1) / zm1, pn…  in PrmHexGetZ()
776       PrmNode pn[4];  in THIFunctionLocal()  local
777       QuadExtract(prm, i, j, pn);  in THIFunctionLocal()
782         PrmHexGetZ(pn, k, zm, zn);  in THIFunctionLocal()
821             const PetscReal   hz    = PetscRealPart(pn[0].h) / (zm - 1.);  in THIFunctionLocal()
833                 rbeta2 += phi[l] * pn[l].beta2;  in THIFunctionLocal()
985       PrmNode     pn[4];  in THIJacobianLocal_2D()  local
987       QuadExtract(prm, i, j, pn);  in THIJacobianLocal_2D()
997           h += phi[l] * pn[l].h;  in THIJacobianLocal_2D()
998           rbeta2 += phi[l] * pn[l].beta2;  in THIJacobianLocal_2D()
1064       PrmNode pn[4];  local
1065       QuadExtract(prm, i, j, pn);
1072         PrmHexGetZ(pn, k, zm, zn);
1143             const PetscReal   hz    = PetscRealPart(pn[0].h) / (zm - 1);
1155                 rbeta2 += phi[l] * pn[l].beta2;