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

Little changes to lbpm_segmented_pp -- still need to test it

Browse Source
This commit is contained in:
James E McClure 2016-01-13 17:56:26 -05:00
parent 6cc6743ee6
commit 41c65125cc
1 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
tests/lbpm_segmented_pp.cpp
View File

@ -147,7 +147,7 @@ int main(int argc, char **argv)
MeanFilter(Averages.SDs); MeanFilter(Averages.SDs);
if (rank==0) printf("Initialized solid phase -- Converting to Signed Distance function \n"); if (rank==0) printf("Initialized solid phase -- Converting to Signed Distance function \n");
SSO(Averages.SDs,id,Dm,25); SSO(Averages.SDs,id,Dm,100);
sprintf(LocalRankFilename,"SignDist.%05i",rank); sprintf(LocalRankFilename,"SignDist.%05i",rank);
FILE *DIST = fopen(LocalRankFilename,"wb"); FILE *DIST = fopen(LocalRankFilename,"wb");
@ -178,7 +178,7 @@ int main(int argc, char **argv)
MeanFilter(Averages.Phase); MeanFilter(Averages.Phase);
if (rank==0) printf("Initialized non-wetting phase -- Converting to Signed Distance function \n"); if (rank==0) printf("Initialized non-wetting phase -- Converting to Signed Distance function \n");
SSO(Averages.Phase,id,Dm,25); SSO(Averages.Phase,id,Dm,100);
for (k=0;k<nz;k++){ for (k=0;k<nz;k++){
for (j=0;j<ny;j++){ for (j=0;j<ny;j++){
@ -217,17 +217,17 @@ int main(int argc, char **argv)
meshData[0].meshName = "domain"; meshData[0].meshName = "domain";
meshData[0].mesh = std::shared_ptr<IO::DomainMesh>( new IO::DomainMesh(Dm.rank_info,Nx-2,Ny-2,Nz-2,Lx,Ly,Lz) ); meshData[0].mesh = std::shared_ptr<IO::DomainMesh>( new IO::DomainMesh(Dm.rank_info,Nx-2,Ny-2,Nz-2,Lx,Ly,Lz) );
std::shared_ptr<IO::Variable> PhaseVar( new IO::Variable() ); std::shared_ptr<IO::Variable> PhaseVar( new IO::Variable() );
std::shared_ptr<IO::Variable> SignDistVar( new IO::Variable() ); std::shared_ptr<IO::Variable> SolidVar( new IO::Variable() );
std::shared_ptr<IO::Variable> BlobIDVar( new IO::Variable() ); std::shared_ptr<IO::Variable> BlobIDVar( new IO::Variable() );
PhaseVar->name = "phase"; PhaseVar->name = "Fluid";
PhaseVar->type = IO::VolumeVariable; PhaseVar->type = IO::VolumeVariable;
PhaseVar->dim = 1; PhaseVar->dim = 1;
PhaseVar->data.resize(Nx-2,Ny-2,Nz-2); PhaseVar->data.resize(Nx-2,Ny-2,Nz-2);
meshData[0].vars.push_back(PhaseVar); meshData[0].vars.push_back(PhaseVar);
SignDistVar->name = "SignDist"; SolidVar->name = "Solid";
SignDistVar->type = IO::VolumeVariable; SolidVar->type = IO::VolumeVariable;
SignDistVar->dim = 1; SolidVar->dim = 1;
SignDistVar->data.resize(Nx-2,Ny-2,Nz-2); SolidVar->data.resize(Nx-2,Ny-2,Nz-2);
meshData[0].vars.push_back(SignDistVar); meshData[0].vars.push_back(SignDistVar);
BlobIDVar->name = "BlobID"; BlobIDVar->name = "BlobID";
BlobIDVar->type = IO::VolumeVariable; BlobIDVar->type = IO::VolumeVariable;
@ -236,7 +236,7 @@ int main(int argc, char **argv)
meshData[0].vars.push_back(BlobIDVar); meshData[0].vars.push_back(BlobIDVar);
fillData.copy(Averages.SDn,PhaseVar->data); fillData.copy(Averages.SDn,PhaseVar->data);
fillData.copy(Averages.SDs,SignDistVar->data); fillData.copy(Averages.SDs,SolidVar->data);
fillData.copy(Averages.Label_NWP,BlobIDVar->data); fillData.copy(Averages.Label_NWP,BlobIDVar->data);
IO::writeData( 0, meshData, 2, comm ); IO::writeData( 0, meshData, 2, comm );