542       PetscInt *findices, *cindices; /* fine is vertices, coarse is particles */  in DMSwarmComputeMassMatrix_Private()  local
546       PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
558                 key.i = cindices[j] * totNc + c + rStart; /* global cols (from Swarm) */  in DMSwarmComputeMassMatrix_Private()
568         PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
591       PetscInt *findices, *cindices;  in DMSwarmComputeMassMatrix_Private()  local
597       PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
603 …       for (PetscInt b = 0; b < bs[i]; ++b, ++off) xr[off] = coordVals[i][cindices[j] * bs[i] + b];  in DMSwarmComputeMassMatrix_Private()
623         for (PetscInt c = 0; c < Nc; ++c) rowIDXs[j * Nc + c] = cindices[j] * totNc + c + rStart;  in DMSwarmComputeMassMatrix_Private()
626       PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrix_Private()
749     PetscInt *cindices;  in DMSwarmComputeMassMatrixSquare_Private()  local
755     PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()
758     for (PetscInt i = 0; i < numCIndices; ++i) dnz[cindices[i]] += numCIndices;  in DMSwarmComputeMassMatrixSquare_Private()
774             key.i = cindices[i] + rStart; /* global rows (from Swarm) */  in DMSwarmComputeMassMatrixSquare_Private()
791     PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()
817       PetscInt *findices, *cindices;  in DMSwarmComputeMassMatrixSquare_Private()  local
823       PetscCall(DMSwarmSortGetPointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()
829 …       for (PetscInt b = 0; b < bs[i]; ++b, ++off) xr[off] = coordVals[i][cindices[p] * bs[i] + b];  in DMSwarmComputeMassMatrixSquare_Private()
846       for (PetscInt p = 0; p < numCIndices; ++p) rowIDXs[p] = cindices[p] + rStart;  in DMSwarmComputeMassMatrixSquare_Private()
859       PetscCall(DMSwarmSortRestorePointsPerCell(dmc, cell, &numCIndices, &cindices));  in DMSwarmComputeMassMatrixSquare_Private()