compiler: more robust safeinv

Author: mloubout

devito/passes/iet/misc.py

@@ -235,9 +235,9 @@ def _(expr, lang):
 
 @_lower_macro_math.register(SafeInv)
 def _(expr, lang):
-    eps = np.finfo(np.float32).resolution**2
+    eps = 2 * np.finfo(expr.base.dtype).resolution**2
     return (('SAFEINV(a, b)',
-             f'(((a) < {eps} || (b) < {eps}) ? (0.0F) : (1.0F / (a)))'),), {}
+             f'(((a + b) < {eps}F) ? (0.0F) : ((1.0F) / (a)))'),), {}
 
 
 @iet_pass

devito/symbolics/extended_dtypes.py

@@ -13,8 +13,8 @@
 limits_mapper = {
     np.int32: Bunch(min=ValueLimit('INT_MIN'), max=ValueLimit('INT_MAX')),
     np.int64: Bunch(min=ValueLimit('LONG_MIN'), max=ValueLimit('LONG_MAX')),
-    np.float32: Bunch(min=-ValueLimit('FLT_MAX'), max=ValueLimit('FLT_MAX')),
-    np.float64: Bunch(min=-ValueLimit('DBL_MAX'), max=ValueLimit('DBL_MAX')),
+    np.float32: Bunch(min=-ValueLimit('FLT_MIN'), max=ValueLimit('FLT_MAX')),
+    np.float64: Bunch(min=-ValueLimit('DBL_MIN'), max=ValueLimit('DBL_MAX')),
 }