|
Georgi Guninski on Fri, 16 Feb 2024 16:56:17 +0100
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
|
Re: Numerical instability in extracting base B digits of the fractional part of real
|
- To: pari-dev@pari.math.u-bordeaux.fr
- Subject: Re: Numerical instability in extracting base B digits of the fractional part of real
- From: Georgi Guninski <gguninski@gmail.com>
- Date: Fri, 16 Feb 2024 17:55:54 +0200
- Delivery-date: Fri, 16 Feb 2024 16:56:17 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20230601; t=1708098975; x=1708703775; darn=pari.math.u-bordeaux.fr; h=content-transfer-encoding:to:subject:message-id:date:from :in-reply-to:references:mime-version:from:to:cc:subject:date :message-id:reply-to; bh=Si5lLnMb71kd8f5vkfC4dBqKBuHhpSqP6QhrQJNfno8=; b=SdEtrTCl708KH4Mpj6nx1MCkE64GUy1FfT97DD0QVzLiu4sOroz2DbzklNjvbB2Z8v b9Yk0GnA2H2DVjkHDAFnOd/HJPRmTUd+W5nbVn46NPaMwCg73tP2Fac6qhREX2qDOI5m AgMTCIEtaGVwZDnHTlnG6s2m1zBE/NtLVgUP/4ceSeBFjKV+F9nd2snBtWRCSboqI/vM btah+ZYTWOSvDzN+oUu1FfWV5wATsBRIL3/wDRdUEq7kw6EbxCXFykcFBYwVC4fS2u// X/pHzXdcfsCwb6ljYK/JFg9bafJ8jYbLpiy1675s0aLgblbSB8S8y/p4bYuqS4Ug0DKN 3i7g==
- In-reply-to: <Zc9EKwc0OJ8gF/b6@seventeen>
- References: <CAGUWgD9jh04yWKo8gnd=oT8sxt---hQiq65tGc0cM6Mxnk3wdw@mail.gmail.com> <Zc9EKwc0OJ8gF/b6@seventeen>
On Fri, Feb 16, 2024 at 1:20 PM Bill Allombert
<Bill.Allombert@math.u-bordeaux.fr> 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.