This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
transpose copy at end of cholesky_decomp
- From: James Bergstra <james dot bergstra at umontreal dot ca>
- To: gsl-discuss <gsl-discuss at sources dot redhat dot com>
- Date: Sat, 25 Mar 2006 20:34:59 -0500
- Subject: transpose copy at end of cholesky_decomp
I am wondering, why does the routine gsl_linalg_cholesky_decomp(gsl_matrix*A)
include nested for loops at the end that transpose-copy the lower triangle of A
into the upper triangle of A?
Could the library expose another routine in which the post-processing is at
least optional?
gsl_linalg_cholesky_decomp_lower(gsl_matrix*A) for example?
--
james bergstra
http://www-etud.iro.umontreal.ca/~bergstrj