[bug-gsl] Re: Bug in gsl_sf_exp_mult_e10_e error estimation (2008-01-17)
|
From: |
Brian Gough <bjg@network-theory.co.uk> |
|
Subject: |
Re: Bug in gsl_sf_exp_mult_e10_e error estimation |
|
Date: |
Thu, 17 Jan 2008 19:53:53 +0000 |
|
To: |
Szymon Jaroszewicz <sj@..........> |
|
Cc: |
bug-gsl@gnu.org |
At Mon, 14 Jan 2008 12:16:10 +0100,
Szymon Jaroszewicz wrote:
>
> There is a bug in error estimation for gsl_sf_exp_mult_e10_e:
>
> gsl_sf_exp_e10_e(10000) gives error estimate 3.911404048e-11
>
> but
>
> gsl_sf_exp_mult_e10_e(10000, 1) gives 3.911012947e-15
>
Thanks for the bug report and example program.
I've confirmed it and logged it in the bug tracker.
It does look like the error from gsl_sf_exp_mult_e10_e
is wrong.
--
Brian Gough
