[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Possible bug in Hessenberg-Triangular Reduction.
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Possible bug in Hessenberg-Triangular Reduction. |
Date: |
Tue, 20 Sep 2011 17:15:10 +0100 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.3 Mule/6.0 (HANACHIRUSATO) |
At Fri, 2 Sep 2011 22:19:36 +0200,
Miguel García Torres wrote:
> I have tried to port the Hessenberg-Triangular Reduction from gsl to Java.
> However I found some wrong outputs. I sent an email to gsl-help and told
> me the it could be a bug in the code:
Thanks for the email. I checked that the output from your code
satisfied the decompositions A=U H V^T and B=U R V^T. A and B were
reproduced correctly in double and extended precision so I think the
results are correct. If you have a case where that is not satisified
please let me know.
Brian Gough