From 43225f798022b3a0d67d20200bd0a7b0e1f6df68 Mon Sep 17 00:00:00 2001
From: Weiqun Zhang <weiqunzhang@lbl.gov>
Date: Mon, 26 Sep 2022 14:47:29 -0700
Subject: [PATCH] Fix MLEBNodeFDLaplacian bottom solver

MLEBNodeFDLaplacian is never singular because it has Dirichlet boundary on
the EB surface.  We did set the singular flag to false, but forgot about the
bottom solver used a different function to query.  This fixes it by
overriding the isBottomSingular function.
---
 Src/LinearSolvers/MLMG/AMReX_MLEBNodeFDLaplacian.H | 1 +
 1 file changed, 1 insertion(+)

diff --git a/Src/LinearSolvers/MLMG/AMReX_MLEBNodeFDLaplacian.H b/Src/LinearSolvers/MLMG/AMReX_MLEBNodeFDLaplacian.H
index 0d294c9da8c..41190f229a8 100644
--- a/Src/LinearSolvers/MLMG/AMReX_MLEBNodeFDLaplacian.H
+++ b/Src/LinearSolvers/MLMG/AMReX_MLEBNodeFDLaplacian.H
@@ -100,6 +100,7 @@ public:
     virtual void fixUpResidualMask (int amrlev, iMultiFab& resmsk) final override;
 
     virtual bool isSingular (int) const final override { return false; }
+    virtual bool isBottomSingular () const final override { return false; }
 
     virtual void compGrad (int amrlev, const Array<MultiFab*,AMREX_SPACEDIM>& grad,
                            MultiFab& sol, Location /*loc*/) const override;