4951   PetscInt *offsets, **closures;  in DMPlexGetFullJoin()  local
4963   PetscCall(PetscCalloc1(numPoints, &closures));  in DMPlexGetFullJoin()
4973     PetscCall(DMPlexGetTransitiveClosure(dm, points[p], PETSC_FALSE, &closureSize, &closures[p]));  in DMPlexGetFullJoin()
4981         if ((pStart > closures[p][i * 2]) || (pEnd <= closures[p][i * 2])) {  in DMPlexGetFullJoin()
4995 …for (joinSize = 0; joinSize < dof; ++joinSize) join[i][joinSize] = closures[0][(offsets[d] + joinS…  in DMPlexGetFullJoin()
5002         const PetscInt point = closures[p][(offsets[p * (depth + 2) + d] + c) * 2];  in DMPlexGetFullJoin()
5018 …ts; ++p) PetscCall(DMPlexRestoreTransitiveClosure(dm, points[p], PETSC_FALSE, NULL, &closures[p]));  in DMPlexGetFullJoin()
5019   PetscCall(PetscFree(closures));  in DMPlexGetFullJoin()
5153   PetscInt *offsets, **closures;  in DMPlexGetFullMeet()  local
5165   PetscCall(PetscMalloc1(numPoints, &closures));  in DMPlexGetFullMeet()
5175     PetscCall(DMPlexGetTransitiveClosure(dm, points[p], PETSC_TRUE, &closureSize, &closures[p]));  in DMPlexGetFullMeet()
5183         if ((pStart > closures[p][i * 2]) || (pEnd <= closures[p][i * 2])) {  in DMPlexGetFullMeet()
5197 …for (meetSize = 0; meetSize < dof; ++meetSize) meet[i][meetSize] = closures[0][(offsets[h] + meetS…  in DMPlexGetFullMeet()
5204         const PetscInt point = closures[p][(offsets[p * (height + 2) + h] + c) * 2];  in DMPlexGetFullMeet()
5220 …nts; ++p) PetscCall(DMPlexRestoreTransitiveClosure(dm, points[p], PETSC_TRUE, NULL, &closures[p]));  in DMPlexGetFullMeet()
5221   PetscCall(PetscFree(closures));  in DMPlexGetFullMeet()