We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 7bd7b32 commit 0c2e12aCopy full SHA for 0c2e12a
src/BaseStar.cpp
@@ -3910,7 +3910,7 @@ double BaseStar::DrawSNKickMagnitude(const double p_Sigma,
3910
kickMagnitude = 1001.0;
3911
double rand = p_Rand; // makes it possible to adjust if p_Rand is too low, to avoid getting stuck
3912
while( utils::Compare(kickMagnitude, 1000.0) > 0 ) { // maximum kick of 1000 km/s, following Disberg & Mandel (2025)
3913
- kickMagnitude = gsl_cdf_lognormal_Pinv(p_Rand, DISBERG_MANDEL_MU, DISBERG_MANDEL_SIGMA);
+ kickMagnitude = gsl_cdf_lognormal_Pinv(rand, DISBERG_MANDEL_MU, DISBERG_MANDEL_SIGMA);
3914
rand = min(rand + p_Rand + 0.0001, 1.0);
3915
}
3916
0 commit comments