702   PetscInt      cdim, Nf, Nfc, cStart, cEnd, totDim, maxC = 0;  in DMSwarmComputeMassMatrixSquare_Private()  local
709   PetscCall(DMGetCoordinateDim(dmf, &cdim));  in DMSwarmComputeMassMatrixSquare_Private()
713   PetscCall(PetscMalloc3(cdim, &v0, cdim * cdim, &J, cdim * cdim, &invJ));  in DMSwarmComputeMassMatrixSquare_Private()
797 …(PetscMalloc4(maxC * totDim, &elemMat, maxC * maxC, &elemMatSq, maxC, &rowIDXs, maxC * cdim, &xi));  in DMSwarmComputeMassMatrixSquare_Private()
831 …cdim, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "The total block size of coordinates is %" PetscInt_FM…  in DMSwarmComputeMassMatrixSquare_Private()
832         CoordinatesRealToRef(cdim, cdim, v0ref, v0, invJ, xr, &xi[p * cdim]);  in DMSwarmComputeMassMatrixSquare_Private()
953   PetscInt      cdim, Nf, Nfc, cStart, cEnd, totDim, maxNpc = 0, totNc = 0;  in DMSwarmComputeGradientMatrix_Private()  local
960   PetscCall(DMGetCoordinateDim(dm, &cdim));  in DMSwarmComputeGradientMatrix_Private()
965   PetscCall(PetscMalloc3(cdim, &v0, cdim * cdim, &J, cdim * cdim, &invJ));  in DMSwarmComputeGradientMatrix_Private()
978   PetscCall(PetscLayoutSetBlockSize(rLayout, cdim));  in DMSwarmComputeGradientMatrix_Private()
1009 …PetscCall(PetscMalloc3(maxNpc * cdim * totDim, &elemMat, maxNpc * cdim, &rowIDXs, maxNpc * cdim, &…  in DMSwarmComputeGradientMatrix_Private()
1031 …cdim, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "The total block size of coordinates is %" PetscInt_FM…  in DMSwarmComputeGradientMatrix_Private()
1032         CoordinatesRealToRef(cdim, cdim, v0ref, v0, invJ, xr, &xi[j * cdim]);  in DMSwarmComputeGradientMatrix_Private()
1036       PetscCall(PetscArrayzero(elemMat, Npc * cdim * totDim));  in DMSwarmComputeGradientMatrix_Private()
1040           for (PetscInt d = 0; d < cdim; ++d) {  in DMSwarmComputeGradientMatrix_Private()
1042 …for (PetscInt e = 0; e < cdim; ++e) xi[d] += invJ[e * cdim + d] * Tcoarse->T[1][(j * numFIndices +…  in DMSwarmComputeGradientMatrix_Private()
1043             elemMat[(j * cdim + d) * numFIndices + i] += xi[d] * (useDeltaFunction ? 1.0 : detJ);  in DMSwarmComputeGradientMatrix_Private()
1048         for (PetscInt d = 0; d < cdim; ++d) rowIDXs[j * cdim + d] = pind[j] * cdim + d + rStart;  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()
1073   PetscInt      cdim, m, n, Np, bs;  in DMCreateGradientMatrix_Swarm()  local
1079   PetscCall(DMGetCoordinateDim(dm, &cdim));  in DMCreateGradientMatrix_Swarm()
1086 …PetscCheck(cdim == bs, comm, PETSC_ERR_ARG_WRONG, "Coordinate dimension %" PetscInt_FMT " != %" Pe…  in DMCreateGradientMatrix_Swarm()