Re: [Help-gsl] GSL linear system solver question

[Brian Gough]

At Thu, 23 Aug 2007 10:46:50 -0500,
Gen Lu wrote:
> I'm looking for a linear system solver written in C/C++.
> Since the systems I wanna solve are very simple, just  general square
> matrices, so I think LU method in GSL is good enough. But I have a
> problem about it, can LU or any other solvers in GSL distinguish the 3
> cases of solutions (unique solution, non-solution and infinite
> solutions)? If so, how can I know which cases the solution is?

The GSL QRPT decomposition and SVD both give the rank of the matrix
and an orthonormal basis for its range, which can be used to test if
there is a solution.

-- 
Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/