From ce1c040d3fef2d9bd4fbc7e0b3ea33f08dc3bdd7 Mon Sep 17 00:00:00 2001
From: Marcel Koch <marcel.koch@uni-muenster.de>
Date: Tue, 23 Oct 2018 09:39:55 +0200
Subject: [PATCH] use Dune fv type if no vectorization is used
---
python/dune/perftool/blockstructured/geometry.py | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/python/dune/perftool/blockstructured/geometry.py b/python/dune/perftool/blockstructured/geometry.py
index 29fc849a..3211d051 100644
--- a/python/dune/perftool/blockstructured/geometry.py
+++ b/python/dune/perftool/blockstructured/geometry.py
@@ -282,7 +282,10 @@ def pymbolic_facet_jacobian_determinant():
# translate a point in the micro element into macro coordinates
def define_point_in_macro(name, point_in_micro):
dim = local_dimension()
- temporary_variable(name, shape=(dim,), managed=True)
+ if get_form_option('vectorization_blockstructured'):
+ temporary_variable(name, shape=(dim,), managed=True)
+ else:
+ temporary_variable(name, shape=(dim,), shape_impl=('fv',))
# point_macro = (point_micro + index_micro) / number_of_blocks
# where index_micro = tensor index of the micro element
--
GitLab