Mix logical and physical coordinates in expression

psydac/api/discretization.py

@@ -396,16 +396,21 @@ def discretize(a, *args, **kwargs):
         kwargs['symbolic_mapping'] = mapping
 
     if isinstance(a, sym_BasicForm):
-        if isinstance(a, sym_Norm):
-            kernel_expr = TerminalExpr(a, domain)
-            if not mapping is None:
-                kernel_expr = tuple(LogicalExpr(i, domain) for i in kernel_expr)
-        else:
-            if not mapping is None:
-                a       = LogicalExpr (a, domain)
-                domain  = domain.logical_domain
-
-            kernel_expr = TerminalExpr(a, domain)
+        # if isinstance(a, sym_Norm):
+        #     kernel_expr = TerminalExpr(a, domain)
+        #     if not mapping is None:
+        #         kernel_expr = tuple(LogicalExpr(i, domain) for i in kernel_expr)
+        # else:
+        #     if not mapping is None:
+        #         a       = LogicalExpr (a, domain)
+        #         domain  = domain.logical_domain
+
+        #     kernel_expr = TerminalExpr(a, domain)
+        if not mapping is None:
+            a       = LogicalExpr (a, domain)
+            domain  = domain.logical_domain
+
+        kernel_expr = TerminalExpr(a, domain)
 
         if len(kernel_expr) > 1:
             return DiscreteSumForm(a, kernel_expr, *args, **kwargs)

pyproject.toml

@@ -32,7 +32,7 @@ dependencies = [
     'pyevtk',
 
     # Our packages from PyPi
-    'sympde==0.16.1',
+    'sympde @ git+https://github.com/pyccel/sympde@logical_n_physical-coordinates',
     'pyccel>=1.7.3',
     'gelato==0.11',