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

fix distance

Browse Source
This commit is contained in:
JamesEMcclure
2019-10-25 06:10:55 -04:00
parent f3f08722d3
commit 90fe77ac6a
1 changed files with 2 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
tests/TestPoiseuille.cpp
View File

@@ -34,15 +34,14 @@ void ParallelPlates(ScaLBL_MRTModel &MRT){
for (int k=0;k<Nz;k++){
for (int j=0;j<Ny;j++){
for (int i=0;i<Nx;i++){
int n=k*Nx*Ny+j*Nx+i;
// Initialize distance to +/- 1
MRT.Distance(i,j,k) = 2.0*double(id_solid(i,j,k))-1.0;
}
}
}
if (rank==0) printf("Initialized solid phase -- Converting to Signed Distance function \n");
if (MRT.rank==0) printf("Initialized solid phase -- Converting to Signed Distance function \n");
CalcDist(MRT.Distance,id_solid,*MRT.Dm);
if (rank == 0) cout << "Domain set." << endl;
if (MRT.rank == 0) cout << "Domain set." << endl;
}
//***************************************************************************************