276 static void PrmHexGetZ(const PrmNode pn[], PetscInt k, PetscInt zm, PetscReal zn[])  in PrmHexGetZ()  argument
278 …const PetscScalar zm1 = zm - 1, znl[8] = {pn[0].b + pn[0].h * (PetscScalar)k / zm1,       pn[1].b …  in PrmHexGetZ()
688   PetscInt  i, j, k, xs, xm, ys, ym, zs, zm, mx, my;  in THIInitial()  local
700   PetscCall(DMDAGetCorners(da3, &zs, &ys, &xs, &zm, &ym, &xm));  in THIInitial()
710       for (k = zs; k < zs + zm; k++) {  in THIInitial()
711 …const PetscScalar zm1 = zm - 1, drivingx = thi->rhog * (prm[i + 1][j].b + prm[i + 1][j].h - prm[i …  in THIInitial()
752   PetscInt  xs, ys, xm, ym, zm, i, j, k, q, l;  in THIFunctionLocal_3D()  local
758   xm = info->zm;  in THIFunctionLocal_3D()
760   zm = info->xm;  in THIFunctionLocal_3D()
774       for (k = 0; k < zm - 1; k++) {  in THIFunctionLocal_3D()
779         PrmHexGetZ(pn, k, zm, zn);  in THIFunctionLocal_3D()
825             const PetscReal   hz    = PetscRealPart(pn[0].h) / (zm - 1.);  in THIFunctionLocal_3D()
860   PetscInt xs, ys, xm, ym, zm, i, j, k;  in THIFunctionLocal_2D()  local
865   xm = info->zm;  in THIFunctionLocal_2D()
867   zm = info->xm;  in THIFunctionLocal_2D()
873       for (k = 0; k < zm; k++) {  in THIFunctionLocal_2D()
874         PetscScalar weight = (k == 0 || k == zm - 1) ? 0.5 / (zm - 1) : 1.0 / (zm - 1);  in THIFunctionLocal_2D()
991   PetscInt    i, j, xs, ys, zs, xm, ym, zm, mx, my, mz;  in THISurfaceStatistics()  local
1000   PetscCall(DMDAGetCorners(da3, &zs, &ys, &xs, &zm, &ym, &xm));  in THISurfaceStatistics()
1001   PetscCheck(zs == 0 && zm == mz, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Unexpected decomposition");  in THISurfaceStatistics()
1005       PetscReal u = PetscRealPart(x[i][j][zm - 1].u);  in THISurfaceStatistics()
1096   PetscInt  xs, ys, xm, ym, zm, i, j, k, q, l, ll;  in THIJacobianLocal_Momentum()  local
1102   xm = info->zm;  in THIJacobianLocal_Momentum()
1104   zm = info->xm;  in THIJacobianLocal_Momentum()
1113       for (k = 0; k < zm - 1; k++) {  in THIJacobianLocal_Momentum()
1119         PrmHexGetZ(pn, k, zm, zn);  in THIJacobianLocal_Momentum()
1167             const PetscReal   hz    = PetscRealPart(pn[0].h) / (zm - 1);  in THIJacobianLocal_Momentum()
1229   PetscInt xs, ys, xm, ym, zm, i, j, k;  in THIJacobianLocal_2D()  local
1234   xm = info->zm;  in THIJacobianLocal_2D()
1236   zm = info->xm;  in THIJacobianLocal_2D()
1238 …PetscCheck(zm <= 1024, ((PetscObject)info->da)->comm, PETSC_ERR_SUP, "Need to allocate more space"…  in THIJacobianLocal_2D()
1247       for (k = 0; k < zm; k++) { /* Coupling to velocity problem */  in THIJacobianLocal_2D()
1252 …zm - 1) ? 0.5 : 0.25, hW = w * (x2[i - 1][j].h + x2[i][j].h) / (zm - 1.), hE = w * (x2[i][j].h + x…  in THIJacobianLocal_2D()
1253                           hN = w * (x2[i][j].h + x2[i][j + 1].h) / (zm - 1.);  in THIJacobianLocal_2D()
1379       PetscInt i, j, k, f, xs, xm, ys, ym, zs, zm;  in THIDAVecView_VTK_XML()  local
1388       zm = range[3];  in THIDAVecView_VTK_XML()
1391       PetscCheck(xm * ym * zm * dof <= nmax, PETSC_COMM_SELF, PETSC_ERR_PLIB, "should not happen");  in THIDAVecView_VTK_XML()
1395         PetscCheck(nn == xm * ym * zm * dof, PETSC_COMM_SELF, PETSC_ERR_PLIB, "corrupt da3 send");  in THIDAVecView_VTK_XML()
1405 …MT " %" PetscInt_FMT " %" PetscInt_FMT " %" PetscInt_FMT "\">\n", zs, zs + zm - 1, ys, ys + ym - 1…  in THIDAVecView_VTK_XML()
1415           for (k = zs; k < zs + zm; k++) {  in THIDAVecView_VTK_XML()