Skip to content
Snippets Groups Projects
Commit 5e7e4f60 authored by Dominic Kempf's avatar Dominic Kempf
Browse files

[bugfix] fix stokes_stress_numdiff

parent f4eb2aa2
No related branches found
No related tags found
No related merge requests found
......@@ -146,7 +146,7 @@ def evaluate_coefficient(element, name, container, restriction, component):
basis = name_basis(leaf_element, restriction)
if isinstance(sub_element, (VectorElement, TensorElement)):
lfs = lfs_child(lfs, idims, shape=shape_as_pymbolic(shape), symmetry=element.symmetry)
lfs = lfs_child(lfs, idims, shape=shape_as_pymbolic(shape), symmetry=element.symmetry())
from dune.perftool.pdelab.argument import pymbolic_coefficient
coeff = pymbolic_coefficient(container, lfs, index)
......@@ -184,7 +184,7 @@ def evaluate_coefficient_gradient(element, name, container, restriction, compone
basis = name_reference_gradient(leaf_element, restriction)
if isinstance(sub_element, (VectorElement, TensorElement)):
lfs = lfs_child(lfs, idims[:-1], shape=shape_as_pymbolic(shape[:-1]), symmetry=element.symmetry)
lfs = lfs_child(lfs, idims[:-1], shape=shape_as_pymbolic(shape[:-1]), symmetry=element.symmetry())
from dune.perftool.pdelab.argument import pymbolic_coefficient
coeff = pymbolic_coefficient(container, lfs, index)
......
......@@ -297,7 +297,7 @@ def determine_accumulation_space(expr, number, measure):
if len(subel.value_shape()) != 0:
from dune.perftool.pdelab.geometry import dimension_iname
idims = tuple(dimension_iname(context='arg', count=i) for i in range(len(subel.value_shape())))
lfs = lfs_child(lfs, idims, shape=subel.value_shape(), symmetry=subel.symmetry)
lfs = lfs_child(lfs, idims, shape=subel.value_shape(), symmetry=subel.symmetry())
subel = subel.sub_elements()[0]
lfsi = Variable(lfs_iname(subel, ma.restriction, count=number))
......
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