At Fri, 12 Jun 2009 18:57:40 -0400, Erwin Kalvelagen wrote: > gsl_linalg_LU_decomp should return an error code when the input matrix is > singular. See any numerical library for examples. I believe ignoring this > makes this routine of little value.
Thanks for your email. The LU decomposition is valid for a singular matrix so we don't return an error in that case. It would make sense to return an error code from gsl_linalg_LU_svx though, I will add one there. -- Brian Gough (GSL Maintainer) Support freedom by joining the FSF http://www.fsf.org/associate/support_freedom/join_fsf?referrer=37 _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
