3183   PetscInt     dim, Nf, field, totDim, cStart, cEnd, cell, ccell, s;  in DMPlexComputeInterpolatorGeneral()  local
3269           for (ccell = 0; ccell < numCoarseCells; ++ccell) {  in DMPlexComputeInterpolatorGeneral()
3273 …PetscCall(DMPlexGetClosureIndices(dmc, csection, globalCSection, coarseCells[ccell].index, PETSC_F…  in DMPlexComputeInterpolatorGeneral()
3279 …PetscCall(DMPlexComputeCellGeometryFEM(dmc, coarseCells[ccell].index, NULL, v0c, Jc, invJc, &detJc…  in DMPlexComputeInterpolatorGeneral()
3280               for (d = 0; d < dim; ++d) pVReal[d] = PetscRealPart(pV[ccell * dim + d]);  in DMPlexComputeInterpolatorGeneral()
3289 …            for (c = 0; c < Nc; ++c) elemMat[j] += T->T[0][j * Nc + c] * qweights[ccell * qNc + c];  in DMPlexComputeInterpolatorGeneral()
3293                   for (c = 0; c < Nc; ++c) elemMat[j] += 1.0 * qweights[ccell * qNc + c];  in DMPlexComputeInterpolatorGeneral()
3301 …PetscCall(DMPlexRestoreClosureIndices(dmc, csection, globalCSection, coarseCells[ccell].index, PET…  in DMPlexComputeInterpolatorGeneral()
3356   PetscInt     dim, Nf, field, totDim, cStart, cEnd, cell, ccell;  in DMPlexComputeMassMatrixGeneral()  local
3432             for (ccell = 0; ccell < numCoarseCells; ++ccell) {  in DMPlexComputeMassMatrixGeneral()
3433 …PetscCall(DMPlexGetClosureIndices(dmc, csection, globalCSection, coarseCells[ccell].index, PETSC_F…  in DMPlexComputeMassMatrixGeneral()
3443 …ll(DMPlexRestoreClosureIndices(dmc, csection, globalCSection, coarseCells[ccell].index, PETSC_FALS…  in DMPlexComputeMassMatrixGeneral()
3505       for (ccell = 0; ccell < numCoarseCells; ++ccell) {  in DMPlexComputeMassMatrixGeneral()
3509 …PetscCall(DMPlexGetClosureIndices(dmc, csection, globalCSection, coarseCells[ccell].index, PETSC_F…  in DMPlexComputeMassMatrixGeneral()
3511 …PetscCall(DMPlexComputeCellGeometryFEM(dmc, coarseCells[ccell].index, NULL, v0c, Jc, invJc, &detJc…  in DMPlexComputeMassMatrixGeneral()
3512         for (d = 0; d < dim; ++d) pVReal[d] = PetscRealPart(pV[ccell * dim + d]);  in DMPlexComputeMassMatrixGeneral()
3525 …) elemMat[j] += T->T[0][j * Nc + c] * Tfine->T[0][(ccell * numFIndices + i) * Nc + c] * qweights[c…  in DMPlexComputeMassMatrixGeneral()
3537               for (c = 0; c < Nc; ++c) elemMat[j] += 1.0 * 1.0 * qweights[ccell * Nc + c] * detJ;  in DMPlexComputeMassMatrixGeneral()
3541 …scInt_FMT " %" PetscInt_FMT " Nc: %" PetscInt_FMT " %" PetscInt_FMT "\n", ccell, Nq, i, numFIndice…  in DMPlexComputeMassMatrixGeneral()
3546 …PetscCall(DMPlexRestoreClosureIndices(dmc, csection, globalCSection, coarseCells[ccell].index, PET…  in DMPlexComputeMassMatrixGeneral()