diff --git a/python/dune/perftool/blockstructured/geometry.py b/python/dune/perftool/blockstructured/geometry.py index dc9e48e1db63bc35ae9a1031aa9fffb33a804d12..f0d1ed6bf78796a150e7a5991217a96cf2c1d318 100644 --- a/python/dune/perftool/blockstructured/geometry.py +++ b/python/dune/perftool/blockstructured/geometry.py @@ -335,9 +335,9 @@ def apply_default_to_global_transformation(name, local): expr = a_pym * local_pym[0] * local_pym[1] + b_pym * local_pym[0] + c_pym * local_pym[1] + corner_0_pym elif dim == 3: a_pym, b_pym, c_pym, d_pym, e_pym, f_pym, g_pym = coeffs_pym - expr = a_pym * local_pym[0] * local_pym[1] * local_pym[2] + b_pym * local_pym[0] * local_pym[1] + \ - c_pym * local_pym[0] * local_pym[2] + d_pym * local_pym[1] * local_pym[2] + \ - e_pym * local_pym[0] + f_pym * local_pym[1] + g_pym * local_pym[2] + corner_0_pym + expr = (a_pym * local_pym[0] * local_pym[1] * local_pym[2] + b_pym * local_pym[0] * local_pym[1] + + c_pym * local_pym[0] * local_pym[2] + d_pym * local_pym[1] * local_pym[2] + + e_pym * local_pym[0] + f_pym * local_pym[1] + g_pym * local_pym[2] + corner_0_pym) else: raise NotImplementedError @@ -362,8 +362,8 @@ def apply_constant_to_global_transformation(name, local): # global[d] = lower_left[d] + local[d] * (upper_right[d] - lower_left[d]) expr = (prim.Subscript(prim.Variable(corners), (0, dim_pym)) + - (prim.Subscript(local, (dim_pym,)) * (prim.Subscript(prim.Variable(corners), (2**dim - 1, dim_pym)) - - prim.Subscript(prim.Variable(corners), (0, dim_pym))))) + prim.Subscript(local, (dim_pym,)) * (prim.Subscript(prim.Variable(corners), (2**dim - 1, dim_pym)) - + prim.Subscript(prim.Variable(corners), (0, dim_pym)))) assignee = prim.Subscript(prim.Variable(name), (dim_pym,)) instruction(assignee=assignee, expression=expr,