Remove dtypes lowering from IET layer

devito/passes/iet/definitions.py

@@ -13,7 +13,7 @@
                        FindNodes, FindSymbols, MapExprStmts, Transformer,
                        make_callable)
 from devito.passes import is_gpu_create
-from devito.passes.iet.dtypes import lower_dtypes
+from devito.passes.iet.dtypes import include_complex
 from devito.passes.iet.engine import iet_pass
 from devito.passes.iet.langbase import LangBB
 from devito.symbolics import (Byref, DefFunction, FieldFromPointer, IndexedPointer,
@@ -466,8 +466,8 @@ def place_casts(self, iet, **kwargs):
         return iet, {}
 
     @iet_pass
-    def make_langtypes(self, iet):
-        iet, metadata = lower_dtypes(iet, self.lang, self.compiler, self.sregistry)
+    def include_complex(self, iet):
+        iet, metadata = include_complex(iet, self.lang, self.compiler)
         return iet, metadata
 
     def process(self, graph):
@@ -476,7 +476,7 @@ def process(self, graph):
         """
         self.place_definitions(graph, globs=set())
         self.place_casts(graph)
-        self.make_langtypes(graph)
+        self.include_complex(graph)
 
 
 class DeviceAwareDataManager(DataManager):
@@ -618,7 +618,7 @@ def process(self, graph):
         self.place_devptr(graph)
         self.place_bundling(graph, writes_input=graph.writes_input)
         self.place_casts(graph)
-        self.make_langtypes(graph)
+        self.include_complex(graph)
 
 
 def make_zero_init(obj, rcompile, sregistry):

devito/passes/iet/dtypes.py

@@ -1,63 +1,21 @@
 import numpy as np
 import ctypes
 
-from devito.ir import FindSymbols, Uxreplace
-from devito.ir.iet.nodes import Dereference
-from devito.tools.utils import as_list
-from devito.types.basic import Symbol
+from devito.ir import FindSymbols
 
-__all__ = ['lower_dtypes']
+__all__ = ['include_complex']
 
 
-def lower_dtypes(iet, lang, compiler, sregistry):
+def include_complex(iet, lang, compiler):
     """
-    Lower language-specific dtypes and add headers for complex arithmetic
-    """
-    # Include complex headers if needed (before we replace complex dtypes)
-    metadata = _complex_includes(iet, lang, compiler)
-
-    body_prefix = []  # Derefs to prepend to the body
-    body_mapper = {}
-    params_mapper = {}
-
-    # Lower scalar float16s to pointers and dereference them
-    if lang.get('half_types') is not None:
-        half, half_p = lang['half_types']  # dtype mappings for half float
-
-        for s in FindSymbols('scalars').visit(iet):
-            if s.dtype != np.float16 or s not in iet.parameters:
-                continue
-
-            ptr = s._rebuild(dtype=half_p, is_const=True)
-            val = Symbol(name=sregistry.make_name(prefix='hf'), dtype=half,
-                         is_const=s.is_const)
-
-            params_mapper[s], body_mapper[s] = ptr, val
-            body_prefix.append(Dereference(val, ptr))  # val = *ptr
-
-    # Lower remaining language-specific dtypes
-    for s in FindSymbols('indexeds|basics|symbolics').visit(iet):
-        if s.dtype in lang['types'] and s not in params_mapper:
-            body_mapper[s] = params_mapper[s] = s._rebuild(dtype=lang['types'][s.dtype])
-
-    # Apply the dtype replacements
-    body = body_prefix + as_list(Uxreplace(body_mapper).visit(iet.body))
-    params = Uxreplace(params_mapper).visit(iet.parameters)
-
-    iet = iet._rebuild(body=body, parameters=params)
-    return iet, metadata
-
-
-def _complex_includes(iet, lang, compiler):
-    """
-    Add headers for complex arithmetic
+    Include complex arithmetic headers for the given language, if needed.
     """
     # Check if there is complex numbers that always take dtype precedence
     types = {f.dtype for f in FindSymbols().visit(iet)
              if not issubclass(f.dtype, ctypes._Pointer)}
 
     if not any(np.issubdtype(d, np.complexfloating) for d in types):
-        return {}
+        return iet, {}
 
     metadata = {}
     lib = (lang['header-complex'],)
@@ -75,4 +33,4 @@ def _complex_includes(iet, lang, compiler):
 
     metadata['includes'] = lib
 
-    return metadata
+    return iet, metadata