From a5ae37b711c77bbb013fe726070dc1d17e4c9c66 Mon Sep 17 00:00:00 2001
From: Dominic Kempf <dominic.kempf@iwr.uni-heidelberg.de>
Date: Fri, 16 Dec 2016 15:31:52 +0100
Subject: [PATCH] Change size to actual size
---
python/dune/perftool/sumfact/sumfact.py | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/python/dune/perftool/sumfact/sumfact.py b/python/dune/perftool/sumfact/sumfact.py
index 86731a14..d363fc0c 100644
--- a/python/dune/perftool/sumfact/sumfact.py
+++ b/python/dune/perftool/sumfact/sumfact.py
@@ -247,7 +247,8 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
if get_option('fastdg'):
ft = get_global_context_value("form_type")
if ft=='residual':
- globalarg(accum, dtype=np.float64, shape=(ansatz_lfs.index,), managed=False)
+ size = basis_functions_per_direction() ** world_dimension()
+ globalarg(accum, dtype=np.float64, shape=(size,), managed=False)
assignee = prim.Subscript(prim.Variable(accum), (test_lfs.index,))
expression = prim.Sum((assignee,result))
instruction(assignee=assignee,
@@ -259,7 +260,8 @@ def generate_accumulation_instruction(visitor, accterm, measure, subdomain_id):
else:
assert ft=='jacobian'
# palpo TODO: think about it
- globalarg(accum, dtype=np.float64, shape=(ansatz_lfs.index, test_lfs.index), managed=False)
+ size = basis_functions_per_direction() ** world_dimension()
+ globalarg(accum, dtype=np.float64, shape=(size, size), managed=False)
assignee = prim.Subscript(prim.Variable(accum), (ansatz_lfs.index, test_lfs.index))
expression = prim.Sum((assignee,result))
instruction(assignee=assignee,
--
GitLab