Hi @ryanjoneil ,
I was checking the timing of your code on different sizes of random matrices, but it faces the "vector subscript out of range" error in some cases during the solve operation. The error occurs at initialize_phase_2 function; Even the changes made by Pigzaum do not fix this issue.
Also, the solver result is incorrect for a matrix with all equal elements!!