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

Merge branch 'feature/blockstructured-fix-yasp-grid' into 'master'

use the backend selector for the jacobian inverse transposed

See merge request dominic/dune-perftool!204
parents 7c7045aa 70cdfc18
No related branches found
No related tags found
No related merge requests found
from dune.perftool.generation import (backend, from dune.perftool.generation import (get_backend,
temporary_variable, temporary_variable,
instruction) instruction)
from dune.perftool.tools import get_pymbolic_basename from dune.perftool.tools import get_pymbolic_basename
from dune.perftool.options import get_option from dune.perftool.options import (get_option,
option_switch)
from dune.perftool.pdelab.geometry import (name_jacobian_determinant, from dune.perftool.pdelab.geometry import (name_jacobian_determinant,
name_jacobian_inverse_transposed,
local_dimension, local_dimension,
apply_to_global_transformation, apply_to_global_transformation,
name_facet_jacobian_determinant) name_facet_jacobian_determinant)
...@@ -20,8 +20,9 @@ def pymbolic_jacobian_determinant(): ...@@ -20,8 +20,9 @@ def pymbolic_jacobian_determinant():
# scale Jacobian according to the order of the blockstructure # scale Jacobian according to the order of the blockstructure
def pymbolic_jacobian_inverse_transposed(i, j, restriction): def pymbolic_jacobian_inverse_transposed(i, j, restriction):
name_jit = get_backend(interface="name_jit", selector=option_switch("constant_transformation_matrix"))(restriction)
return prim.Product((get_option("number_of_blocks"), return prim.Product((get_option("number_of_blocks"),
prim.Subscript(prim.Variable(name_jacobian_inverse_transposed(restriction)), (j, i)))) prim.Subscript(prim.Variable(name_jit), (j, i))))
# scale determinant according to the order of the blockstructure # scale determinant according to the order of the blockstructure
......
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