BUG: RosinRammler: correct the order of scale and shape parameters (fixes #3342)

Co-authored-by:Kutalmis Bercin <kutalmis.bercin@esi-group.com>
This commit is contained in:
Jan Gärtner
2025-05-09 13:14:53 +01:00
committed by Kutalmis Bercin
parent e5f147a211
commit c6220c162e

View File

@ -100,13 +100,13 @@ Foam::scalar Foam::distributionModels::RosinRammler::sample() const
Foam::scalar Foam::distributionModels::RosinRammler::meanValue() const
{
// (C:Eq. 5)
const scalar a = scalar(1)/lambda_ + scalar(1);
const scalar qMax = pow(maxValue_/n_, lambda_);
const scalar qMin = pow(minValue_/n_, lambda_);
const scalar a = scalar(1)/n_ + scalar(1);
const scalar qMax = pow(maxValue_/lambda_, n_);
const scalar qMin = pow(minValue_/lambda_, n_);
const scalar gMax = Math::incGamma_P(a, qMax);
const scalar gMin = Math::incGamma_P(a, qMin);
return n_/(exp(-qMin) - exp(-qMax))*(gMax - gMin);
return lambda_/(exp(-qMin) - exp(-qMax))*(gMax - gMin);
}