| 235 | } |
| 236 | |
| 237 | double MultiPhase::phaseCharge(size_t p) const |
| 238 | { |
| 239 | double phasesum = 0.0; |
| 240 | size_t nsp = m_phase[p]->nSpecies(); |
| 241 | for (size_t ik = 0; ik < nsp; ik++) { |
| 242 | size_t k = speciesIndex(ik, p); |
| 243 | phasesum += m_phase[p]->charge(ik)*m_moleFractions[k]; |
| 244 | } |
| 245 | return Faraday*phasesum*m_moles[p]; |
| 246 | } |
| 247 | |
| 248 | void MultiPhase::getChemPotentials(double* mu) const |
| 249 | { |
no test coverage detected