Pseudoinverse from a Bunch-Kaufman factorization - View it on GitHub