2161     PetscReal *pcoord = &coords[p * dim];  in ComputeFieldAtParticles_Coulomb()  local
2171       for (PetscInt d = 0; d < dim; ++d) rpq[d] = pcoord[d] - qcoord[d];  in ComputeFieldAtParticles_Coulomb()
2250   PetscReal      *pcoord, *refcoord;  in ComputeFieldAtParticles_Primal()  local
2261   PetscCall(DMGetWorkArray(dm, maxNcp * dim, MPIU_REAL, &pcoord));  in ComputeFieldAtParticles_Primal()
2273       for (PetscInt d = 0; d < dim; ++d) pcoord[cp * dim + d] = coords[points[cp] * dim + d];  in ComputeFieldAtParticles_Primal()
2278 …CoordinatesRealToRef(dim, dim, x0, chunkgeom->v, chunkgeom->invJ, &pcoord[dim * i], &refcoord[dim …  in ComputeFieldAtParticles_Primal()
2294   PetscCall(DMRestoreWorkArray(dm, maxNcp * dim, MPIU_REAL, &pcoord));  in ComputeFieldAtParticles_Primal()
2380   PetscReal      *pcoord, *refcoord;  in ComputeFieldAtParticles_Mixed()  local
2391   PetscCall(DMGetWorkArray(dm, maxNcp * dim, MPIU_REAL, &pcoord));  in ComputeFieldAtParticles_Mixed()
2400       for (PetscInt d = 0; d < dim; ++d) pcoord[cp * dim + d] = coords[points[cp] * dim + d];  in ComputeFieldAtParticles_Mixed()
2405 …CoordinatesRealToRef(dim, dim, x0, chunkgeom->v, chunkgeom->invJ, &pcoord[dim * i], &refcoord[dim …  in ComputeFieldAtParticles_Mixed()
2421   PetscCall(DMRestoreWorkArray(dm, maxNcp * dim, MPIU_REAL, &pcoord));  in ComputeFieldAtParticles_Mixed()