Re: Numerical instability in extracting base B digits of the fractional part of real

Georgi Guninski <[email protected]> Fri, 16 Feb 2024 17:55:54 +0200

Newsgroups	gmane.comp.mathematics.pari.devel
Message-ID	<CAGUWgD9R94BGzucES2RMAkakRr8z0yygJrAm-ontN-uWe0Pn+g@mail.gmail.com>

On Fri, Feb 16, 2024 at 1:20 PM Bill Allombert
<[email protected]> wrote:
>
> ? a=.33
> %6 = 0.33000000000000000000000000000000000000
> ? a*10-3
> %7 = 0.30000000000000000000000000000000000000
> ? b=a*10-3
> %8 = 0.30000000000000000000000000000000000000
> ? b*10-3
> %9 = -2.350988701644575016E-38
>

Thanks. Why I don't see -2.350988701644575016E-38 in my session?
It must be there, since you unrolled the first two iterations of
the loop.