This is the mail archive of the
gcc@gcc.gnu.org
mailing list for the GCC project.
Re: real.c implementation
- From: Roger Sayle <roger at eyesopen dot com>
- To: Richard Henderson <rth at redhat dot com>
- Cc: <gcc at gcc dot gnu dot org>
- Date: Fri, 18 Oct 2002 20:26:08 -0600 (MDT)
- Subject: Re: real.c implementation
Brad Lucier writes:
> Richard Henderson writes:
> > Perhaps I'm being naieve. If this is incorrect, stop me now.
> > An explanation of why this is wrong would also be appreciated.
>
> I don't know. What I do know is that having correct intermediate
> results to 2k+2 bits after doing all the conversion stuff does guarantee
> the correct target result to k bits. Something similar will come up
> with sqrt.
A reference on this topic recommended by Brad that I found interesting
(even if I didn't understand much of it) is available on the net:
"High Precision Division and Square Root", Alan H. Karp and Peter
Markstein, HP Labs Technical Report 93-93-42, June 1993.
http://www.hpl.hp.com/techreports/93/HPL-93-42.pdf
Roger
--