Skip to content

Implement matrix free block preconditioner operators

Dominic Kempf requested to merge feature/block-preconditioner into master

This MR allows to generate local operators necessary for the block jacobi/GS preconditioners from exadune. The jacobian form is modified to only assemble into diagonal blocks or offdiagonal blocks before the action operator is applied to get the jacobian_apply form. I have verified by looking through code, but I do not have a test yet (that requires an actual solver setup, which might be a bit of work.)

This fixes #114 (closed)

Merge request reports

Loading