2144     PetscInt numColIndices = 0;  in DMPlexComputeInterpolatorTree()  local
2172 … localCoarse, closureSize, numRowIndices, closure, NULL, NULL, NULL, &numColIndices, NULL, NULL, n…  in DMPlexComputeInterpolatorTree()
2174 …if (!numColIndices) { /* there are no hanging constraint modifications, so the matrix is just the …  in DMPlexComputeInterpolatorTree()
2175         numColIndices = numRowIndices;  in DMPlexComputeInterpolatorTree()
2187         matSize = numRowIndices * numColIndices;  in DMPlexComputeInterpolatorTree()
2204           numColIndices += aLocalDof;  in DMPlexComputeInterpolatorTree()
2216           matSize = numColIndices * dof;  in DMPlexComputeInterpolatorTree()
2219         numColIndices = dof;  in DMPlexComputeInterpolatorTree()
2224 …PetscCall(PetscSectionSetDof(rootIndicesSec, p, numColIndices ? numColIndices + 2 * numFields : 0)…  in DMPlexComputeInterpolatorTree()
2236       PetscInt     numRowIndices = 0, numColIndices, matSize, dof;  in DMPlexComputeInterpolatorTree()  local
2242       PetscCall(PetscSectionGetDof(rootIndicesSec, p, &numColIndices));  in DMPlexComputeInterpolatorTree()
2243       if (!numColIndices) continue;  in DMPlexComputeInterpolatorTree()
2250       numColIndices -= 2 * numFields;  in DMPlexComputeInterpolatorTree()
2266 …PetscCheck(numIndices == numColIndices, PETSC_COMM_SELF, PETSC_ERR_PLIB, "mismatching constraint i…  in DMPlexComputeInterpolatorTree()
2267           for (i = 0; i < numColIndices; i++) pInd[i] = indices[i];  in DMPlexComputeInterpolatorTree()
2269             pInd[numColIndices + i]             = offsets[i + 1];  in DMPlexComputeInterpolatorTree()
2270             pInd[numColIndices + numFields + i] = offsets[i + 1];  in DMPlexComputeInterpolatorTree()
2328             for (i = 0; i < numRowIndices * numColIndices; i++) pMat[i] = pMatModified[i];  in DMPlexComputeInterpolatorTree()
2333 …Offsets[f]; j < newOffsets[f + 1]; j++, count++) pMat[count] = pMatModified[i * numColIndices + j];  in DMPlexComputeInterpolatorTree()
2337 …  PetscCall(DMRestoreWorkArray(coarse, numRowIndices * numColIndices, MPIU_SCALAR, &pMatModified));  in DMPlexComputeInterpolatorTree()
2339 …        PetscCall(DMRestoreWorkArray(coarse, numRowIndices * numColIndices, MPIU_SCALAR, &pMatIn));  in DMPlexComputeInterpolatorTree()
2342               pInd[numColIndices + f]             = offsets[f + 1];  in DMPlexComputeInterpolatorTree()
2343               pInd[numColIndices + numFields + f] = newOffsets[f + 1];  in DMPlexComputeInterpolatorTree()
2371         PetscCall(DMGetWorkArray(coarse, numColIndices, MPIU_INT, &colIndices));  in DMPlexComputeInterpolatorTree()
2415             pInd[numColIndices + f]             = offsets[f + 1];  in DMPlexComputeInterpolatorTree()
2416             pInd[numColIndices + numFields + f] = newOffsets[f + 1];  in DMPlexComputeInterpolatorTree()
2426           PetscCall(MatGetValues(cMat, numRowIndices, rowIndices, numColIndices, colIndices, pMat));  in DMPlexComputeInterpolatorTree()
2434         PetscCall(DMRestoreWorkArray(coarse, numColIndices, MPIU_INT, &colIndices));  in DMPlexComputeInterpolatorTree()
2447             pInd[numColIndices + f]             = offsets[f + 1];  in DMPlexComputeInterpolatorTree()
2448             pInd[numColIndices + numFields + f] = offsets[f + 1];  in DMPlexComputeInterpolatorTree()
2512       PetscInt numColIndices, pIndOff, *pInd;  in DMPlexComputeInterpolatorTree()  local
2522       PetscCall(PetscSectionGetDof(leafIndicesSec, p, &numColIndices));  in DMPlexComputeInterpolatorTree()
2524       numColIndices -= 2 * numFields;  in DMPlexComputeInterpolatorTree()
2525 …PetscCheck(numColIndices > 0, PETSC_COMM_SELF, PETSC_ERR_PLIB, "global fine dof with no dofs to in…  in DMPlexComputeInterpolatorTree()
2539           newOffsets[f + 1]  = pInd[numColIndices + numFields + f];  in DMPlexComputeInterpolatorTree()
2562         for (i = 0; i < numColIndices; i++) {  in DMPlexComputeInterpolatorTree()
2681       PetscInt numColIndices, pIndOff, *pInd;  in DMPlexComputeInterpolatorTree()  local
2690       PetscCall(PetscSectionGetDof(leafIndicesSec, p, &numColIndices));  in DMPlexComputeInterpolatorTree()
2692       numColIndices -= 2 * numFields;  in DMPlexComputeInterpolatorTree()
2707           rowOffsets[f + 1]  = pInd[numColIndices + f];  in DMPlexComputeInterpolatorTree()
2708           newOffsets[f + 1]  = pInd[numColIndices + numFields + f];  in DMPlexComputeInterpolatorTree()
2736 …PetscCall(MatSetValues(mat, gDof, rowIndices, numColIndices, pInd, refPointFieldMats[childId - pRe…  in DMPlexComputeInterpolatorTree()
2758 …          PetscCall(MatSetValues(mat, gDof, rowIndices, numColIndices, pInd, pMat, INSERT_VALUES));  in DMPlexComputeInterpolatorTree()
2782             PetscInt numCols   = numColIndices;  in DMPlexComputeInterpolatorTree()
2783             PetscInt numInRows = matSize / numColIndices;  in DMPlexComputeInterpolatorTree()