2191   Vec        rhoRhs;      // Weak charge density, \int phi_i rho  in ComputeFieldAtParticles_Primal()  local
2205   PetscCall(DMGetGlobalVector(dm, &rhoRhs));  in ComputeFieldAtParticles_Primal()
2206   PetscCall(PetscObjectSetName((PetscObject)rhoRhs, "Weak charge density"));  in ComputeFieldAtParticles_Primal()
2215   PetscCall(MatMultTranspose(M_p, f, rhoRhs));  in ComputeFieldAtParticles_Primal()
2222   PetscCall(KSPSolve(ksp, rhoRhs, rho));  in ComputeFieldAtParticles_Primal()
2224   PetscCall(VecScale(rhoRhs, -1.0));  in ComputeFieldAtParticles_Primal()
2226   PetscCall(VecViewFromOptions(rhoRhs, NULL, "-rho_view"));  in ComputeFieldAtParticles_Primal()
2232   PetscCall(SNESSolve(snes, rhoRhs, phi));  in ComputeFieldAtParticles_Primal()
2233   PetscCall(DMRestoreGlobalVector(dm, &rhoRhs));  in ComputeFieldAtParticles_Primal()
2312 …Vec        rhoRhs, rhoRhsFull;   // Weak charge density, \int phi_i rho, and embedding in mixed pr…  in ComputeFieldAtParticles_Mixed()  local
2326   PetscCall(DMGetGlobalVector(ctx->dmPot, &rhoRhs));  in ComputeFieldAtParticles_Mixed()
2327   PetscCall(PetscObjectSetName((PetscObject)rhoRhs, "Weak charge density"));  in ComputeFieldAtParticles_Mixed()
2338   PetscCall(MatMultTranspose(M_p, f, rhoRhs));  in ComputeFieldAtParticles_Mixed()
2345   PetscCall(KSPSolve(ksp, rhoRhs, rho));  in ComputeFieldAtParticles_Mixed()
2347   PetscCall(VecISCopy(rhoRhsFull, ctx->isPot, SCATTER_FORWARD, rhoRhs));  in ComputeFieldAtParticles_Mixed()
2350   PetscCall(VecViewFromOptions(rhoRhs, NULL, "-rho_view"));  in ComputeFieldAtParticles_Mixed()
2353   PetscCall(DMRestoreGlobalVector(ctx->dmPot, &rhoRhs));  in ComputeFieldAtParticles_Mixed()