1425       PetscInt cDof, cOff, numCols, r, i;  in DMPlexReferenceTreeGetChildrenMatrices()  local
1438       numCols = 0;  in DMPlexReferenceTreeGetChildrenMatrices()
1452         for (j = 0; j < aDof; j++) cols[numCols++] = aOff + (perm ? perm[j] : j);  in DMPlexReferenceTreeGetChildrenMatrices()
1454       refPointFieldN[p - pRefStart][f] = numCols;  in DMPlexReferenceTreeGetChildrenMatrices()
1455       PetscCall(PetscMalloc1(cDof * numCols, &refPointFieldMats[p - pRefStart][f]));  in DMPlexReferenceTreeGetChildrenMatrices()
1456 …  PetscCall(MatGetValues(refCmat, cDof, rows, numCols, cols, refPointFieldMats[p - pRefStart][f]));  in DMPlexReferenceTreeGetChildrenMatrices()
1475 …for (j = 0; j < aDof; j++) refPointFieldMats[p - pRefStart][f][k * numCols + colOff + j] *= flip[j…  in DMPlexReferenceTreeGetChildrenMatrices()
1638         PetscInt            cDof, cOff, numCols, numFillCols, i, r, matOffset, offset;  in DMPlexComputeAnchorMatrix_Tree_FromReference()  local
1667         numCols     = refPointFieldN[childid - pRefStart][f];  in DMPlexComputeAnchorMatrix_Tree_FromReference()
1705 …inVal += pointMat[r * numCols + offset + col] * vals[aMatOffset + aNumFillCols * k + j] * (flip ? …  in DMPlexComputeAnchorMatrix_Tree_FromReference()
1729 …vals[matOffset + numFillCols * r + k + col] += pointMat[r * numCols + offset + j] * (flip ? flip[c…  in DMPlexComputeAnchorMatrix_Tree_FromReference()
2732               PetscInt numCols = newOffsets[f + 1] - newOffsets[f];  in DMPlexComputeInterpolatorTree()  local
2733 …PetscCall(MatSetValues(mat, numRows, &rowIndices[offsets[f]], numCols, &pInd[newOffsets[f]], refPo…  in DMPlexComputeInterpolatorTree()
2750               PetscInt     numCols   = newOffsets[f + 1] - newOffsets[f];  in DMPlexComputeInterpolatorTree()  local
2754 …PetscCall(MatSetValues(mat, numRows, &rowIndices[offsets[f]], numCols, &pInd[newOffsets[f]], inMat…  in DMPlexComputeInterpolatorTree()
2755               count += numCols * numInRows;  in DMPlexComputeInterpolatorTree()
2765               PetscInt     numCols   = newOffsets[f + 1] - newOffsets[f];  in DMPlexComputeInterpolatorTree()  local
2771                 for (j = 0; j < numCols; j++) {  in DMPlexComputeInterpolatorTree()
2773 … k++) val += refPointFieldMats[childId - pRefStart][f][i * numInRows + k] * inMat[k * numCols + j];  in DMPlexComputeInterpolatorTree()
2774                   pointWork[i * numCols + j] = val;  in DMPlexComputeInterpolatorTree()
2777 …PetscCall(MatSetValues(mat, numRows, &rowIndices[offsets[f]], numCols, &pInd[newOffsets[f]], point…  in DMPlexComputeInterpolatorTree()
2778               count += numCols * numInRows;  in DMPlexComputeInterpolatorTree()
2782             PetscInt numCols   = numColIndices;  in DMPlexComputeInterpolatorTree()  local
2787               for (j = 0; j < numCols; j++) {  in DMPlexComputeInterpolatorTree()
2789 …; k++) val += refPointFieldMats[childId - pRefStart][0][i * numInRows + k] * pMat[k * numCols + j];  in DMPlexComputeInterpolatorTree()
2790                 pointWork[i * numCols + j] = val;  in DMPlexComputeInterpolatorTree()
2793 …        PetscCall(MatSetValues(mat, numRows, rowIndices, numCols, pInd, pointWork, INSERT_VALUES));  in DMPlexComputeInterpolatorTree()
3184       PetscInt cDof, cOff, numCols, r;  in DMPlexReferenceTreeGetChildrenMatrices_Injection()  local
3195       numCols = 0;  in DMPlexReferenceTreeGetChildrenMatrices_Injection()
3207         for (j = 0; j < aDof; j++) cols[numCols++] = aOff + j;  in DMPlexReferenceTreeGetChildrenMatrices_Injection()
3209       PetscCall(PetscMalloc1(cDof * numCols, &refPointFieldMats[p - pRefStart][f]));  in DMPlexReferenceTreeGetChildrenMatrices_Injection()
3211       PetscCall(MatGetValues(inj, numCols, cols, cDof, rows, refPointFieldMats[p - pRefStart][f]));  in DMPlexReferenceTreeGetChildrenMatrices_Injection()
4016           PetscInt           numCols  = newOffsets[f + 1] - newOffsets[f];  in DMPlexTransferVecTree_Interpolate()  local
4022 …ointFieldN %" PetscInt_FMT " maxDof %" PetscInt_FMT "\n",childId,numRows,numCols,refPointFieldN[ch…  in DMPlexTransferVecTree_Interpolate()
4026             for (j = 0; j < numCols; j++) val += childMat[i * numCols + j] * cVal[j];  in DMPlexTransferVecTree_Interpolate()