From d91adc45a599fc127059e18702ae4aebd7703c8e Mon Sep 17 00:00:00 2001
From: Zeno Rogue <zeno@attnam.com>
Date: Mon, 19 Aug 2019 13:58:19 +0200
Subject: [PATCH] product:: in_underlying functions now simply call f when not
 in product

---
 nonisotropic.cpp | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/nonisotropic.cpp b/nonisotropic.cpp
index 286d16a3..bec722d7 100644
--- a/nonisotropic.cpp
+++ b/nonisotropic.cpp
@@ -719,11 +719,13 @@ EX namespace product {
     }
   
   EX void in_underlying_map(const reaction_t& f) {
-    ((hrmap_product*)currentmap)->in_underlying(f);
+    if(!prod) f();
+    else ((hrmap_product*)currentmap)->in_underlying(f);
     }
 
   #if HDR
   template<class T> auto in_underlying_geometry(const T& f) {
+    if(!prod) return f();
     dynamicval<eGeometry> g(geometry, underlying);
     dynamicval<geometry_information*> gc(cgip, underlying_cgip);
     return f();