221   PetscInt     dim, p, pStart, pEnd, pRefStart, pRefEnd, d, offset, parentSize, *parents, *childIDs;  in DMPlexCreateReferenceTree_Union()  local
234   PetscCall(DMPlexGetChart(Kref, &pRefStart, &pRefEnd));  in DMPlexCreateReferenceTree_Union()
236   PetscCall(PetscSectionSetChart(unionSection, 0, (pEnd - pStart) + (pRefEnd - pRefStart)));  in DMPlexCreateReferenceTree_Union()
239   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexCreateReferenceTree_Union()
243     if (qSize > 1) PetscCall(PetscSectionSetDof(unionSection, p - pRefStart + (pEnd - pStart), 1));  in DMPlexCreateReferenceTree_Union()
245   PetscCall(PetscMalloc1(pEnd - pStart + pRefEnd - pRefStart, &permvals));  in DMPlexCreateReferenceTree_Union()
257 …PetscCall(ISCreateGeneral(comm, (pEnd - pStart) + (pRefEnd - pRefStart), permvals, PETSC_OWN_POINT…  in DMPlexCreateReferenceTree_Union()
286   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexCreateReferenceTree_Union()
290     PetscCall(PetscSectionGetDof(unionSection, p - pRefStart + (pEnd - pStart), &uDof));  in DMPlexCreateReferenceTree_Union()
291     PetscCall(PetscSectionGetOffset(unionSection, p - pRefStart + (pEnd - pStart), &uOff));  in DMPlexCreateReferenceTree_Union()
318   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexCreateReferenceTree_Union()
325     PetscCall(PetscSectionGetDof(unionSection, p - pRefStart + (pEnd - pStart), &uDof));  in DMPlexCreateReferenceTree_Union()
326     PetscCall(PetscSectionGetOffset(unionSection, p - pRefStart + (pEnd - pStart), &uOff));  in DMPlexCreateReferenceTree_Union()
333         PetscCall(PetscSectionGetDof(unionSection, e - pRefStart + (pEnd - pStart), &eDof));  in DMPlexCreateReferenceTree_Union()
335           PetscCall(PetscSectionGetOffset(unionSection, e - pRefStart + (pEnd - pStart), &eOff));  in DMPlexCreateReferenceTree_Union()
370       PetscCall(PetscSectionGetDof(unionSection, v - pRefStart + (pEnd - pStart), &vDof));  in DMPlexCreateReferenceTree_Union()
371       PetscCall(PetscSectionGetOffset(unionSection, v - pRefStart + (pEnd - pStart), &vOff));  in DMPlexCreateReferenceTree_Union()
386   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexCreateReferenceTree_Union()
389     PetscCall(PetscSectionGetDof(unionSection, p - pRefStart + (pEnd - pStart), &uDof));  in DMPlexCreateReferenceTree_Union()
390     PetscCall(PetscSectionGetOffset(unionSection, p - pRefStart + (pEnd - pStart), &uOff));  in DMPlexCreateReferenceTree_Union()
396   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexCreateReferenceTree_Union()
399     PetscCall(PetscSectionGetDof(unionSection, p - pRefStart + (pEnd - pStart), &uDof));  in DMPlexCreateReferenceTree_Union()
400     PetscCall(PetscSectionGetOffset(unionSection, p - pRefStart + (pEnd - pStart), &uOff));  in DMPlexCreateReferenceTree_Union()
1391 …PetscInt            numFields, maxFields, f, pRefStart, pRefEnd, p, *rows, *cols, maxDof, maxAnDof…  in DMPlexReferenceTreeGetChildrenMatrices()  local
1407   PetscCall(PetscSectionGetChart(refConSec, &pRefStart, &pRefEnd));  in DMPlexReferenceTreeGetChildrenMatrices()
1408   PetscCall(PetscMalloc1(pRefEnd - pRefStart, &refPointFieldMats));  in DMPlexReferenceTreeGetChildrenMatrices()
1409   PetscCall(PetscMalloc1(pRefEnd - pRefStart, &refPointFieldN));  in DMPlexReferenceTreeGetChildrenMatrices()
1414   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexReferenceTreeGetChildrenMatrices()
1421     PetscCall(PetscMalloc1(maxFields, &refPointFieldMats[p - pRefStart]));  in DMPlexReferenceTreeGetChildrenMatrices()
1422     PetscCall(PetscCalloc1(maxFields, &refPointFieldN[p - pRefStart]));  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()
1502   PetscInt       numFields, maxFields, pRefStart, pRefEnd, p, f;  in DMPlexReferenceTreeRestoreChildrenMatrices()  local
1514   PetscCall(PetscSectionGetChart(refConSec, &pRefStart, &pRefEnd));  in DMPlexReferenceTreeRestoreChildrenMatrices()
1515   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexReferenceTreeRestoreChildrenMatrices()
1531       PetscCall(PetscFree(refPointFieldMats[p - pRefStart][f]));  in DMPlexReferenceTreeRestoreChildrenMatrices()
1533     PetscCall(PetscFree(refPointFieldMats[p - pRefStart]));  in DMPlexReferenceTreeRestoreChildrenMatrices()
1534     PetscCall(PetscFree(refPointFieldN[p - pRefStart]));  in DMPlexReferenceTreeRestoreChildrenMatrices()
1546 …PetscInt        numFields, maxFields, f, pRefStart, pRefEnd, p, maxDof, maxAnDof, *perm, *iperm, p…  in DMPlexComputeAnchorMatrix_Tree_FromReference()  local
1565   PetscCall(PetscSectionGetChart(refConSec, &pRefStart, &pRefEnd));  in DMPlexComputeAnchorMatrix_Tree_FromReference()
1666         pointMat    = refPointFieldMats[childid - pRefStart][f];  in DMPlexComputeAnchorMatrix_Tree_FromReference()
1667         numCols     = refPointFieldN[childid - pRefStart][f];  in DMPlexComputeAnchorMatrix_Tree_FromReference()
2497     PetscInt       pRefStart, pRefEnd, maxConDof, maxColumns, leafStart, leafEnd;  in DMPlexComputeInterpolatorTree()  local
2675     PetscCall(PetscSectionGetChart(refConSec, &pRefStart, &pRefEnd));  in DMPlexComputeInterpolatorTree()
2733 …offsets[f]], numCols, &pInd[newOffsets[f]], refPointFieldMats[childId - pRefStart][f], INSERT_VALU…  in DMPlexComputeInterpolatorTree()
2736 …mat, gDof, rowIndices, numColIndices, pInd, refPointFieldMats[childId - pRefStart][0], INSERT_VALU…  in DMPlexComputeInterpolatorTree()
2769 …PetscCheck(refPointFieldN[childId - pRefStart][f] == numInRows, PETSC_COMM_SELF, PETSC_ERR_PLIB, "…  in DMPlexComputeInterpolatorTree()
2773 …for (k = 0; k < numInRows; k++) val += refPointFieldMats[childId - pRefStart][f][i * numInRows + k…  in DMPlexComputeInterpolatorTree()
2785 …PetscCheck(refPointFieldN[childId - pRefStart][0] == numInRows, PETSC_COMM_SELF, PETSC_ERR_PLIB, "…  in DMPlexComputeInterpolatorTree()
2789 …for (k = 0; k < numInRows; k++) val += refPointFieldMats[childId - pRefStart][0][i * numInRows + k…  in DMPlexComputeInterpolatorTree()
3160   PetscInt       numFields, f, pRefStart, pRefEnd, p, *rows, *cols, maxDof;  in DMPlexReferenceTreeGetChildrenMatrices_Injection()  local
3169   PetscCall(PetscSectionGetChart(refConSec, &pRefStart, &pRefEnd));  in DMPlexReferenceTreeGetChildrenMatrices_Injection()
3170   PetscCall(PetscMalloc1(pRefEnd - pRefStart, &refPointFieldMats));  in DMPlexReferenceTreeGetChildrenMatrices_Injection()
3174   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexReferenceTreeGetChildrenMatrices_Injection()
3182     PetscCall(PetscMalloc1(numFields, &refPointFieldMats[p - pRefStart]));  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()
3224   PetscInt       numFields, pRefStart, pRefEnd, p, f;  in DMPlexReferenceTreeRestoreChildrenMatrices_Injection()  local
3234   PetscCall(PetscSectionGetChart(refConSec, &pRefStart, &pRefEnd));  in DMPlexReferenceTreeRestoreChildrenMatrices_Injection()
3235   for (p = pRefStart; p < pRefEnd; p++) {  in DMPlexReferenceTreeRestoreChildrenMatrices_Injection()
3252       PetscCall(PetscFree(refPointFieldMats[p - pRefStart][f]));  in DMPlexReferenceTreeRestoreChildrenMatrices_Injection()
3254     PetscCall(PetscFree(refPointFieldMats[p - pRefStart]));  in DMPlexReferenceTreeRestoreChildrenMatrices_Injection()
3517   PetscInt      *rootIndices = NULL, *parentIndices, pRefStart, pRefEnd;  in DMPlexComputeInjectorTree()  local
3533   PetscCall(PetscSectionGetChart(cSecRef, &pRefStart, &pRefEnd));  in DMPlexComputeInjectorTree()
3728           PetscScalar    *childMat   = &childrenMats[childId - pRefStart][f][0];  in DMPlexComputeInjectorTree()
3950     PetscInt       pRefStart, pRefEnd, leafStart, leafEnd;  in DMPlexTransferVecTree_Interpolate()  local
3961     PetscCall(PetscSectionGetChart(refConSec, &pRefStart, &pRefEnd));  in DMPlexTransferVecTree_Interpolate()
3994 …fsets[f + 1] = newOffsets[f] + ((childId == -1) ? rowDof : refPointFieldN[childId - pRefStart][f]);  in DMPlexTransferVecTree_Interpolate()
4001         newOffsets[1] = (childId == -1) ? lDof : refPointFieldN[childId - pRefStart][0];  in DMPlexTransferVecTree_Interpolate()
4014           const PetscScalar *childMat = refPointFieldMats[childId - pRefStart][f];  in DMPlexTransferVecTree_Interpolate()
4022 …xDof %" PetscInt_FMT "\n",childId,numRows,numCols,refPointFieldN[childId - pRefStart][f], maxDof));  in DMPlexTransferVecTree_Interpolate()
4066   PetscInt      *parentIndices, pRefStart, pRefEnd;  in DMPlexTransferVecTree_Inject()  local
4083   PetscCall(PetscSectionGetChart(cSecRef, &pRefStart, &pRefEnd));  in DMPlexTransferVecTree_Inject()
4183           PetscScalar       *childMat = &childrenMats[childId - pRefStart][f][0];  in DMPlexTransferVecTree_Inject()