1570     PetscInt       c, f, fx, fy, fz, e, ex, ey, ez;  in DMPlexCreateCubeMesh_Internal()  local
1579       for (fy = 0; fy < numYEdges; ++fy) {  in DMPlexCreateCubeMesh_Internal()
1581           PetscInt cell  = (fz * numYEdges + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1582           PetscInt faceB = firstZFace + (fy * numXEdges + fx) * numZVertices + fz;  in DMPlexCreateCubeMesh_Internal()
1583 …    PetscInt faceT = firstZFace + (fy * numXEdges + fx) * numZVertices + ((fz + 1) % numZVertices);  in DMPlexCreateCubeMesh_Internal()
1584           PetscInt faceF = firstYFace + (fz * numXEdges + fx) * numYVertices + fy;  in DMPlexCreateCubeMesh_Internal()
1585 …    PetscInt faceK = firstYFace + (fz * numXEdges + fx) * numYVertices + ((fy + 1) % numYVertices);  in DMPlexCreateCubeMesh_Internal()
1586           PetscInt faceL = firstXFace + (fz * numYEdges + fy) * numXVertices + fx;  in DMPlexCreateCubeMesh_Internal()
1587 …    PetscInt faceR = firstXFace + (fz * numYEdges + fy) * numXVertices + ((fx + 1) % numXVertices);  in DMPlexCreateCubeMesh_Internal()
1602 …if (bdY != DM_BOUNDARY_NONE && fy == numYEdges - 1 && cutLabel) PetscCall(DMLabelSetValue(cutLabel…  in DMPlexCreateCubeMesh_Internal()
1609       for (fy = 0; fy < numYEdges; ++fy) {  in DMPlexCreateCubeMesh_Internal()
1611           PetscInt face    = firstXFace + (fz * numYEdges + fy) * numXVertices + fx;  in DMPlexCreateCubeMesh_Internal()
1612           PetscInt edgeL   = firstZEdge + (fy * numXVertices + fx) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1613 …  PetscInt edgeR   = firstZEdge + (((fy + 1) % numYVertices) * numXVertices + fx) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1614           PetscInt edgeB   = firstYEdge + (fz * numXVertices + fx) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1615 …  PetscInt edgeT   = firstYEdge + (((fz + 1) % numZVertices) * numXVertices + fx) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1643         for (fy = 0; fy < numYVertices; ++fy) {  in DMPlexCreateCubeMesh_Internal()
1644           PetscInt face    = firstYFace + (fz * numXEdges + fx) * numYVertices + fy;  in DMPlexCreateCubeMesh_Internal()
1645           PetscInt edgeL   = firstZEdge + (fy * numXVertices + fx) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1646 …  PetscInt edgeR   = firstZEdge + (fy * numXVertices + ((fx + 1) % numXVertices)) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1647           PetscInt edgeB   = firstXEdge + (fz * numYVertices + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1648 …  PetscInt edgeT   = firstXEdge + (((fz + 1) % numZVertices) * numYVertices + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1655               if (fy == numYVertices - 1) {  in DMPlexCreateCubeMesh_Internal()
1658               } else if (fy == 0) {  in DMPlexCreateCubeMesh_Internal()
1674     for (fy = 0; fy < numYEdges; ++fy) {  in DMPlexCreateCubeMesh_Internal()
1677           PetscInt face    = firstZFace + (fy * numXEdges + fx) * numZVertices + fz;  in DMPlexCreateCubeMesh_Internal()
1678           PetscInt edgeL   = firstYEdge + (fz * numXVertices + fx) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1679 …  PetscInt edgeR   = firstYEdge + (fz * numXVertices + ((fx + 1) % numXVertices)) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1680           PetscInt edgeB   = firstXEdge + (fz * numYVertices + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1681 …  PetscInt edgeT   = firstXEdge + (fz * numYVertices + ((fy + 1) % numYVertices)) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1687               edgeR += numYEdges - 1 - 2 * fy;  in DMPlexCreateCubeMesh_Internal()
1690             if (bdY == DM_BOUNDARY_TWIST && fy == numYEdges - 1) {  in DMPlexCreateCubeMesh_Internal()
1695 …if (bdY != DM_BOUNDARY_NONE && fy == numYEdges - 1 && cutLabel) PetscCall(DMLabelSetValue(cutLabel…  in DMPlexCreateCubeMesh_Internal()