From b5572d3f7f1774cc1e56da73c3c2836f25a93b9c Mon Sep 17 00:00:00 2001
From: Hamidreza Norouzi <hamid.r.norouzi@gmail.com>
Date: Mon, 3 Apr 2023 05:48:54 -0700
Subject: [PATCH] bug fix for push_back

---
 DEMSystems/domainDistribute/domainDistribute.cpp    | 4 ++--
 src/phasicFlow/containers/VectorHD/VectorDual.hpp   | 6 +++++-
 src/phasicFlow/containers/VectorHD/VectorSingle.hpp | 5 ++++-
 3 files changed, 11 insertions(+), 4 deletions(-)

diff --git a/DEMSystems/domainDistribute/domainDistribute.cpp b/DEMSystems/domainDistribute/domainDistribute.cpp
index b888f2d6..b6985243 100644
--- a/DEMSystems/domainDistribute/domainDistribute.cpp
+++ b/DEMSystems/domainDistribute/domainDistribute.cpp
@@ -52,11 +52,13 @@ bool pFlow::domainDistribute::locateParticles(
 	
 	range activeRange = mask.activeRange();
 
+
 	for(int32 di=0; di<numDomains_; di++)
 	{
 		particlesInDomains_[di].clear();
 	}
 
+	
 	for(int32 i=activeRange.first; i<activeRange.second; i++)
 	{
 		if(mask(i))
@@ -76,8 +78,6 @@ bool pFlow::domainDistribute::locateParticles(
 		numParInDomain_[di] = particlesInDomains_[di].size();
 	}
 
-	output<<" numParInDomain_ "<< numParInDomain_<<endl;
-
 	return true;
 }
 
diff --git a/src/phasicFlow/containers/VectorHD/VectorDual.hpp b/src/phasicFlow/containers/VectorHD/VectorDual.hpp
index c4f29fe1..0839601d 100644
--- a/src/phasicFlow/containers/VectorHD/VectorDual.hpp
+++ b/src/phasicFlow/containers/VectorHD/VectorDual.hpp
@@ -744,7 +744,11 @@ public:
 			syncToHost();
 			modifyOnHost();
 			
-			if(size_ == capacity_) changeSize(capacity_);
+			if(size_ == capacity_) 
+			{
+				changeSize(evalCapacity(capacity_), true);
+
+			}
 			data()[size_++] = val;
 			subViewsUpdated_ = false;
 			
diff --git a/src/phasicFlow/containers/VectorHD/VectorSingle.hpp b/src/phasicFlow/containers/VectorHD/VectorSingle.hpp
index f2e35fa6..cbe173aa 100644
--- a/src/phasicFlow/containers/VectorHD/VectorSingle.hpp
+++ b/src/phasicFlow/containers/VectorHD/VectorSingle.hpp
@@ -710,7 +710,10 @@ public:
 			void>::type
 		push_back(const T& val)
 		{
-			if(size_ == capacity_) changeSize(capacity_);
+			if(size_ == capacity_) 
+			{
+				changeSize(evalCapacity(capacity_), true);
+			}
 			data()[size_++] = val;
 			subViewUpdated_ = false;
 		}