983   Vec        rhoRhs;      // Weak charge density, \int phi_i rho  in ComputeFieldAtParticles_Primal()  local
997   PetscCall(DMGetGlobalVector(dm, &rhoRhs));  in ComputeFieldAtParticles_Primal()
998   PetscCall(PetscObjectSetName((PetscObject)rhoRhs, "Weak charge density"));  in ComputeFieldAtParticles_Primal()
1007   PetscCall(MatMultTranspose(M_p, f, rhoRhs));  in ComputeFieldAtParticles_Primal()
1021     PetscCall(VecGetLocalSize(rhoRhs, &n));  in ComputeFieldAtParticles_Primal()
1022     PetscCall(VecGetArrayWrite(rhoRhs, &a));  in ComputeFieldAtParticles_Primal()
1029     PetscCall(VecRestoreArrayWrite(rhoRhs, &a));  in ComputeFieldAtParticles_Primal()
1036   PetscCall(KSPSolve(ksp, rhoRhs, rho));  in ComputeFieldAtParticles_Primal()
1038   PetscCall(VecScale(rhoRhs, -1.0));  in ComputeFieldAtParticles_Primal()
1040   PetscCall(VecViewFromOptions(rhoRhs, NULL, "-rho_view"));  in ComputeFieldAtParticles_Primal()
1046   PetscCall(SNESSolve(snes, rhoRhs, phi));  in ComputeFieldAtParticles_Primal()
1047   PetscCall(DMRestoreGlobalVector(dm, &rhoRhs));  in ComputeFieldAtParticles_Primal()