From 33ae4ed2f213782546c76a38a90d5a1f7da8d32e Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Ren=C3=A9=20He=C3=9F?= <rene.hess@iwr.uni-heidelberg.de>
Date: Tue, 17 Jan 2017 10:01:31 +0100
Subject: [PATCH] Bugfix for vectorization and direct input

---
 python/dune/perftool/sumfact/sumfact.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/python/dune/perftool/sumfact/sumfact.py b/python/dune/perftool/sumfact/sumfact.py
index e2b1b6ed..ea22fd5d 100644
--- a/python/dune/perftool/sumfact/sumfact.py
+++ b/python/dune/perftool/sumfact/sumfact.py
@@ -550,12 +550,13 @@ def sum_factorization_kernel(a_matrices,
         if l == 0 and direct_input is not None:
             # See comment bellow
             input_inames = _permute_backward(input_inames, perm)
+            inp_shape = _permute_backward(inp_shape, perm)
 
             globalarg(direct_input, dtype=np.float64, shape=inp_shape)
             if a_matrix.vectorized:
                 input_summand = prim.Call(prim.Variable("Vec4d"),
                                           (prim.Subscript(prim.Variable(direct_input),
-                                                          input_inames)))
+                                                          input_inames),))
             else:
                 input_summand = prim.Subscript(prim.Variable(direct_input),
                                                input_inames + vec_iname)
-- 
GitLab