From 637858558cd7722cbbdef6f46e9d59bacae7f255 Mon Sep 17 00:00:00 2001
From: James E McClure <mcclurej@vt.edu>
Date: Thu, 12 Nov 2015 19:59:40 -0500
Subject: [PATCH] pdated sign in distance calculation (no change in phase)

---
 common/Domain.h     | 9 ++++-----
 common/TwoPhase.cpp | 2 +-
 2 files changed, 5 insertions(+), 6 deletions(-)

diff --git a/common/Domain.h b/common/Domain.h
index b4b773f5..f977ed27 100755
--- a/common/Domain.h
+++ b/common/Domain.h
@@ -16,11 +16,8 @@
 #include "common/MPI_Helpers.h"
 #include "common/Communication.h"
 
-
 using namespace std;
 
-
-
 struct Domain{
     // Default constructor
 	Domain(int nx, int ny, int nz, int rnk, int npx, int npy, int npz, 
@@ -184,8 +181,10 @@ inline void SSO(DoubleArray &Distance, char *ID, Domain &Dm, int timesteps){
             for (j=1;j<Dm.Ny-1;j++){
                 for (i=1;i<Dm.Nx-1;i++){
 
-                    sign = Distance(i,j,k) / fabs(Distance(i,j,k));
-
+                	int n = k*Dm.Nx*Dm.Ny + j*Dm.Nx + i;
+//                    sign = Distance(i,j,k) / fabs(Distance(i,j,k));
+                	sign = -1;
+                	if (ID[n] == 0) sign = 1;
               /*
                     if (!(i+1<Nx))     nx=0.5*Distance(i,j,k);
                     else             nx=0.5*Distance(i+1,j,k);;
diff --git a/common/TwoPhase.cpp b/common/TwoPhase.cpp
index 9850b801..a9bae7b2 100644
--- a/common/TwoPhase.cpp
+++ b/common/TwoPhase.cpp
@@ -203,7 +203,7 @@ void TwoPhase::ColorToSignedDistance(double Beta, DoubleArray &ColorData, Double
 		}
 	}
 
-	SSO(DistData,Dm.id,Dm,40);
+	SSO(DistData,TempID,Dm,40);
 
     for (int k=0; k<Nz; k++){
 	  for (int j=0; j<Ny; j++){