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

aprilnovak · Jun 12, 2019 · cd6ce5a · cd6ce5a
1 parent 110ba23
commit cd6ce5a
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_;