2942   PetscInt     cTotDim = 0, rTotDim = 0;  in DMPlexComputeInterpolatorNested()  local
3001     cTotDim += cNb;  in DMPlexComputeInterpolatorNested()
3003   PetscCall(PetscMalloc1(rTotDim * cTotDim, &elemMat));  in DMPlexComputeInterpolatorNested()
3004   PetscCall(PetscArrayzero(elemMat, rTotDim * cTotDim));  in DMPlexComputeInterpolatorNested()
3066 …for (c = 0; c < Nc; ++c) elemMat[(offsetI + i) * cTotDim + offsetJ + j] += T->T[0][k * cpdim * NcJ…  in DMPlexComputeInterpolatorNested()
3087 …for (c = 0; c < Nc; ++c) elemMat[(offsetI + i) * cTotDim + offsetJ + j] += 1.0 * qweights[p * qNc …  in DMPlexComputeInterpolatorNested()
3098   if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(0, name, rTotDim, cTotDim, elemMat));  in DMPlexComputeInterpolatorNested()
3111 … PetscCall(PetscCalloc3(rTotDim * cTotDim, &vals, cTotDim, &cellCIndices, rTotDim, &cellFIndices));  in DMPlexComputeInterpolatorNested()
3116 …PetscCall(MatSetValues(preallocator, rTotDim, cellFIndices, cTotDim, cellCIndices, vals, INSERT_VA…  in DMPlexComputeInterpolatorNested()
3587 …PetscInt     cTotDim, fTotDim = 0, Nf, f, field, cStart, cEnd, c, dim, d, startC, endC, offsetC, o…  in DMPlexComputeInjectorFEM()  local
3631   PetscCall(PetscDSGetTotalDimension(prob, &cTotDim));  in DMPlexComputeInjectorFEM()
3632   PetscCall(PetscMalloc1(cTotDim, &cmap));  in DMPlexComputeInjectorFEM()
3705   PetscCall(PetscMalloc2(cTotDim, &cellCIndices, fTotDim, &cellFIndices));  in DMPlexComputeInjectorFEM()
3711     for (d = 0; d < cTotDim; ++d) {  in DMPlexComputeInjectorFEM()