PETSc: add higher level functions to perform global linear algebra operations that performs `NumberingSubset` checks
Currently those operations are defined at the ThirdParty
module level, where NumberingSubset
aren't yet defined. The idea would be to provide a thin wrapper in Core
, which would check in debug mode if the matrices involved match.
This would enable early detection of something wrong in models.
It is closely related to #1889, for which I considered doing this, dismissed it in one commit yesterday... before changing my mind while updating the tests about matrix operations (I figured out how to implement it quickly and efficiently).