Add interface to material density method in OpenMC. Refs enrico-dev#17

aprilnovak · Jun 12, 2019 · 9251407 · 9251407
1 parent 31f4898
commit 9251407
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 0 deletions.
diff --git a/include/enrico/cell_instance.h b/include/enrico/cell_instance.h
@@ -34,6 +34,10 @@ class CellInstance {
  //! \return cell temperature in [K]
  double get_temperature() const;
 
+ //! Get the density of this cell
+ //! \return cell density in [g/cm^3]
+ double get_density() const;
+
  //! Check for equality
  bool operator==(const CellInstance& other) const;
 

diff --git a/src/cell_instance.cpp b/src/cell_instance.cpp
@@ -53,6 +53,13 @@ void CellInstance::set_density(double rho) const
  err_chk(openmc_material_set_density(material_index_, rho, "g/cm3"));
 }
 
+double CellInstance::get_density() const
+{
+ double rho;
+ err_chk(openmc_material_get_density(material_index));
+ return rho;
+}
+
 bool CellInstance::operator==(const CellInstance& other) const
 {
  return index_ == other.index_ && instance_ == other.instance_;