Georgi Guninski on Thu, 22 Feb 2024 08:13:58 +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: Thu, 22 Feb 2024 09:13:31 +0200
- Delivery-date: Thu, 22 Feb 2024 08:13:58 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20230601; t=1708586036; x=1709190836; darn=pari.math.u-bordeaux.fr; h=to:subject:message-id:date:from:in-reply-to:references:mime-version :from:to:cc:subject:date:message-id:reply-to; bh=HmvFVSpkvMy5UOjdJKrERAIm2uZyZKqpIHxf/azwk+A=; b=H/KGgq2YZvIyXkNpybuuPaqoC/DQhc757e5PcduxRQfSQ6/imFYuf3pUWDU2uO9LpS jtvgF8vK84Fy03doe9qTezE2kCHJPM4CI3OiqRiZtvjdXJEUkwxBG+j1hKb9gPynlQQT cqNoOOQXgCL/6Ql+W4RglzJKAMIIvaKgDzxPPx11Of5c+tSua3WHU7PPrTI0oCvTEsNN nKOwIsZWC+We5P30QtRwQy97k/8oCCQi5kZYOGGomKqMYXRkFb439ybu6xadXQhEOMtR MXOIizJCvlc5mN2uGRf+iFpLk3dOqc0U4o/B23CsEW6MXtZE0AEckdAZYcT/ZEZk+fAB KKiA==
- In-reply-to: <CAGUWgD8tZP-8T1jStNjMCAqG59QNRfNAJaMCNQYq+q_BV42_hA@mail.gmail.com>
- References: <CAGUWgD9jh04yWKo8gnd=oT8sxt---hQiq65tGc0cM6Mxnk3wdw@mail.gmail.com> <Zc9EKwc0OJ8gF/b6@seventeen> <CAGUWgD9R94BGzucES2RMAkakRr8z0yygJrAm-ontN-uWe0Pn+g@mail.gmail.com> <Zc+InbDj44/YHxrK@seventeen> <CAGUWgD8tZP-8T1jStNjMCAqG59QNRfNAJaMCNQYq+q_BV42_hA@mail.gmail.com>
Some observations:
1. With precision 100 I get correct result.
2. The issue is triggered by two calls of C=A*B, starting with A=0.3
works correctly.