eebasso
diff --git a/‎Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_2D_K.H‎
Lines changed: 2 additions & 0 deletions b/‎Src/LinearSolvers/MLMG/AMReX_MLEBABecLap_2D_K.H‎
Lines changed: 2 additions & 0 deletions
@@ -153,6 +153,7 @@ void mlebabeclap_adotx_centroid (Box const& box, Array4<Real> const& y,
                                                          on_y_face, domlo_y, domhi_y);
 
                 feb *= ba(i,j,k) * beb(i,j,k,n);
+                feb *= Real(0.0);
             }
 
 
@@ -291,6 +292,7 @@ void mlebabeclap_adotx (Box const& box, Array4<Real> const& y,
                 Real dphidn = (phib-phig) / dg;
 
                 feb = dphidn * ba(i,j,k) * beb(i,j,k,n);
+                feb *= Real(0.0);
             }
 
             y(i,j,k,n) = alpha*a(i,j,k)*x(i,j,k,n) + (Real(1.0)/kappa) *
Original file line number	Diff line number	Diff line change
`@@ -153,6 +153,7 @@ void mlebabeclap_adotx_centroid (Box const& box, Array4<Real> const& y,`
`153`	`153`	`on_y_face, domlo_y, domhi_y);`
`154`	`154`
`155`	`155`	`feb = ba(i,j,k) beb(i,j,k,n);`
	`156`	`+ feb *= Real(0.0);`
`156`	`157`	`}`
`157`	`158`
`158`	`159`
`@@ -291,6 +292,7 @@ void mlebabeclap_adotx (Box const& box, Array4<Real> const& y,`
`291`	`292`	`Real dphidn = (phib-phig) / dg;`
`292`	`293`
`293`	`294`	`feb = dphidn * ba(i,j,k) * beb(i,j,k,n);`
	`295`	`+ feb *= Real(0.0);`
`294`	`296`	`}`
`295`	`297`
`296`	`298`	`y(i,j,k,n) = alphaa(i,j,k)x(i,j,k,n) + (Real(1.0)/kappa) *`