Merge pull request #2560 from devitocodes/mloubout/reduce_tiles

compiler: Fix reduction  multi_partile

Author: mloubout

devito/passes/clusters/blocking.py

@@ -508,7 +508,9 @@ def __init__(self, par_tile):
             # The user has supplied one specific par-tile per blocked nest
             self.umt_small = par_tile
             self.umt_sparse = par_tile
-            self.umt_reduce = par_tile
+            # Reduction tiles have to be unbounded.
+            self.umt_reduce = UnboundedMultiTuple(
+                *[UnboundTuple(*p, 1) for p in par_tile])
         else:
             # Special case 1: a smaller par-tile to avoid under-utilizing
             # computational resources when the iteration spaces are too small

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_MIN'), max=ValueLimit('FLT_MAX')),
-    np.float64: Bunch(min=-ValueLimit('DBL_MIN'), max=ValueLimit('DBL_MAX')),
+    np.float32: Bunch(min=-ValueLimit('FLT_MAX'), max=ValueLimit('FLT_MAX')),
+    np.float64: Bunch(min=-ValueLimit('DBL_MAX'), max=ValueLimit('DBL_MAX')),
 }