crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 21:25:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Add Material::volume() accessor function

Browse source
This commit is contained in:
Paul Romano 2019-07-03 23:08:08 -05:00
parent 1c97438be0
commit bbf529bef0
2 changed files with 18 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

4
include/openmc/material.h
View file

@ -80,6 +80,10 @@ public:
//! \param[in] units Units of density
void set_density(double density, gsl::cstring_span units);
//! Get volume of material
//! \return Volume in [cm^3]
double volume() const;
//! Write material data to HDF5
void to_hdf5(hid_t group) const;

22
src/material.cpp
View file

@ -888,6 +888,15 @@ void Material::set_density(double density, gsl::cstring_span units)
}
}
double Material::volume() const
{
if (volume_ < 0.0) {
throw std::runtime_error{"Volume for material with ID="
+ std::to_string(id_) + " not set."};
}
return volume_;
}
void Material::to_hdf5(hid_t group) const
{
hid_t material_group = create_group(group, "material " + std::to_string(id_));
@ -1250,16 +1259,13 @@ extern "C" int
openmc_material_get_volume(int32_t index, double* volume)
{
if (index >= 0 && index < model::materials.size()) {
auto& m = model::materials[index];
if (m->volume_ >= 0.0) {
*volume = m->volume_;
return 0;
} else {
std::stringstream msg;
msg << "Volume for material with ID=" << m->id_ << " not set.";
set_errmsg(msg);
try {
*volume = model::materials[index]->volume();
} catch (const std::exception& e) {
set_errmsg(e.what());
return OPENMC_E_UNASSIGNED;
}
return 0;
} else {
set_errmsg("Index in materials array is out of bounds.");
return OPENMC_E_OUT_OF_BOUNDS;