1011     PetscReal      v0[3], J[9], invJ[9], detJ, detJp = 2. / (gmax[0] - gmin[0]), xr[3], den;  in DMSwarmComputeLocalSize()  local
1014     PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, NULL, v0, J, invJ, &detJ));  in DMSwarmComputeLocalSize()
1017       CoordinatesRefToReal(dim, dim, xi0, v0, J, &xq[q * dim], xr);  in DMSwarmComputeLocalSize()
1167       PetscReal v0[3], J[9], invJ[9], detJ;  in DMSwarmInitializeCoordinates()  local
1171       PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ));  in DMSwarmInitializeCoordinates()
1177         CoordinatesRefToReal(dim, dim, xi0, v0, J, xref, &x[p * dim]);  in DMSwarmInitializeCoordinates()
1213 PetscErrorCode DMSwarmInitializeVelocities(DM sw, PetscProbFn *sampler, const PetscReal v0[])  in DMSwarmInitializeVelocities()  argument
1228   if (v0[0] == 0.) {  in DMSwarmInitializeVelocities()
1237       for (d = 0; d < dim; ++d) v[p * dim + d] = (v0[s] / v0[0]) * PetscRealPart(vel[d]);  in DMSwarmInitializeVelocities()
1252       for (d = 0; d < dim; ++d) v[p * dim + d] = (v0[s] / v0[0]) * vel[d];  in DMSwarmInitializeVelocities()
1274 PetscErrorCode DMSwarmInitializeVelocitiesFromOptions(DM sw, const PetscReal v0[])  in DMSwarmInitializeVelocitiesFromOptions()  argument
1296   PetscCall(DMSwarmInitializeVelocities(sw, sampler, v0));  in DMSwarmInitializeVelocitiesFromOptions()
1310   PetscReal       *xi, *v0, *J, *invJ, detJ = 1.0, v0ref[3] = {-1.0, -1.0, -1.0};  in DMProjectFieldLocal_Swarm()  local
1327   PetscCall(PetscMalloc3(dim, &v0, dim * dim, &J, dim * dim, &invJ));  in DMProjectFieldLocal_Swarm()
1332   fegeom.v        = v0;  in DMProjectFieldLocal_Swarm()
1368       PetscCall(DMPlexComputeCellGeometryFEM(dmIn, cell, NULL, v0, J, invJ, &detJ));  in DMProjectFieldLocal_Swarm()
1413   PetscCall(PetscFree3(v0, J, invJ));  in DMProjectFieldLocal_Swarm()