7848   PetscInt        newNumPoints = 0, newNumIndices = 0;  in DMPlexAnchorsGetSubMatModification()  local
7902           newNumIndices += aDof;  in DMPlexAnchorsGetSubMatModification()
7913         newNumIndices += bSecDof;  in DMPlexAnchorsGetSubMatModification()
7933   if (outNumIndices) *outNumIndices = newNumIndices;  in DMPlexAnchorsGetSubMatModification()
7946 …newNumIndices, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Invalid size for closure %" PetscInt_FMT " should…  in DMPlexAnchorsGetSubMatModification()
7990     PetscCall(DMGetWorkArray(dm, newNumIndices, MPIU_INT, &newIndices));  in DMPlexAnchorsGetSubMatModification()
7991     PetscCall(DMGetWorkArray(dm, newNumIndices, MPIU_INT, &tmpNewIndices));  in DMPlexAnchorsGetSubMatModification()
7994     for (PetscInt i = 0; i < newNumIndices; i++) newIndices[i] = -1;  in DMPlexAnchorsGetSubMatModification()
7999     PetscCall(DMGetWorkArray(dm, numIndices * newNumIndices, MPIU_SCALAR, &modMat));  in DMPlexAnchorsGetSubMatModification()
8000     PetscCall(DMGetWorkArray(dm, numIndices * newNumIndices, MPIU_SCALAR, &tmpMat));  in DMPlexAnchorsGetSubMatModification()
8001     PetscCall(PetscArrayzero(modMat, newNumIndices * numIndices));  in DMPlexAnchorsGetSubMatModification()
8034 …for (PetscInt e = 0; e < aDof; e++) modMat[(o + d) * newNumIndices + oNew + e] = tmpMat[d * aDof +…  in DMPlexAnchorsGetSubMatModification()
8041           for (PetscInt d = 0; d < bSecDof; d++) modMat[(o + d) * newNumIndices + oNew + d] = 1;  in DMPlexAnchorsGetSubMatModification()
8050     PetscCall(DMRestoreWorkArray(dm, numIndices * newNumIndices, MPIU_SCALAR, &tmpMat));  in DMPlexAnchorsGetSubMatModification()
8051     PetscCall(DMRestoreWorkArray(dm, newNumIndices, MPIU_INT, &tmpNewIndices));  in DMPlexAnchorsGetSubMatModification()
8052     PetscCall(DMRestoreWorkArray(dm, newNumIndices, MPIU_INT, &newIndices));  in DMPlexAnchorsGetSubMatModification()
8071   PetscInt     newNumIndices = -1;  in DMPlexAnchorsModifyMat_Internal()  local
8076 …ion(dm, section, numPoints, numIndices, points, perms, outNumPoints, &newNumIndices, outPoints, of…  in DMPlexAnchorsModifyMat_Internal()
8077   if (outNumIndices) *outNumIndices = newNumIndices;  in DMPlexAnchorsModifyMat_Internal()
8088       PetscCall(PetscBLASIntCast(newNumIndices, &M));  in DMPlexAnchorsModifyMat_Internal()
8091       PetscCall(DMGetWorkArray(dm, numRows * newNumIndices, MPIU_SCALAR, &newNewValues));  in DMPlexAnchorsModifyMat_Internal()
8094       numCols   = newNumIndices;  in DMPlexAnchorsModifyMat_Internal()
8106       PetscCall(PetscBLASIntCast(newNumIndices, &N));  in DMPlexAnchorsModifyMat_Internal()
8108       PetscCall(DMGetWorkArray(dm, newNumIndices * numCols, MPIU_SCALAR, &newNewValues));  in DMPlexAnchorsModifyMat_Internal()
8111 …if (newValues != values) PetscCall(DMRestoreWorkArray(dm, numIndices * newNumIndices, MPIU_SCALAR,…  in DMPlexAnchorsModifyMat_Internal()
8115     PetscCall(DMRestoreWorkArray(dm, numIndices * newNumIndices, MPIU_SCALAR, &modMat));  in DMPlexAnchorsModifyMat_Internal()