diff --git a/src/thermo/Phase.cpp b/src/thermo/Phase.cpp index a39cea996..85798d5ec 100644 --- a/src/thermo/Phase.cpp +++ b/src/thermo/Phase.cpp @@ -657,7 +657,7 @@ doublereal Phase::sum_xlogx() const { double sumxlogx = 0; for (size_t k = 0; k < m_kk; k++) { - sumxlogx += m_ym[k] * std::log(m_ym[k] + Tiny); + sumxlogx += m_ym[k] * std::log(std::max(m_ym[k], SmallNumber)); } return m_mmw * sumxlogx + std::log(m_mmw); }