impls/plex/plexgeometry.c

982   PetscGridHash   lbox;  in DMPlexComputeGridHash_Internal()  local
992   PetscCall(DMPlexCreateGridHash(dm, &lbox));  in DMPlexComputeGridHash_Internal()
1002     PetscCall(PetscGridHashSetGrid(lbox, n, NULL));  in DMPlexComputeGridHash_Internal()
1004 … %" PetscInt_FMT "\n  h %g %g %g\n", (double)lbox->lower[0], (double)lbox->lower[1], cdim > 2 ? (d…  in DMPlexComputeGridHash_Internal()
1005 …lbox->upper[0], (double)lbox->upper[1], cdim > 2 ? (double)lbox->upper[2] : 0, n[0], n[1], cdim > …  in DMPlexComputeGridHash_Internal()
1013   PetscCall(DMLabelCreate(PETSC_COMM_SELF, "cells", &lbox->cellsSparse));  in DMPlexComputeGridHash_Internal()
1014   PetscCall(DMLabelCreateIndex(lbox->cellsSparse, cStart, cEnd));  in DMPlexComputeGridHash_Internal()
1019     const PetscReal   *h                 = lbox->h;  in DMPlexComputeGridHash_Internal()
1031     PetscCall(PetscGridHashGetEnclosingBox(lbox, numCoords / cdim, coords, dboxes, boxes));  in DMPlexComputeGridHash_Internal()
1046       lp[d] = lbox->lower[d] + dlim[d * 2 + 0] * h[d];  in DMPlexComputeGridHash_Internal()
1072           const PetscInt box    = (k * lbox->n[1] + j) * lbox->n[0] + i;  in DMPlexComputeGridHash_Internal()
1109             PetscCall(DMLabelSetValue(lbox->cellsSparse, c, box));  in DMPlexComputeGridHash_Internal()
1134                 PetscCall(DMLabelSetValue(lbox->cellsSparse, c, box));  in DMPlexComputeGridHash_Internal()
1149                 PetscCall(DMLabelSetValue(lbox->cellsSparse, c, box));  in DMPlexComputeGridHash_Internal()
1166               PetscCall(DMLabelSetValue(lbox->cellsSparse, c, box));  in DMPlexComputeGridHash_Internal()
1186               PetscCall(DMLabelSetValue(lbox->cellsSparse, c, box));  in DMPlexComputeGridHash_Internal()
1208               PetscCall(DMLabelSetValue(lbox->cellsSparse, c, box));  in DMPlexComputeGridHash_Internal()
1219         lp[0]             = lbox->lower[0] + dlim[0 * 2 + 0] * h[0];  in DMPlexComputeGridHash_Internal()
1227       lp[1]             = lbox->lower[1] + dlim[1 * 2 + 0] * h[1];  in DMPlexComputeGridHash_Internal()
1238   if (debug) PetscCall(DMLabelView(lbox->cellsSparse, PETSC_VIEWER_STDOUT_SELF));  in DMPlexComputeGridHash_Internal()
1239   PetscCall(DMLabelConvertToSection(lbox->cellsSparse, &lbox->cellSection, &lbox->cells));  in DMPlexComputeGridHash_Internal()
1240   PetscCall(DMLabelDestroy(&lbox->cellsSparse));  in DMPlexComputeGridHash_Internal()
1241   *localBox = lbox;  in DMPlexComputeGridHash_Internal()
1302     if (!mesh->lbox) {  in DMLocatePoints_Plex()
1304       PetscCall(DMPlexComputeGridHash_Internal(dm, &mesh->lbox));  in DMLocatePoints_Plex()
1310     PetscCall(ISGetIndices(mesh->lbox->cells, &boxCells));  in DMLocatePoints_Plex()
1358 …PetscCall(PetscGridHashGetEnclosingBoxQuery(mesh->lbox, mesh->lbox->cellSection, 1, point, dbin, &…  in DMLocatePoints_Plex()
1362         PetscCall(PetscSectionGetDof(mesh->lbox->cellSection, bin, &numCells));  in DMLocatePoints_Plex()
1363         PetscCall(PetscSectionGetOffset(mesh->lbox->cellSection, bin, &cellOffset));  in DMLocatePoints_Plex()
1397   if (hash) PetscCall(ISRestoreIndices(mesh->lbox->cells, &boxCells));  in DMLocatePoints_Plex()
1405         PetscCall(PetscGridHashGetEnclosingBox(mesh->lbox, 1, point, dbin, &bin));  in DMLocatePoints_Plex()
1406         PetscCall(PetscSectionGetDof(mesh->lbox->cellSection, bin, &numCells));  in DMLocatePoints_Plex()
1407         PetscCall(PetscSectionGetOffset(mesh->lbox->cellSection, bin, &cellOffset));  in DMLocatePoints_Plex()