1650   PetscFEGeom    *chunkgeom = NULL;  in ComputeFieldAtParticles_Primal()  local
1671       PetscCall(PetscFEGeomGetChunk(ctx->fegeom, c - cStart, c - cStart + 1, &chunkgeom));  in ComputeFieldAtParticles_Primal()
1674 …CoordinatesRealToRef(dim, dim, x0, chunkgeom->v, chunkgeom->invJ, &pcoord[dim * i], &refcoord[dim …  in ComputeFieldAtParticles_Primal()
1684 …PetscCall(PetscFEFreeInterpolateGradient_Static(fe, basisDer, clPhi, dim, chunkgeom->invJ, NULL, c…  in ComputeFieldAtParticles_Primal()
1696   PetscCall(PetscFEGeomRestoreChunk(ctx->fegeom, 0, 1, &chunkgeom));  in ComputeFieldAtParticles_Primal()
1776   PetscFEGeom    *chunkgeom = NULL;  in ComputeFieldAtParticles_Mixed()  local
1797       PetscCall(PetscFEGeomGetChunk(ctx->fegeom, c - cStart, c - cStart + 1, &chunkgeom));  in ComputeFieldAtParticles_Mixed()
1800 …CoordinatesRealToRef(dim, dim, x0, chunkgeom->v, chunkgeom->invJ, &pcoord[dim * i], &refcoord[dim …  in ComputeFieldAtParticles_Mixed()
1809       PetscCall(PetscFEInterpolateAtPoints_Static(fe, tab, clPhi, chunkgeom, cp, &E[p * dim]));  in ComputeFieldAtParticles_Mixed()
1810       PetscCall(PetscFEPushforward(fe, chunkgeom, 1, &E[p * dim]));  in ComputeFieldAtParticles_Mixed()
1822   PetscCall(PetscFEGeomRestoreChunk(ctx->fegeom, 0, 1, &chunkgeom));  in ComputeFieldAtParticles_Mixed()