7854   PetscScalar    *modMat         = NULL;  in DMPlexAnchorsGetSubMatModification()  local
7999     PetscCall(DMGetWorkArray(dm, numIndices * newNumIndices, MPIU_SCALAR, &modMat));  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()
8048     *outMat = modMat;  in DMPlexAnchorsGetSubMatModification()
8070   PetscScalar *modMat        = NULL;  in DMPlexAnchorsModifyMat_Internal()  local
8076 …ces, points, perms, outNumPoints, &newNumIndices, outPoints, offsets, outValues ? &modMat : NULL));  in DMPlexAnchorsModifyMat_Internal()
8078   if (modMat) {  in DMPlexAnchorsModifyMat_Internal()
8093 …PetscCallBLAS("BLASgemm", BLASgemm_("N", "N", &M, &N, &K, &a, modMat, &M, newValues, &K, &b, newNe…  in DMPlexAnchorsModifyMat_Internal()
8110 …PetscCallBLAS("BLASgemm", BLASgemm_("N", "T", &M, &N, &K, &a, newValues, &M, modMat, &N, &b, newNe…  in DMPlexAnchorsModifyMat_Internal()
8115     PetscCall(DMRestoreWorkArray(dm, numIndices * newNumIndices, MPIU_SCALAR, &modMat));  in DMPlexAnchorsModifyMat_Internal()