bugfix in Jacobi - still not quite right

AMReX-Codes · Nov 27, 2024 · 4ed14fe · 4ed14fe
1 parent 1dd9b6f
commit 4ed14fe
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/ExampleCodes/LinearSolvers/GMRES/Poisson/AMReX_GMRES_Poisson.H b/ExampleCodes/LinearSolvers/GMRES/Poisson/AMReX_GMRES_Poisson.H
@@ -169,8 +169,8 @@ void GMRESPOISSON::precond (MultiFab& lhs, MultiFab const& rhs) const
 
         const GpuArray<Real, AMREX_SPACEDIM> dx = m_geom.CellSizeArray();
 
-        amrex::Real fac = 1.;
-        for (int d=0; d<AMREX_SPACEDIM; ++d) fac += 2./(dx[d]*dx[d]);
+        amrex::Real fac = 0.;
+        for (int d=0; d<AMREX_SPACEDIM; ++d) fac -= 2./(dx[d]*dx[d]);
 
         for ( MFIter mfi(lhs,TilingIfNotGPU()); mfi.isValid(); ++mfi ) {