Ejemplo n.º 1
0
doublereal MaskellSolidSolnPhase::entropy_mole() const
{
    _updateThermo();
    const doublereal s0 = GasConstant * mean_X(m_s0_R);
    const doublereal r = moleFraction(product_species_index);
    const doublereal fmval = fm(r);
    const doublereal rfm = r * fmval;
    return s0 + GasConstant * (xlogx(1-rfm) - xlogx(rfm) - xlogx(1-r-rfm) - xlogx((1-fmval)*r) - xlogx(1-r) - xlogx(r));
}
Ejemplo n.º 2
0
double Utils::info(int_array &counts)
{

    int total = 0;
    double x = 0;
    for (auto count : counts) {
        x -= xlogx(count);
        total += count;
    }
    return x + xlogx(total);
}