339     PetscInt        Ncp, cp;  in ComputeFieldAtParticles_Primal()  local
344     for (cp = 0; cp < Ncp; ++cp)  in ComputeFieldAtParticles_Primal()
345       for (d = 0; d < dim; ++d) pcoord[cp * dim + d] = coords[points[cp] * dim + d];  in ComputeFieldAtParticles_Primal()
350     for (cp = 0; cp < Ncp; ++cp) {  in ComputeFieldAtParticles_Primal()
352       const PetscInt   p        = points[cp];  in ComputeFieldAtParticles_Primal()
355 …PetscCall(PetscFEFreeInterpolateGradient_Static(fe, basisDer, clPhi, dim, invJ, NULL, cp, &E[p * d…  in ComputeFieldAtParticles_Primal()
456     PetscInt        Ncp, cp;  in ComputeFieldAtParticles_Mixed()  local
461     for (cp = 0; cp < Ncp; ++cp)  in ComputeFieldAtParticles_Mixed()
462       for (d = 0; d < dim; ++d) pcoord[cp * dim + d] = coords[points[cp] * dim + d];  in ComputeFieldAtParticles_Mixed()
467     for (cp = 0; cp < Ncp; ++cp) {  in ComputeFieldAtParticles_Mixed()
468       const PetscInt p = points[cp];  in ComputeFieldAtParticles_Mixed()
473       PetscCall(PetscFEInterpolateAtPoints_Static(fe, tab, clPhi, &feGeometry, cp, &E[p * dim]));  in ComputeFieldAtParticles_Mixed()
717     for (PetscInt cp = 0; cp < Ncp; ++cp) {  in RHSObjectiveF()  local
718       const PetscInt  p     = points[cp];  in RHSObjectiveF()