void ThermoPhase::getActivities(doublereal* a) const
{
getActivityConcentrations(a);
for (size_t k = 0; k < nSpecies(); k++) {
a[k] /= standardConcentration(k);
}
}
void SurfPhase::getStandardVolumes(doublereal* vol) const
{
_updateThermo();
for (size_t k = 0; k < m_kk; k++) {
vol[k] = 1.0/standardConcentration(k);
}
}
void IdealMolalSoln::getActivityConcentrations(doublereal* c) const
{
if (m_formGC != 1) {
double c_solvent = standardConcentration();
getActivities(c);
for (size_t k = 0; k < m_kk; k++) {
c[k] *= c_solvent;
}
} else {
getActivities(c);
for (size_t k = 0; k < m_kk; k++) {
double c0 = standardConcentration(k);
c[k] *= c0;
}
}
}
doublereal ThermoPhase::logStandardConc(size_t k) const
{
return log(standardConcentration(k));
}
/*
* Returns the natural logarithm of the standard
* concentration of the kth species
*/
doublereal IdealSolnGasVPSS::logStandardConc(int k) const {
double c = standardConcentration(k);
double lc = std::log(c);
return lc;
}