[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] AX=B octave vs gsl
|
From: |
Awhan Patnaik |
|
Subject: |
Re: [Help-gsl] AX=B octave vs gsl |
|
Date: |
Wed, 15 Sep 2010 08:02:59 +0530 |
On Sat, Sep 11, 2010 at 2:14 PM, Brian Gough <address@hidden> wrote:
> At Wed, 8 Sep 2010 08:01:23 +0530,
> Awhan Patnaik wrote:
>> What can I do to get as much accuracy as the Octave method?
>
> Did you compute the norm from the printed output data or the data in
> memory? Note that it is only printed to 7 digits, so you would lose
> ~11 digits of precision in that case.
>
That is precisely where I was going wrong. The loss of significant
digits while transferring data to Octave through a data file resulted
in the loss of accuracy. When the norm is computed in the C++ code all
is right. Many thanks to Francesco Abbate for helping me in this.