1567   Vec        rhoRhs;      // Weak charge density, \int phi_i rho  in ComputeFieldAtParticles_Primal()  local
1581   PetscCall(DMGetGlobalVector(dm, &rhoRhs));  in ComputeFieldAtParticles_Primal()
1582   PetscCall(PetscObjectSetName((PetscObject)rhoRhs, "Weak charge density"));  in ComputeFieldAtParticles_Primal()
1591   PetscCall(MatMultTranspose(M_p, f, rhoRhs));  in ComputeFieldAtParticles_Primal()
1606     PetscCall(VecGetLocalSize(rhoRhs, &n));  in ComputeFieldAtParticles_Primal()
1607     PetscCall(VecGetArrayWrite(rhoRhs, &a));  in ComputeFieldAtParticles_Primal()
1614     PetscCall(VecRestoreArrayWrite(rhoRhs, &a));  in ComputeFieldAtParticles_Primal()
1621   PetscCall(KSPSolve(ksp, rhoRhs, rho));  in ComputeFieldAtParticles_Primal()
1623   PetscCall(VecScale(rhoRhs, -1.0));  in ComputeFieldAtParticles_Primal()
1625   PetscCall(VecViewFromOptions(rhoRhs, NULL, "-rho_view"));  in ComputeFieldAtParticles_Primal()
1631   PetscCall(SNESSolve(snes, rhoRhs, phi));  in ComputeFieldAtParticles_Primal()
1632   PetscCall(DMRestoreGlobalVector(dm, &rhoRhs));  in ComputeFieldAtParticles_Primal()
1708 …Vec        rhoRhs, rhoRhsFull;   // Weak charge density, \int phi_i rho, and embedding in mixed pr…  in ComputeFieldAtParticles_Mixed()  local
1721   PetscCall(DMGetGlobalVector(ctx->dmPot, &rhoRhs));  in ComputeFieldAtParticles_Mixed()
1722   PetscCall(PetscObjectSetName((PetscObject)rhoRhs, "Weak charge density"));  in ComputeFieldAtParticles_Mixed()
1733   PetscCall(MatMultTranspose(M_p, f, rhoRhs));  in ComputeFieldAtParticles_Mixed()
1740   PetscCall(KSPSolve(ksp, rhoRhs, rho));  in ComputeFieldAtParticles_Mixed()
1742   PetscCall(VecISCopy(rhoRhsFull, ctx->isPot, SCATTER_FORWARD, rhoRhs));  in ComputeFieldAtParticles_Mixed()
1745   PetscCall(VecViewFromOptions(rhoRhs, NULL, "-rho_view"));  in ComputeFieldAtParticles_Mixed()
1748   PetscCall(DMRestoreGlobalVector(ctx->dmPot, &rhoRhs));  in ComputeFieldAtParticles_Mixed()