From ee79d6a1ffb3f3c1cd4ea536d3ed88d5486852b8 Mon Sep 17 00:00:00 2001
From: James E McClure <jemcclur@gmail.com>
Date: Sat, 16 Feb 2019 11:07:57 -0500
Subject: [PATCH] adding checkerboard inlet option

---
 tests/lbpm_serial_decomp.cpp | 28 +++++++++++++++++++++++++++-
 1 file changed, 27 insertions(+), 1 deletion(-)

diff --git a/tests/lbpm_serial_decomp.cpp b/tests/lbpm_serial_decomp.cpp
index da016fb1..b284ecfd 100644
--- a/tests/lbpm_serial_decomp.cpp
+++ b/tests/lbpm_serial_decomp.cpp
@@ -53,7 +53,9 @@ int main(int argc, char **argv)
 	int64_t i,j,k,n;
 	int BC=0;
 	int64_t xStart,yStart,zStart;
-	//  char fluidValue,solidValue;
+	int64_t checkerSize;
+	int64_t CropInlet[3];
+	//  char fluidValue,solidValue;	
 
 	xStart=yStart=zStart=0;
 	// read the input database 
@@ -72,6 +74,12 @@ int main(int argc, char **argv)
 		yStart = offset[1];
 		zStart = offset[2];
 	}
+	if (domain_db->keyExists( "CropInlet" )){
+		CropInlet = domain_db->getVector<int>( "CropInlet" );
+	}
+	if (domain_db->keyExists( "checkerSize" )){
+		checkerSize = domain_db->getScalar<int64_t>( "checkerSize" );
+	}
 	auto ReadValues = domain_db->getVector<char>( "ReadValues" );
 	auto WriteValues = domain_db->getVector<char>( "WriteValues" );
 	auto ReadType = domain_db->getScalar<std::string>( "ReadType" );
@@ -136,6 +144,24 @@ int main(int argc, char **argv)
 		printf("Read segmented data from %s \n",Filename.c_str());
 	}
 
+	if (CropInlet[0] > 0){
+		// use checkerboard pattern
+		for (int k = 0; k<SIZE[2]; k++){
+			for (int j = 0; j<SIZE[1]; j++){
+				for (int i = 0; i<CropInlet[0]; i++){
+					if ( ((j+k)/checkerSize)%2 == 0){
+						// solid checkers
+						SegData[k*SIZE[0]*SIZE[1]+i*SIZE[0]+i] = 0;
+					}
+					else{
+						// void checkers
+						SegData[k*SIZE[0]*SIZE[1]+i*SIZE[0]+i] = 2;
+					}
+				}
+			}
+		}
+	}
+
 	// Get the rank info
 	int64_t N = (nx+2)*(ny+2)*(nz+2);