This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Problem with matrices bigger than 2x2 using gsl_linalg_LU_invert
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Philippe Huber <huberph at infomaniak dot ch>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Mon, 19 Aug 2002 19:31:50 +0100 (BST)
- Subject: Re: Problem with matrices bigger than 2x2 using gsl_linalg_LU_invert
- References: <001101c2476d$18a7c250$7e28c281@osiris>
Philippe Huber writes:
> First of all, thank you for having created such a great library. I have
> a small problem with gsl_linalg_LU_invert. I am using a 4x4 matrix which
> is not singular. I made a LU decomposition using gsl_linalg_LU_decomp.
> Then, I send it to gsl_linalg_LU_invert and it crashes. I tried with a
> 2x2 matrix and it works perfectly well. I verified that the program also
> crashes with a 3x3 matrix. Is it a bug or did I make something wrong.
> Please apologize if I don't send my code, it is well too big.
Please send a small example which reproduces the bug, with details of
your operating system, compiler, etc.
regards
Brian Gough