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

tweak morphlogy routine

Browse Source
This commit is contained in:
James E McClure 2018-12-17 13:00:56 -05:00
parent 8cf445c787
commit 7b575c4d81
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
analysis/morphology.cpp
View File

@ -341,9 +341,8 @@ double MorphGrow(DoubleArray &BoundaryDist, DoubleArray &Dist, Array<char> &id,
for (int k=0; k<Nz; k++){
for (int j=0; j<Ny; j++){
for (int i=0; i<Nx; i++){
//double walldist=BoundaryDist(i,j,k);
//double wallweight = 1.f / (1+exp(-5.f*(walldist-1.f)));
double wallweight = 1.f;
double walldist=BoundaryDist(i,j,k);
double wallweight = 1.f / (1+exp(-5.f*(walldist-1.f)));
if (fabs(wallweight*morph_delta) > MAX_DISPLACEMENT) MAX_DISPLACEMENT= fabs(wallweight*morph_delta);
if (Dist(i,j,k) - wallweight*morph_delta > 0.0){