Skip to content

Keep any (near) null space attached to submatrices in PCFIELDSPLIT even if the...

Barry Smith requested to merge barry/2024-02-13/fix-nullspace-loss-pcfieldsplit into main

the nonzero structure of the original matrix changes.

    

    Also applied to PCBJACOBI and PCASM

    

    Reported-by: "Jeremy Theler (External)" <jeremy.theler-ext@ansys.com>

    

    The long-term way to handle this is to have MatGetSubMatr*() reuse the same submatrix objects even if the nonzero structure changes. This would just require a bit of recoding in multiple places.

Merge request reports