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

Only writing ID.dat (8bit)

Browse Source
This commit is contained in:
James E McClure 2015-11-12 21:31:13 -05:00
parent e869ef7151
commit 56d28abed0
1 changed files with 6 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
tests/ColorToBinary.cpp
View File

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