From 172f66a16906f8b72380f0c60c58ce7cdfd5cd31 Mon Sep 17 00:00:00 2001
From: Rex Zhe Li <zhe.rex.li@gmail.com>
Date: Mon, 13 Sep 2021 04:45:37 -0400
Subject: [PATCH] PoissonSolver using MSE as stopping criterion should work now

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

diff --git a/models/PoissonSolver.cpp b/models/PoissonSolver.cpp
index a9513f25..3317aa95 100644
--- a/models/PoissonSolver.cpp
+++ b/models/PoissonSolver.cpp
@@ -150,6 +150,7 @@ void ScaLBL_Poisson::SetDomain(){
 	N = Nx*Ny*Nz;
 	Distance.resize(Nx,Ny,Nz);
 	Psi_host.resize(Nx,Ny,Nz);
+	Psi_previous.resize(Nx,Ny,Nz);
 
 	for (int i=0; i<Nx*Ny*Nz; i++) Dm->id[i] = 1;               // initialize this way
 	//Averages = std::shared_ptr<TwoPhase> ( new TwoPhase(Dm) ); // TwoPhase analysis object
@@ -558,7 +559,7 @@ void ScaLBL_Poisson::Run(double *ChargeDensity, int timestep_from_Study){
 		//************************************************************************/
 
         // Check convergence of steady-state solution
-        if (timesetp==2){
+        if (timestep==2){
             //save electric potential for convergence check
             ScaLBL_CopyToHost(Psi_previous.data(),Psi,sizeof(double)*Nx*Ny*Nz);
         }