Solve linear system repeatedly without allocation

Maybe this helps: Non-allocating matrix inversion - #7 by stevengj