From: Brian Gough <bjg@network-theory.co.uk>
To: Frank Reininghaus <frank78ac@googlemail.com>
Cc: gsl-discuss@sourceware.org
Subject: Re: Examining bug #21838
Date: Mon, 18 Feb 2008 09:56:00 -0000 [thread overview]
Message-ID: <87wsp2pozy.wl%bjg@network-theory.co.uk> (raw)
In-Reply-To: <47B4C827.1040004@googlemail.com>
At Fri, 15 Feb 2008 00:00:55 +0100,
Frank Reininghaus wrote:
> Maybe these (at least for me) surprising observations could be used to
> isolate the cause of the bug. Does anyone have a good idea?
It seems that in double precision the error term (delta_frac - 1) is
exactly zero at some point, so the continued fraction terminates. It
looks like this happens before the desired accuracy is actually
reached though, so it seems to be undesirable. In extended precision
that doesn't happen, the loop continues up to its maximum iteration
limit and NAN is returned to signal that.
If you add #include <gsl/gsl_ieee_utils.h>
and gsl_ieee_env_setup();
then run with GSL_IEEE_MODE=double-precision ./a.out
you can see the difference.
prev parent reply other threads:[~2008-02-18 9:56 UTC|newest]
Thread overview: 5+ messages / expand[flat|nested] mbox.gz Atom feed top
2008-02-14 23:01 Frank Reininghaus
2008-02-15 9:07 ` Jason Stover
2008-02-20 18:04 ` Brian Gough
2008-02-16 18:22 ` Brian Gough
2008-02-18 9:56 ` Brian Gough [this message]
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=87wsp2pozy.wl%bjg@network-theory.co.uk \
--to=bjg@network-theory.co.uk \
--cc=frank78ac@googlemail.com \
--cc=gsl-discuss@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).