2251   PetscFEGeom    *chunkgeom = NULL;  in ComputeFieldAtParticles_Primal()  local
2275       PetscCall(PetscFEGeomGetChunk(ctx->fegeom, c - cStart, c - cStart + 1, &chunkgeom));  in ComputeFieldAtParticles_Primal()
2278 …CoordinatesRealToRef(dim, dim, x0, chunkgeom->v, chunkgeom->invJ, &pcoord[dim * i], &refcoord[dim …  in ComputeFieldAtParticles_Primal()
2288 …PetscCall(PetscFEFreeInterpolateGradient_Static(fe, basisDer, clPhi, dim, chunkgeom->invJ, NULL, c…  in ComputeFieldAtParticles_Primal()
2300   PetscCall(PetscFEGeomRestoreChunk(ctx->fegeom, 0, 1, &chunkgeom));  in ComputeFieldAtParticles_Primal()
2381   PetscFEGeom    *chunkgeom = NULL;  in ComputeFieldAtParticles_Mixed()  local
2402       PetscCall(PetscFEGeomGetChunk(ctx->fegeom, c - cStart, c - cStart + 1, &chunkgeom));  in ComputeFieldAtParticles_Mixed()
2405 …CoordinatesRealToRef(dim, dim, x0, chunkgeom->v, chunkgeom->invJ, &pcoord[dim * i], &refcoord[dim …  in ComputeFieldAtParticles_Mixed()
2414       PetscCall(PetscFEInterpolateAtPoints_Static(fe, tab, clPhi, chunkgeom, cp, &E[p * dim]));  in ComputeFieldAtParticles_Mixed()
2415       PetscCall(PetscFEPushforward(fe, chunkgeom, 1, &E[p * dim]));  in ComputeFieldAtParticles_Mixed()
2427   PetscCall(PetscFEGeomRestoreChunk(ctx->fegeom, 0, 1, &chunkgeom));  in ComputeFieldAtParticles_Mixed()