From a26a241b696da612f7fb36eb7e883424782936ce Mon Sep 17 00:00:00 2001
From: James E McClure <mcclurej@vt.edu>
Date: Wed, 4 Feb 2015 13:01:10 -0500
Subject: [PATCH] slight implementation change to make pgi compilers on titna
 happy

---
 common/TwoPhase.h              | 3 ++-
 tests/lbpm_color_simulator.cpp | 6 +++---
 2 files changed, 5 insertions(+), 4 deletions(-)

diff --git a/common/TwoPhase.h b/common/TwoPhase.h
index f7cc9426..62a97223 100644
--- a/common/TwoPhase.h
+++ b/common/TwoPhase.h
@@ -8,7 +8,6 @@ class TwoPhase{
 	//...........................................................................
 	int n_nw_pts,n_ns_pts,n_ws_pts,n_nws_pts,n_local_sol_pts,n_local_nws_pts;
 	int n_nw_tris,n_ns_tris,n_ws_tris,n_nws_seg,n_local_sol_tris;
-	int cube[8][3] = {{0,0,0},{1,0,0},{0,1,0},{1,1,0},{0,0,1},{1,0,1},{0,1,1},{1,1,1}};
 	//...........................................................................
 	int nc;
 	int kstart,kfinish;
@@ -345,6 +344,8 @@ void TwoPhase::UpdateMeshValues(){
 void TwoPhase::ComputeLocal(){
 	int i,j,k,n;
 	double delphi;
+	int cube[8][3] = {{0,0,0},{1,0,0},{0,1,0},{1,1,0},{0,0,1},{1,0,1},{0,1,1},{1,1,1}};
+
 	for (int c=0;c<ncubes;c++){
 		// Get cube from the list
 		i = cubeList(0,c);
diff --git a/tests/lbpm_color_simulator.cpp b/tests/lbpm_color_simulator.cpp
index a920182a..9ad2311a 100644
--- a/tests/lbpm_color_simulator.cpp
+++ b/tests/lbpm_color_simulator.cpp
@@ -1547,7 +1547,7 @@ int main(int argc, char **argv)
 		// Timestep completed!
 		timestep++;
 		//...................................................................
-		if (timestep%1000 == 995){
+		if (timestep%5000 == 995){
 			//...........................................................................
 			// Copy the phase indicator field for the earlier timestep
 			DeviceBarrier();
@@ -1555,7 +1555,7 @@ int main(int argc, char **argv)
 			Averages.ColorToSignedDistance(beta,Averages.Phase.data,Averages.Phase_tplus.data);
 			//...........................................................................
 		}
-		if (timestep%1000 == 0){
+		if (timestep%5000 == 0){
 			//...........................................................................
 			// Copy the data for for the analysis timestep
 			//...........................................................................
@@ -1570,7 +1570,7 @@ int main(int argc, char **argv)
 			CopyToHost(Averages.Vel_z.data,&Velocity[2*N],N*sizeof(double));
 			MPI_Barrier(MPI_COMM_WORLD);
 		}
-		if (timestep%1000 == 5){
+		if (timestep%5000 == 5){
 			//...........................................................................
 			// Copy the phase indicator field for the later timestep
 			DeviceBarrier();