research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[Transport] Fix segfault when species transport data is missing

Browse source 
Resolves #312.
This commit is contained in:
Ray Speth 2015-11-23 10:56:02 -05:00
parent a7db9a7411
commit bab07860e3
1 changed files with 17 additions and 12 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

29
src/transport/GasTransport.cpp
View file

@ -455,23 +455,28 @@ void GasTransport::getTransportData()
{
for (size_t k = 0; k < m_thermo->nSpecies(); k++) {
shared_ptr<Species> s = m_thermo->species(m_thermo->speciesName(k));
const GasTransportData& sptran =
dynamic_cast<GasTransportData&>(*s->transport.get());
if (sptran.geometry == "atom") {
const GasTransportData* sptran =
dynamic_cast<GasTransportData*>(s->transport.get());
if (!sptran) {
throw CanteraError("GasTransport::getTransportData",
"Missing gas-phase transport data for species '{}'.", s->name);
}
if (sptran->geometry == "atom") {
m_crot[k] = 0.0;
} else if (sptran.geometry == "linear") {
} else if (sptran->geometry == "linear") {
m_crot[k] = 1.0;
} else if (sptran.geometry == "nonlinear") {
} else if (sptran->geometry == "nonlinear") {
m_crot[k] = 1.5;
}
m_sigma[k] = sptran.diameter;
m_eps[k] = sptran.well_depth;
m_dipole(k,k) = sptran.dipole;
m_polar[k] = (sptran.dipole > 0);
m_alpha[k] = sptran.polarizability;
m_zrot[k] = sptran.rotational_relaxation;
m_w_ac[k] = sptran.acentric_factor;
m_sigma[k] = sptran->diameter;
m_eps[k] = sptran->well_depth;
m_dipole(k,k) = sptran->dipole;
m_polar[k] = (sptran->dipole > 0);
m_alpha[k] = sptran->polarizability;
m_zrot[k] = sptran->rotational_relaxation;
m_w_ac[k] = sptran->acentric_factor;
}
}