40   PetscScalar *closure;  in LoadData2D()  local
43   PetscCall(PetscMalloc1(clSize, &closure));  in LoadData2D()
47       PetscCall(PetscArrayzero(closure, clSize));  in LoadData2D()
54 …for (c = 0; c < user->Nc[f]; ++c) closure[o++] = ((kj + joff) * (Ni * user->k[f] + 1) + ki + ioff)…  in LoadData2D()
58       PetscCall(DMPlexVecSetClosure(dm, NULL, u, j * Ni + i, closure, INSERT_VALUES));  in LoadData2D()
61   PetscCall(PetscFree(closure));  in LoadData2D()
68   PetscScalar *closure;  in LoadData3D()  local
71   PetscCall(PetscMalloc1(clSize, &closure));  in LoadData3D()
76         PetscCall(PetscArrayzero(closure, clSize));  in LoadData3D()
84 …for (c = 0; c < user->Nc[f]; ++c) closure[o++] = (((kk + koff) * (Nj * user->k[f] + 1) + kj + joff…  in LoadData3D()
89         PetscCall(DMPlexVecSetClosure(dm, NULL, u, (k * Nj + j) * Ni + i, closure, INSERT_VALUES));  in LoadData3D()
93   PetscCall(PetscFree(closure));  in LoadData3D()
126     PetscScalar *closure = NULL;  in ReadData2D()  local
129     PetscCall(DMPlexVecGetClosure(dm, NULL, u, cell, &closureSize, &closure));  in ReadData2D()
138 …PetscCall(PetscPrintf(PETSC_COMM_SELF, "%2.0f", (double)PetscRealPart(closure[(kj * (user->k[f] + …  in ReadData2D()
146     PetscCall(DMPlexVecRestoreClosure(dm, NULL, u, cell, &closureSize, &closure));  in ReadData2D()
159     PetscScalar *closure = NULL;  in ReadData3D()  local
162     PetscCall(DMPlexVecGetClosure(dm, NULL, u, cell, &closureSize, &closure));  in ReadData3D()
172 …PetscCall(PetscPrintf(PETSC_COMM_SELF, "%2.0f", (double)PetscRealPart(closure[((kk * (user->k[f] +…  in ReadData3D()
182     PetscCall(DMPlexVecRestoreClosure(dm, NULL, u, cell, &closureSize, &closure));  in ReadData3D()