From 7dda5edc78e2269fc0d9e703b62730c5ebb0dc00 Mon Sep 17 00:00:00 2001
From: James E McClure <mcclurej@vt.edu>
Date: Thu, 12 Nov 2015 15:23:02 -0500
Subject: [PATCH] Adding error check to TestSegDist

---
 tests/TestSegDist.cpp | 11 ++++-------
 1 file changed, 4 insertions(+), 7 deletions(-)

diff --git a/tests/TestSegDist.cpp b/tests/TestSegDist.cpp
index 4dde2e01..c59d0b9f 100644
--- a/tests/TestSegDist.cpp
+++ b/tests/TestSegDist.cpp
@@ -71,14 +71,17 @@ int main(int argc, char **argv)
 	for (k=1;k<nz-1;k++){
 		for (j=1;j<ny-1;j++){
 			for (i=1;i<nx-1;i++){
+
+				// True signed distance
 				x = (nx-2)*Dm.iproc+i-1;
 				y = (ny-2)*Dm.jproc+j-1;
 				z = (nz-2)*Dm.kproc+k-1;
+				TrueDist(i,j,k) = sqrt((x-Cx)*(x-Cx)+(y-Cy)*(y-Cy)+(z-Cz)*(z-Cz)) - BubbleRadius;
 
 				n = k*nx*ny+j*nx+i;
 
 				// Initialize phase positions
-				if ( (x-Cx)*(x-Cx)+(y-Cy)*(y-Cy)+(z-Cz)*(z-Cz) < BubbleRadius*BubbleRadius){
+				if (TrueDist(i,j,k) < 0.0){
 					id[n] = 0;
 				}
 				else{
@@ -100,12 +103,6 @@ int main(int argc, char **argv)
 				n=k*nx*ny+j*nx+i;
 				// Initialize distance to +/- 1
 				Distance(i,j,k) = 2.0*id[n]-1.0;
-
-				// True signed distance
-				x = (nx-2)*Dm.iproc+i-1;
-				y = (ny-2)*Dm.jproc+j-1;
-				z = (nz-2)*Dm.kproc+k-1;
-				TrueDist(i,j,k) = sqrt((x-Cx)*(x-Cx)+(y-Cy)*(y-Cy)+(z-Cz)*(z-Cz)) - BubbleRadius;
 			}
 		}
 	}