Re: [Bug-gsl] Possible bug in cholesky.c: gsl_linalg_cholesky_decomp()

[Brian Gough]

Stephen Milborrow writes:
 > gsl_linalg_cholesky_decomp() appears to have a bug which causes it take
 > square roots of negative numbers when passed a non positive-definite matrix.
 > I believe this is contrary to the intent of the library which is that such
 > errors should be handled by GSL_ERROR.

Thanks for the bug report.  As you say, the error checking should go
before the sqrt().

-- 
best regards,

Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/