From e7c227bd9fdf1923c6ccff74d94614b0450924e2 Mon Sep 17 00:00:00 2001
From: Zhe Li <zzl109@gadi-login-04.gadi.nci.org.au>
Date: Wed, 13 Apr 2022 16:56:26 +1000
Subject: [PATCH] fix bugs in initializing electric potential; the Psi on solid
 was accidentally overwritten before.

---
 models/PoissonSolver.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/models/PoissonSolver.cpp b/models/PoissonSolver.cpp
index ab349b1d..5550ad02 100644
--- a/models/PoissonSolver.cpp
+++ b/models/PoissonSolver.cpp
@@ -437,7 +437,8 @@ void ScaLBL_Poisson::Potential_Init(double *psi_init){
                             idx = NLABELS;
                         }
                     }
-                    psi_init[n] = AFFINITY;
+                    int idx=Map(i,j,k);
+                    if (!(idx<0))  psi_init[n] = AFFINITY;
                 }
             }
         }