1668     for (PetscInt cp = 0; cp < Ncp; ++cp)  in ComputeFieldAtParticles_Primal()  local
1669       for (PetscInt d = 0; d < dim; ++d) pcoord[cp * dim + d] = coords[points[cp] * dim + d];  in ComputeFieldAtParticles_Primal()
1679     for (PetscInt cp = 0; cp < Ncp; ++cp) {  in ComputeFieldAtParticles_Primal()  local
1681       const PetscInt   p        = points[cp];  in ComputeFieldAtParticles_Primal()
1684 …EFreeInterpolateGradient_Static(fe, basisDer, clPhi, dim, chunkgeom->invJ, NULL, cp, &E[p * dim]));  in ComputeFieldAtParticles_Primal()
1794     for (PetscInt cp = 0; cp < Ncp; ++cp)  in ComputeFieldAtParticles_Mixed()  local
1795       for (PetscInt d = 0; d < dim; ++d) pcoord[cp * dim + d] = coords[points[cp] * dim + d];  in ComputeFieldAtParticles_Mixed()
1805     for (PetscInt cp = 0; cp < Ncp; ++cp) {  in ComputeFieldAtParticles_Mixed()  local
1806       const PetscInt p = points[cp];  in ComputeFieldAtParticles_Mixed()
1809       PetscCall(PetscFEInterpolateAtPoints_Static(fe, tab, clPhi, chunkgeom, cp, &E[p * dim]));  in ComputeFieldAtParticles_Mixed()
2068     for (PetscInt cp = 0; cp < Ncp; ++cp) {  in RHSObjectiveF()  local
2069       const PetscInt  p  = points[cp];  in RHSObjectiveF()