Fixup normal directions

python/dune/perftool/sumfact/geometry.py

@@ -167,8 +167,8 @@ def pymbolic_unit_outer_normal(visitor_indices):
     assert isinstance(index, int)
     if get_option("diagonal_transformation_matrix"):
         from dune.perftool.sumfact.switch import get_facedir, get_facemod
-        if index == get_facedir(Restriction.POSITIVE):
-            if get_facemod(Restriction.POSITIVE):
+        if index == get_facedir(Restriction.NEGATIVE):
+            if get_facemod(Restriction.NEGATIVE):
                 return 1, None
             else:
                 return -1, None
@@ -179,7 +179,7 @@ def pymbolic_unit_outer_normal(visitor_indices):
         return _norm(), visitor_indices
 
 
-def pymbolic_unit_outer_normal(visitor_indices):
+def pymbolic_unit_inner_normal(visitor_indices):
     index, = visitor_indices
     assert isinstance(index, int)
     if get_option("diagonal_transformation_matrix"):