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/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
