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

tweak morph

Browse Source
This commit is contained in:
James E McClure
2018-12-17 16:07:41 -05:00
parent 133d650153
commit 8ae563b566
1 changed files with 3 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
analysis/morphology.cpp
View File

@@ -333,8 +333,9 @@ double MorphGrow(DoubleArray &BoundaryDist, DoubleArray &Dist, Array<char> &id,
double GrowthEstimate = 0.0; double GrowthEstimate = 0.0;
int COUNT_FOR_LOOP = 0; int COUNT_FOR_LOOP = 0;
double ERROR = 100.0;
if (rank == 0) printf("Estimate delta for growth=%f \n",TargetGrowth); if (rank == 0) printf("Estimate delta for growth=%f \n",TargetGrowth);
while ( fabs((GrowthEstimate - TargetGrowth)/TargetGrowth) < 0.01 && COUNT_FOR_LOOP < 10 ){ while ( ERROR > 0.01 && COUNT_FOR_LOOP < 10 ){
COUNT_FOR_LOOP++; COUNT_FOR_LOOP++;
count = 0.0; count = 0.0;
double MAX_DISPLACEMENT = 0.0; double MAX_DISPLACEMENT = 0.0;
@@ -354,6 +355,7 @@ double MorphGrow(DoubleArray &BoundaryDist, DoubleArray &Dist, Array<char> &id,
count=sumReduce( Dm->Comm, count); count=sumReduce( Dm->Comm, count);
MAX_DISPLACEMENT = maxReduce( Dm->Comm, MAX_DISPLACEMENT); MAX_DISPLACEMENT = maxReduce( Dm->Comm, MAX_DISPLACEMENT);
GrowthEstimate = count - count_original; GrowthEstimate = count - count_original;
ERROR = fabs((GrowthEstimate - TargetGrowth)/TargetGrowth);
if (rank == 0) printf(" delta=%f, growth=%f, max. displacement = %f \n",morph_delta, GrowthEstimate, MAX_DISPLACEMENT); if (rank == 0) printf(" delta=%f, growth=%f, max. displacement = %f \n",morph_delta, GrowthEstimate, MAX_DISPLACEMENT);
// Now adjust morph_delta // Now adjust morph_delta