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

Update Color to Binary

Browse Source
This commit is contained in:
James E McClure
2015-11-12 20:40:05 -05:00
parent 637858558c
commit c4af8c2e47
1 changed files with 13 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
tests/ColorToBinary.cpp
View File

@@ -11,8 +11,9 @@
using namespace std;
inline void ReadFromRank(char *FILENAME, DoubleArray &Phase, int nx, int ny, int nz, int iproc, int jproc, int kproc)
inline void ReadFromRank(char *FILENAME, DoubleArray &Phase, DoubleArray &Pressure, DoubleArray &Vel_x,
DoubleArray &Vel_y, DoubleArray &Vel_z, int nx, int ny, int nz, int iproc, int
jproc, int kproc)
{
int i,j,k,q,n,N;
int iglobal,jglobal,kglobal;
@@ -99,10 +100,10 @@ inline void ReadFromRank(char *FILENAME, DoubleArray &Phase, int nx, int ny, int
kglobal = kproc*(nz-2)+k;
//........................................................................
Phase(iglobal,jglobal,kglobal) = (denA-denB)/(denA+denB);
//Pressure(iglobal,jglobal,kglobal) = value;
//Vel_x(iglobal,jglobal,kglobal) = vx;
//Vel_y(iglobal,jglobal,kglobal) = vy;
//Vel_z(iglobal,jglobal,kglobal) = vz;
Pressure(iglobal,jglobal,kglobal) = value;
Vel_x(iglobal,jglobal,kglobal) = vx;
Vel_y(iglobal,jglobal,kglobal) = vy;
Vel_z(iglobal,jglobal,kglobal) = vz;
//........................................................................
}
}
@@ -169,6 +170,10 @@ int main(int argc, char **argv)
printf("Full domain size: %i x %i x %i \n", Nx,Ny,Nz);
DoubleArray Phase(Nx,Ny,Nz);
DoubleArray Press(Nx,Ny,Nz);
DoubleArray Vel_x(Nx,Ny,Nz);
DoubleArray Vel_y(Nx,Ny,Nz);
DoubleArray Vel_z(Nx,Ny,Nz);
DoubleArray SignDist(Nx,Ny,Nz);
// Filenames used
@@ -208,8 +213,8 @@ int main(int argc, char **argv)
}
sprintf(LocalRankFilename,"%s%s","Restart.",LocalRankString);
ReadFromRank(LocalRankFilename,Phase,nx,ny,nz,iproc,jproc,kproc);
ReadFromRank(LocalRankFilename,Phase,Press,Vel_x,Vel_y,Vel_z,
nx,ny,nz,iproc,jproc,kproc);
}
}
}