543       PetscInt  numFIndices, numCIndices;  in DMSwarmComputeMassMatrix_Private()  local
545 …exGetClosureIndices(dmf, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeMassMatrix_Private()
551         for (PetscInt i = 0; i < numFIndices; ++i) {  in DMSwarmComputeMassMatrix_Private()
570 …storeClosureIndices(dmf, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeMassMatrix_Private()
592       PetscInt  numFIndices, numCIndices;  in DMSwarmComputeMassMatrix_Private()  local
596 …exGetClosureIndices(dmf, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeMassMatrix_Private()
612       for (PetscInt i = 0; i < numFIndices / Nc; ++i) {  in DMSwarmComputeMassMatrix_Private()
617 …elemMat[(j * totNc + c) * numFIndices + i * Nc + c] += Tcoarse->T[0][(j * numFIndices + i * Nc + c…  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()
627 …storeClosureIndices(dmf, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeMassMatrix_Private()
818       PetscInt  numFIndices, numCIndices;  in DMSwarmComputeMassMatrixSquare_Private()  local
822 …exGetClosureIndices(dmf, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeMassMatrixSquare_Private()
837       for (PetscInt i = 0; i < numFIndices; ++i) {  in DMSwarmComputeMassMatrixSquare_Private()
841 …elemMat[p * numFIndices + i] += Tcoarse->T[0][(p * numFIndices + i) * Nc + c] * (useDeltaFunction …  in DMSwarmComputeMassMatrixSquare_Private()
854         PetscCall(PetscBLASIntCast(numFIndices, &blask));  in DMSwarmComputeMassMatrixSquare_Private()
860 …storeClosureIndices(dmf, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeMassMatrixSquare_Private()
1018       PetscInt  numFIndices, Npc;  in DMSwarmComputeGradientMatrix_Private()  local
1022 …lexGetClosureIndices(dm, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeGradientMatrix_Private()
1037       for (PetscInt i = 0; i < numFIndices; ++i) {  in DMSwarmComputeGradientMatrix_Private()
1042 …e = 0; e < cdim; ++e) xi[d] += invJ[e * cdim + d] * Tcoarse->T[1][(j * numFIndices + i) * cdim + e…  in DMSwarmComputeGradientMatrix_Private()
1043             elemMat[(j * cdim + d) * numFIndices + i] += xi[d] * (useDeltaFunction ? 1.0 : detJ);  in DMSwarmComputeGradientMatrix_Private()
1049       if (0) PetscCall(DMPrintCellMatrix(cell, name, Npc * cdim, numFIndices, elemMat));  in DMSwarmComputeGradientMatrix_Private()
1050 …    PetscCall(MatSetValues(derv, Npc * cdim, rowIDXs, numFIndices, findices, elemMat, ADD_VALUES));  in DMSwarmComputeGradientMatrix_Private()
1052 …estoreClosureIndices(dm, fsection, globalFSection, cell, PETSC_FALSE, &numFIndices, &findices, NUL…  in DMSwarmComputeGradientMatrix_Private()