Fix bug in MLEBABecLap getEBFluxes(). (#3217)

Bug introduced when adding anisotropy to MLEBABecLap, bareascaling goes on the denominator when calculating the EBfluxes.
AMReX-Codes · Mar 27, 2023 · 4fab19a · 4fab19a
1 parent 7aad6ca
commit 4fab19a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_2D_K.H b/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_2D_K.H
@@ -358,7 +358,7 @@ void mlebabeclap_ebflux (int i, int j, int k, int n,
             +       (            - gy*sy - gx*gy*sx*sy) * x(i ,jj,k,n)
             +       (                    + gx*gy*sx*sy) * x(ii,jj,k,n) ;
 
-        Real dphidn = (phib-phig)/dg * bareascaling;
+        Real dphidn = (phib-phig)/(dg * bareascaling);
         feb(i,j,k,n) = -beb(i,j,k,n) * dphidn;
     }
 }