In MLEBABecLap::compGrad, add missing call to addInhomogNeumannFlux (A…

…MReX-Codes#3441) For inhomogeneous Neumann BC, we convert the inhomogeneous part into RHS and treat the BC as homogeneous Neumann during the solve. Thus, the flux is zero there before inhomogeneous flux is added.
IlianCS · Jul 23, 2023 · 6faf7ee · 6faf7ee
1 parent e41dbfa
commit 6faf7ee
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp b/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap.cpp
@@ -837,6 +837,8 @@ MLEBABecLap::compGrad (int amrlev, const Array<MultiFab*,AMREX_SPACEDIM>& grad,
             );
         }
     }
+
+    addInhomogNeumannFlux(amrlev, grad, sol, false);
 }
 
 void