While testing the linear algebra functions provided in gsl I came across the following.
With appropriate setup the call to, gsl_linalg_LU_invert(squareMatrix, perm, squareInverse); gave the inverse of the transpose of squareMatrix. I verified this using octave. Note that squareMatrix is a simple 3x3 matrix of Gaussian variates. Is this the expected behaviour? Regards, Philip _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
