From 4fab19a02359d05f011b5a5ee4aacfda4b86f640 Mon Sep 17 00:00:00 2001
From: Lucas Esclapez <13371051+esclapez@users.noreply.github.com>
Date: Mon, 27 Mar 2023 08:56:27 -0700
Subject: [PATCH] Fix bug in MLEBABecLap getEBFluxes(). (#3217)

Bug introduced when adding anisotropy to MLEBABecLap, bareascaling goes
on the denominator when calculating the EBfluxes.
---
 Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_2D_K.H | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_2D_K.H b/Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_2D_K.H
index 3369b8c4b34..f2bac77be77 100644
--- 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;
     }
 }