Abstract.In this paper, we derive a representation for the Drazin inverse of a block matrix M=( A B C D ) under the assumptions 𝐴𝐴𝜋𝐵 = 0, 𝐶𝐴𝜋𝐵 ...
Seems like GE on a 2x2 matrix should be pretty darn fast.... especially if you use scalar * vector type operations.