Does the Boost UBLAS library have a built-in solver for solving systems of equations? The documentation implies that all the ublas solver routines require the matrix to already be in triangular form.
But, if a matrix is not in triangular form, is there anything in ublas that can reduce the matrix and then back-substitute, to solve a system of equations?