We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 2b99c72 commit b159260Copy full SHA for b159260
src/library/HRGBase/IdealGasFunctions.cpp
@@ -1719,7 +1719,7 @@ namespace thermalfist {
1719
for (int i = 1; i <= order; ++i) {
1720
double k = static_cast<double>(i);
1721
if (m == 0.) {
1722
- ret += sign * k * 2. * T / k * (mu * mu - 4. * mu * T / k + 6. * T / k * T / k) / k * cfug;
+ ret += sign * 2. * T * (mu * mu - 4. * mu * T / k + 6. * T / k * T / k) / k * cfug;
1723
} else {
1724
ret += sign * k * (
1725
(m * (m * m + mu * mu - 4. * mu * T / k + 6. * T * T / k / k) * xMath::BesselKexp(0, i*moverT)
0 commit comments