Set the departure function in both ij and ji pairs

This commit is contained in:
Ian Bell
2017-11-17 21:24:20 -07:00
parent 3ef8042ce8
commit c065cacd8a

View File

@@ -315,6 +315,7 @@ double HelmholtzEOSMixtureBackend::get_binary_interaction_double(const std::size
void HelmholtzEOSMixtureBackend::set_binary_interaction_string(const std::size_t i, const std::size_t j, const std::string &parameter, const std::string & value){
if (parameter == "function"){
residual_helmholtz->Excess.DepartureFunctionMatrix[i][j].reset(get_departure_function(value));
residual_helmholtz->Excess.DepartureFunctionMatrix[j][i].reset(get_departure_function(value));
}
else{
throw ValueError(format("Cannot process this string parameter [%s] in set_binary_interaction_string", parameter.c_str()));