throw if requested storage region is empty

Author: philippwindischhofer

eisvogel/core/Geometry.hh

@@ -13,6 +13,8 @@ public:
   CylinderRegion();
   CylinderRegion(scalar_t r_min, scalar_t r_max, scalar_t z_min, scalar_t z_max);
   void SetRegion(scalar_t r_min, scalar_t r_max, scalar_t z_min, scalar_t z_max);
+
+  bool IsEmpty();
 
   scalar_t GetRMin() const {return m_r_min;};
   scalar_t GetRMax() const {return m_r_max;};

eisvogel/core/impl/Geometry.hxx

@@ -26,6 +26,10 @@ void CylinderRegion::throw_if_unphysical() {
   }
 }
 
+bool CylinderRegion::IsEmpty() {
+  return ((m_r_min == m_r_max) || (m_z_min == m_z_max));
+}
+
 CylinderGeometry::CylinderGeometry(scalar_t r_max, scalar_t z_min, scalar_t z_max,
 				   std::function<scalar_t(scalar_t r, scalar_t z)> eps_func) :
   m_r_max(r_max), m_z_min(z_min), m_z_max(z_max), m_eps_func(eps_func) {

eisvogel/meep/impl/MEEPCylindricalGreensFunctionCalculator.hxx

@@ -613,9 +613,13 @@ namespace GreensFunctionCalculator::MEEP {
     m_t_end(t_end), m_geom(geom), m_request_to_store(request_to_store), m_antenna(antenna) {
 
     // Sanity checks
-    if(!m_geom.contains(request_to_store)) {
+    if(!m_geom.contains(m_request_to_store)) {
       throw std::runtime_error("Error: requested storage region contains undefined domain!");
     }
+
+    if(m_request_to_store.IsEmpty()) {
+      throw std::runtime_error("Error: requested storage region is empty!");
+    }
   }
 
   void CylindricalGreensFunctionCalculator::merge_mpi_chunks(std::filesystem::path outdir, const std::vector<std::filesystem::path>& indirs) {