3037       PetscInt     NcJ = 0, cpdim = 0, j, qNc;  in DMPlexComputeInterpolatorNested()  local
3047         PetscCall(PetscFEGetDimension(fe, &cpdim));  in DMPlexComputeInterpolatorNested()
3057               for (j = 0; j < cpdim; ++j) {  in DMPlexComputeInterpolatorNested()
3066 … < Nc; ++c) elemMat[(offsetI + i) * cTotDim + offsetJ + j] += T->T[0][k * cpdim * NcJ + j * Nc + c…  in DMPlexComputeInterpolatorNested()
3077         cpdim = 1;  in DMPlexComputeInterpolatorNested()
3086               for (j = 0; j < cpdim; ++j) {  in DMPlexComputeInterpolatorNested()
3093       offsetJ += cpdim;  in DMPlexComputeInterpolatorNested()
3245           PetscInt           numCoarseCells, cpdim, row = findices[i + off], q, c, j;  in DMPlexComputeInterpolatorGeneral()  local
3274             if (id == PETSCFE_CLASSID) PetscCall(PetscFEGetDimension((PetscFE)obj, &cpdim));  in DMPlexComputeInterpolatorGeneral()
3275             else cpdim = 1;  in DMPlexComputeInterpolatorGeneral()
3286                 PetscCall(PetscArrayzero(elemMat, cpdim));  in DMPlexComputeInterpolatorGeneral()
3288                 for (j = 0; j < cpdim; ++j) {  in DMPlexComputeInterpolatorGeneral()
3292                 for (j = 0; j < cpdim; ++j) {  in DMPlexComputeInterpolatorGeneral()
3300             PetscCall(MatSetValues(interp, 1, &row, cpdim, &cindices[off], elemMat, INSERT_VALUES));  in DMPlexComputeInterpolatorGeneral()
3481       PetscInt           numCoarseCells, cpdim, q, c, j;  in DMPlexComputeMassMatrixGeneral()  local
3519           PetscCall(PetscFEGetDimension(fe, &cpdim));  in DMPlexComputeMassMatrixGeneral()
3523             PetscCall(PetscArrayzero(elemMat, cpdim));  in DMPlexComputeMassMatrixGeneral()
3524             for (j = 0; j < cpdim; ++j) {  in DMPlexComputeMassMatrixGeneral()
3529 …Indices == cpdim, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Number of element matrix columns %" PetscInt_F…  in DMPlexComputeMassMatrixGeneral()
3533           cpdim = 1;  in DMPlexComputeMassMatrixGeneral()
3535             PetscCall(PetscArrayzero(elemMat, cpdim));  in DMPlexComputeMassMatrixGeneral()
3536             for (j = 0; j < cpdim; ++j) {  in DMPlexComputeMassMatrixGeneral()
3542 …Indices == cpdim, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Number of element matrix columns %" PetscInt_F…  in DMPlexComputeMassMatrixGeneral()
3637     PetscInt       order = -1, NcF, NcC, fpdim, cpdim;  in DMPlexComputeInjectorFEM()  local
3647       PetscCall(PetscDualSpaceGetDimension(QC, &cpdim));  in DMPlexComputeInjectorFEM()
3655       PetscCall(PetscDualSpaceGetDimension(QC, &cpdim));  in DMPlexComputeInjectorFEM()
3658     for (c = 0; c < cpdim; ++c) {  in DMPlexComputeInjectorFEM()
3692     offsetC += cpdim;  in DMPlexComputeInjectorFEM()