tweak distance in solid

This commit is contained in:
James E McClure 2018-10-22 23:31:53 -04:00
parent bce4069111
commit 4933f84632

View File

@ -573,8 +573,12 @@ void ScaLBL_ColorModel::MorphInit(const double beta, const double morph_delta){
// temp -- distance based on analytical form McClure, Prins et al, Comp. Phys. Comm.
temp = -factor*log((1.0+value)/(1.0-value));
/// use this approximation close to the object
if (Averages->SDs(i,j,k) < 0.f ) phase_distance(i,j,k) = (-1.f)*Averages->SDs(i,j,k);
else if (fabs(phase_distance(i,j,k)) < 2.f ) phase_distance(i,j,k) = temp;
if (fabs(phase_distance(i,j,k)) < 2.f ){
if (Averages->SDs(i,j,k) < 0.f )
phase_distance(i,j,k) = (-1.f)*Averages->SDs(i,j,k);
else
phase_distance(i,j,k) = temp;
}
}
}
}