multiply grad(phi) = E fluxes by -1 for the gmres case

ajnonaka · Dec 17, 2024 · cc9d724 · cc9d724
1 parent 98584df
commit cc9d724
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/Source/ablastr/fields/PoissonSolver.H b/Source/ablastr/fields/PoissonSolver.H
@@ -547,6 +547,9 @@ computePhi (
 #endif
                 );
             linop->compGrad(0, e_field[0], *phi[0], amrex::MLLinOp::Location::CellCenter);
+            for(int dir=0; dir<AMREX_SPACEDIM; ++dir) {
+                e_field[0][dir]->mult(-1.);
+            }
         }
 
         rho[lev]->mult(-ablastr::constant::SI::ep0);  // Multiply rho by epsilon again