[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Output of the gsl_linalg_complex_LU_solve
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Output of the gsl_linalg_complex_LU_solve |
Date: |
Fri, 23 Mar 2007 15:19:49 +0000 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/21.3 Mule/5.0 (SAKAKI) |
At Thu, 22 Mar 2007 15:56:16 -0400,
Ho, Don T. wrote:
>
> The result of solving for A*x = B using the gsl_linalg_complex_LU_sovle
> is wrong.
>
> But if I do the inverse of A using gsl_linalg_complex_LU_invert() and
> then take B*(A)^-1 then the output x is correct.
Thanks for your email, I'm afraid we can't investigate any bug reports
without an example program that shows the problem.
--
best regards,
Brian Gough
(GSL Maintainer)
Network Theory Ltd,
Commercial support for GSL --- http://www.network-theory.com/gsl/