working on scaling of delta

This commit is contained in:
James E McClure 2018-12-17 14:44:13 -05:00
parent 61d3eaf78e
commit 7cde2ca358

View File

@ -357,8 +357,8 @@ double MorphGrow(DoubleArray &BoundaryDist, DoubleArray &Dist, Array<char> &id,
if (rank == 0) printf(" delta=%f, growth=%f, max. displacement = %f \n",morph_delta, GrowthEstimate, MAX_DISPLACEMENT);
// Now adjust morph_delta
morph_delta *= TargetGrowth/GrowthEstimate;
MAX_DISPLACEMENT *= TargetGrowth/GrowthEstimate;
morph_delta *= max(TargetGrowth/GrowthEstimate,1.25);
MAX_DISPLACEMENT *= max(TargetGrowth/GrowthEstimate,1.25);
if (MAX_DISPLACEMENT > 2.0 ){
morph_delta /= 0.5*MAX_DISPLACEMENT;
//if (COUNT_FOR_LOOP > 2) COUNT_FOR_LOOP = 100;