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

updated morph condition

Browse Source
This commit is contained in:
James E McClure 2018-12-23 14:13:15 -05:00
parent 84801a625a
commit d17d154c9c
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
models/ColorModel.cpp
View File

@ -653,11 +653,11 @@ void ScaLBL_ColorModel::Run(){
if (MORPH_ADAPT ){ if (MORPH_ADAPT ){
if (rank==0) printf("***Morphological step with target volume change %f ***\n", delta_volume_target); if (rank==0) printf("***Morphological step with target volume change %f ***\n", delta_volume_target);
//double delta_volume_target = volB - (volA + volB)*TARGET_SATURATION; // change in volume to A //double delta_volume_target = volB - (volA + volB)*TARGET_SATURATION; // change in volume to A
delta_volume = MorphInit(beta,delta_volume_target); delta_volume += MorphInit(beta,delta_volume_target-delta_volume);
//if ( (delta_volume_target - delta_volume )/delta_volume_target < 0.1 ){ if ( (delta_volume - delta_volume_target)/delta_volume_target > 0.0 ){
MORPH_ADAPT = false; MORPH_ADAPT = false;
delta_volume = 0.0; delta_volume = 0.0;
//} }
/*if ((delta_volume_target - delta_volume) / delta_volume > 0.f){ /*if ((delta_volume_target - delta_volume) / delta_volume > 0.f){
morph_delta *= 1.01*min((delta_volume_target - delta_volume) / delta_volume, 2.0); morph_delta *= 1.01*min((delta_volume_target - delta_volume) / delta_volume, 2.0);
if (morph_delta > 1.f) morph_delta = 1.f; if (morph_delta > 1.f) morph_delta = 1.f;