Skip to content
Snippets Groups Projects
Commit ce1c040d authored by Marcel Koch's avatar Marcel Koch
Browse files

use Dune fv type if no vectorization is used

parent 98f42427
No related branches found
No related tags found
No related merge requests found
...@@ -282,7 +282,10 @@ def pymbolic_facet_jacobian_determinant(): ...@@ -282,7 +282,10 @@ def pymbolic_facet_jacobian_determinant():
# translate a point in the micro element into macro coordinates # translate a point in the micro element into macro coordinates
def define_point_in_macro(name, point_in_micro): def define_point_in_macro(name, point_in_micro):
dim = local_dimension() 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 # point_macro = (point_micro + index_micro) / number_of_blocks
# where index_micro = tensor index of the micro element # where index_micro = tensor index of the micro element
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment