543       PetscInt  numFIndices, numCIndices;  in DMSwarmComputeMassMatrix_Private()  local
546       PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
547       maxC = PetscMax(maxC, numCIndices);  in DMSwarmComputeMassMatrix_Private()
555             for (PetscInt j = 0; j < numCIndices; ++j) {  in DMSwarmComputeMassMatrix_Private()
568         PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
592       PetscInt  numFIndices, numCIndices;  in DMSwarmComputeMassMatrix_Private()  local
597       PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
598       for (PetscInt j = 0; j < numCIndices; ++j) {  in DMSwarmComputeMassMatrix_Private()
608 …if (id == PETSCFE_CLASSID) PetscCall(PetscFECreateTabulation((PetscFE)obj, 1, numCIndices, xi, 0, …  in DMSwarmComputeMassMatrix_Private()
609       else PetscCall(PetscFVCreateTabulation((PetscFV)obj, 1, numCIndices, xi, 0, &Tcoarse));  in DMSwarmComputeMassMatrix_Private()
611       PetscCall(PetscArrayzero(elemMat, numCIndices * Nc * totDim));  in DMSwarmComputeMassMatrix_Private()
613         for (PetscInt j = 0; j < numCIndices; ++j) {  in DMSwarmComputeMassMatrix_Private()
621       for (PetscInt j = 0; j < numCIndices; ++j)  in DMSwarmComputeMassMatrix_Private()
624       if (0) PetscCall(DMPrintCellMatrix(cell, name, numCIndices * Nc, numFIndices, elemMat));  in DMSwarmComputeMassMatrix_Private()
625 …PetscCall(MatSetValues(mass, numCIndices * Nc, rowIDXs, numFIndices, findices, elemMat, ADD_VALUES…  in DMSwarmComputeMassMatrix_Private()
626       PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
750     PetscInt  numCIndices;  in DMSwarmComputeMassMatrixSquare_Private()  local
755     PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()
756     maxC = PetscMax(maxC, numCIndices);  in DMSwarmComputeMassMatrixSquare_Private()
758     for (PetscInt i = 0; i < numCIndices; ++i) dnz[cindices[i]] += numCIndices;  in DMSwarmComputeMassMatrixSquare_Private()
773           for (i = 0; i < numCIndices; ++i) {  in DMSwarmComputeMassMatrixSquare_Private()
791     PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()
818       PetscInt  numFIndices, numCIndices;  in DMSwarmComputeMassMatrixSquare_Private()  local
823       PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()
824       for (PetscInt p = 0; p < numCIndices; ++p) {  in DMSwarmComputeMassMatrixSquare_Private()
834       PetscCall(PetscFECreateTabulation((PetscFE)obj, 1, numCIndices, xi, 0, &Tcoarse));  in DMSwarmComputeMassMatrixSquare_Private()
836       PetscCall(PetscArrayzero(elemMat, numCIndices * totDim));  in DMSwarmComputeMassMatrixSquare_Private()
838         for (PetscInt p = 0; p < numCIndices; ++p) {  in DMSwarmComputeMassMatrixSquare_Private()
846       for (PetscInt p = 0; p < numCIndices; ++p) rowIDXs[p] = cindices[p] + rStart;  in DMSwarmComputeMassMatrixSquare_Private()
847       if (0) PetscCall(DMPrintCellMatrix(cell, name, 1, numCIndices, elemMat));  in DMSwarmComputeMassMatrixSquare_Private()
849       if (numCIndices) {  in DMSwarmComputeMassMatrixSquare_Private()
853         PetscCall(PetscBLASIntCast(numCIndices, &blasn));  in DMSwarmComputeMassMatrixSquare_Private()
857 …  PetscCall(MatSetValues(mass, numCIndices, rowIDXs, numCIndices, rowIDXs, elemMatSq, ADD_VALUES));  in DMSwarmComputeMassMatrixSquare_Private()
859       PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()