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

change morph_delta adapt

Browse Source
This commit is contained in:
James McClure 2018-11-21 14:25:06 -05:00
parent 4dedfca4b0
commit 5f1eb0ec29
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
models/ColorModel.cpp
View File

@ -635,13 +635,13 @@ void ScaLBL_ColorModel::Run(){
// update the volume
volA += delta_volume;
volB -= delta_volume;
//update size of morphological operation -- don't change sign of morph_delta
if (delta_volume_target / delta_volume > 0.f)
morph_delta *= 0.2*delta_volume_target / delta_volume;
if (morph_delta > 1.f) morph_delta = 1.f;
if (morph_delta < -1.f) morph_delta = -1.f;
if (fabs(morph_delta) < 0.05 ) morph_delta = 0.05*(morph_delta)/fabs(morph_delta); // set minimum
if (rank==0) printf(" Adjust morph delta: %f \n", morph_delta);
if (delta_volume_target / delta_volume > 0.f){
morph_delta *= 1.01*min((delta_volume_target - delta_volume) / delta_volume, 2.f);
if (morph_delta > 1.f) morph_delta = 1.f;
if (morph_delta < -1.f) morph_delta = -1.f;
if (fabs(morph_delta) < 0.05 ) morph_delta = 0.05*(morph_delta)/fabs(morph_delta); // set minimum
if (rank==0) printf(" Adjust morph delta: %f \n", morph_delta);
}
//MORPH_ADAPT = false;
if (volB/(volA + volB) > TARGET_SATURATION){
MORPH_ADAPT = false;