OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

pdated sign in distance calculation (no change in phase)

Browse Source
This commit is contained in:
James E McClure 2015-11-12 19:59:40 -05:00
parent c0d2c4dc3d
commit 637858558c
2 changed files with 5 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
common/Domain.h
View File

@ -16,11 +16,8 @@
#include "common/MPI_Helpers.h" #include "common/MPI_Helpers.h"
#include "common/Communication.h" #include "common/Communication.h"
using namespace std; using namespace std;
struct Domain{ struct Domain{
// Default constructor // Default constructor
Domain(int nx, int ny, int nz, int rnk, int npx, int npy, int npz, 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 (j=1;j<Dm.Ny-1;j++){
for (i=1;i<Dm.Nx-1;i++){ 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); if (!(i+1<Nx)) nx=0.5*Distance(i,j,k);
else nx=0.5*Distance(i+1,j,k);; else nx=0.5*Distance(i+1,j,k);;

2
common/TwoPhase.cpp
View File

@ -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 k=0; k<Nz; k++){
for (int j=0; j<Ny; j++){ for (int j=0; j<Ny; j++){