Vectorized block diagonal mat inversion

[OmarDuran]

Proposed changes

This PR refactor and introduce a vectorized solution for the function invert_diagonal_blocks_python in file matrix_operations.py. The function is 2x faster than numba function for a large block diagonal matrix with irregular block sizes.

Types of changes

What types of changes does this PR introduce to PorePy?
Put an x in the boxes that apply.

• Minor change (e.g., dependency bumps, broken links).
• Bugfix (non-breaking change which fixes an issue).
• New feature (non-breaking change which adds functionality).
• Breaking change (fix or feature that would cause existing functionality to not work as expected).
• Testing (contribution related to testing of existing or new functionality).
• Documentation (contribution related to adding, improving, or fixing documentation).
• Maintenance (e.g., improve logic and performance, remove obsolete code).
• Other:

Checklist

Put an x in the boxes that apply or explain briefly why the box is not relevant.

• The documentation is up-to-date.
• Static typing is included in the update.
• This PR does not duplicate existing functionality.
• The update is covered by the test suite (including tests added in the PR).
• If new skipped tests have been introduced in this PR, pytest was run with the --run-skipped flag.

[keileg]

Control question: You have verified that the inverted matrices are correct?

Given the performance results, I think we should make this into the default inverter. This will require changing parameter values in mpfa, mpsa and biot. A repository-wide search for 'inverter' may be needed.