478   PetscScalar  *elemMat;  in DMSwarmComputeMassMatrix_Private()  local
577   PetscCall(PetscMalloc3(maxC * totNc * totDim, &elemMat, maxC * totNc, &rowIDXs, maxC * dim, &xi));  in DMSwarmComputeMassMatrix_Private()
611       PetscCall(PetscArrayzero(elemMat, numCIndices * Nc * totDim));  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()
632   PetscCall(PetscFree3(elemMat, rowIDXs, xi));  in DMSwarmComputeMassMatrix_Private()
701   PetscScalar  *elemMat, *elemMatSq;  in DMSwarmComputeMassMatrixSquare_Private()  local
797 …PetscCall(PetscMalloc4(maxC * totDim, &elemMat, maxC * maxC, &elemMatSq, maxC, &rowIDXs, maxC * cd…  in DMSwarmComputeMassMatrixSquare_Private()
836       PetscCall(PetscArrayzero(elemMat, numCIndices * totDim));  in DMSwarmComputeMassMatrixSquare_Private()
841 …elemMat[p * numFIndices + i] += Tcoarse->T[0][(p * numFIndices + i) * Nc + c] * (useDeltaFunction …  in DMSwarmComputeMassMatrixSquare_Private()
847       if (0) PetscCall(DMPrintCellMatrix(cell, name, 1, numCIndices, elemMat));  in DMSwarmComputeMassMatrixSquare_Private()
855 …S("BLASgemm", BLASgemm_("T", "N", &blasn, &blasn, &blask, &one, elemMat, &blask, elemMat, &blask, …  in DMSwarmComputeMassMatrixSquare_Private()
864   PetscCall(PetscFree4(elemMat, elemMatSq, rowIDXs, xi));  in DMSwarmComputeMassMatrixSquare_Private()
952   PetscScalar  *elemMat;  in DMSwarmComputeGradientMatrix_Private()  local
1009 …PetscCall(PetscMalloc3(maxNpc * cdim * totDim, &elemMat, maxNpc * cdim, &rowIDXs, maxNpc * cdim, &…  in DMSwarmComputeGradientMatrix_Private()
1036       PetscCall(PetscArrayzero(elemMat, Npc * cdim * totDim));  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()
1057   PetscCall(PetscFree3(elemMat, rowIDXs, xi));  in DMSwarmComputeGradientMatrix_Private()