This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Trouble with Cholesky
Brian Gough wrote:
> The matrix needs to be decomposed into Cholesky form -- if you
> add gsl_linalg_cholesky_decomp (m); before calling the solver it
> should work.
Ahhhh! I thought that would be handled internally. I was already
wondering why m is preserved (as opposed to the way the HH solver
handles this), but I didn't draw the right conclusion. Indeed, the docu
does not tell me that I have first to triangularize m...
Thanks a million!!
Frank
--
Dr Frank Boehme | Email: f.boehme@cs.ucc.ie
National University of Ireland, Cork | phone: +353-21-4903163
Dept of Computer Science | fax: +353-21-4903113
Cork, Ireland |