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

Apply block preconditioner trafo

parent 2d8bcff3
No related branches found
No related tags found
No related merge requests found
...@@ -819,6 +819,13 @@ def generate_jacobian_kernels(form, original_form): ...@@ -819,6 +819,13 @@ def generate_jacobian_kernels(form, original_form):
from dune.perftool.ufl.preprocess import preprocess_form from dune.perftool.ufl.preprocess import preprocess_form
jacform = preprocess_form(jacform).preprocessed_form jacform = preprocess_form(jacform).preprocessed_form
if get_form_option("block_preconditioner_diagonal"):
from dune.perftool.ufl.transformations.blockpreconditioner import diagonal_block_jacobian
jacform = diagonal_block_jacobian(jacform)
if get_form_option("block_preconditioner_offdiagonal"):
from dune.perftool.ufl.transformations.blockpreconditioner import offdiagonal_block_jacobian
jacform = offdiagonal_block_jacobian(jacform)
operator_kernels = {} operator_kernels = {}
with global_context(form_type="jacobian"): with global_context(form_type="jacobian"):
if get_form_option("generate_jacobians"): if get_form_option("generate_jacobians"):
......
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