|
Bill Allombert on Fri, 03 Feb 2023 20:17:42 +0100
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
|
Re: Fortran Pari precision
|
- To: pari-users@pari.math.u-bordeaux.fr
- Subject: Re: Fortran Pari precision
- From: Bill Allombert <Bill.Allombert@math.u-bordeaux.fr>
- Date: Fri, 3 Feb 2023 20:16:13 +0100
- Arc-authentication-results: i=1; smail; arc=none
- Arc-message-signature: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1675451769; c=relaxed/relaxed; bh=7sWLb0ANtFgR9xxBTiPO5+N+pcsntCwSfNq2SqPSK10=; h=DKIM-Signature:Date:From:To:Subject:Message-ID:Mail-Followup-To: References:MIME-Version:Content-Type:Content-Disposition: In-Reply-To; b=NThDFRq89FWiBWgs3Hi5n1CDmTfPet/AdQ9bgKX/b97g8Gk/CfnWkZozjp78sYzGalyrbZXwA49rHnkoN5bUr9tP/AyL3klVVqrPFGft5r7F44IlxiUl4VGiZl6jxJgGntQkgn+qhPHV6tfr9OI0ni/5HH2iA7fi/Jo7KoC9dBWCgX+yK5lTaXjgQTsBaXUAlyQD/+wohNuEWCKM6MCBNETTVDP+kYtj9QwEkCY1hZQcwnxIrDhC049G3UjbKo4fuBxy3/Fl2PNBM8UNg1ScftTvYtWROmGjiO4KJ/cGgKGSt5ksJDnPGqA+ljsPXXRAlsHAo9Z27iNgAEFvHQiq6EBHb/w+cw/yi86G+x56dN6XCMVdmTiAEkIN/yF4EEzANeIUs5GT/N3JEsHNJu6BnGS96fxF7kmVsioFEb3eAPSnqpsjWJawUuWMy9XRRR+/ABGNtTbIbtQA7vRYKYV4cWKgfMtzy06nfo47SqxQ/UNu7e00WAa39ES/h/1o6rKz4C9hz1/RiejTuKJczCcp2b9668pF8p3F8bJ8DWzAMCexbrJIdKbjfGu0G1I+zB7ixIOfFKOTvbCfGhqcEWnC88xQrPhcz1eiQmshJ4ShR7TYgTfUXCUfHRkw0i73nMpTogDzS8l193pUlKB+4aX+tHpBVmblWHJ7uZg7VdCS7oc=
- Arc-seal: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1675451769; cv=none; b=CORJbQiWEz/1aNzF38ygrX9YeD89XtKbvAWLv43A8YN7AWN6HTvKZnVdh9LhBxr6v1cIiZ+S1cONaXKZoHhZ/NOyPFgW0RGCMY1qKOCjdK2Nukb9C61dbZUWn1HQjuyKtmWENFGx7aH6BHi/24si7zuhPoeaTpRkQ860d136FTSuPF4qKJJ+7vl00VVnly0PB0497q1h8sSQI0gPj7kSt6JKURJxOP3obEVFVcr2T9WUX9Uo8UY9CAZ0Vsb7wnmXV2TOVL19zuU7M2lkesMBpvBqm3HzoHsc1KdjtZvV8uPlouNW/fvGDmcmQqPU344FEFKJlfwRTC31BHd1FFSFLxQse3k8OTbxw1HNeh37AG2SUFEcTTagV6R2AHOoIFqlWzhPu1yT2AFXgY96lmu/EsCcM2U9v1AKHJeIvz/XX+11f7CgpscF9gG82Z7TLhiNzN+kojSJDqF5M0x0cci3mk68F3khWiynu3zQM4cJ9yoTTFacqKP7xf2v9qCCn8Y2aqxIwIjmjIZqIAKTWBwZNnewYtTzqOqgZz9cK0UtNpmBcNTcw4711SHgnSuMTm/opp49vSCd4MYOg6ubbpKgWpzO2TI5sqH28i2wEWMNgzQStNbqV8CzxyaXW9KiIHlRAksa7zVEaaiM6ucswi2ZJJ5a7ZZn6CkQZcSe39xSvTw=
- Authentication-results: smail; arc=none
- Delivery-date: Fri, 03 Feb 2023 20:17:42 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=math.u-bordeaux.fr; s=2022; t=1675451769; bh=7sWLb0ANtFgR9xxBTiPO5+N+pcsntCwSfNq2SqPSK10=; h=Date:From:To:Subject:References:In-Reply-To:From; b=OP+YGTtxkqm3lGCG64q8YIpgZsL8F/FLJSt6QYiv62z1fwHD079rIGsgEOIv1u91F yzd6mkl5haRu0rWeRIEuzm9RVBnvYDlGBtGyeFxWa27IIwQrFNXQhSAlRKgUgSSKLy fME9F8Y5NFW0UGvq8QSGWNOQR+9Tsk1peGc5MJJXwxhA4deA8ivkwtlX9yQRRMViIO mZ2b4wZBMa32Hc11B3Ic2SzIMbx4qIyRvXLCMfycU6UOp3MXqF3jmfHAY4VA5/YAXt 51mQjl+DaVuk3P6emqR/A0sUkFYrDt0FDvkpC91Y7xN1vgyXbRB+Qg5oHUmVwZ4zAA N3gpiKPHkvvzaF6ZDlpQ6jT0hAbduthCrTU3bDo8OjS1QDsgNJluAac1pCCv0Ce+q/ TYnTRfR4Hw4axupMLm+qaUDh6bFrZmrcGvdW8N4TWTdL9kXDQn20Il0K8uf26mg8lR 8hz0DkI/MWoNkfEQ3JAXVVstCOxOUbj9zBZeFFoCcH/3dLXZ4wiA0C5TNAcy26I9Iq g9eSJDHuhzIbnWisqNVw7M0iWenvxS54QaMty395eglAY6s4nLuewyWDYoIHD+WHRa 6JwccRV6xoMiSe70cucrRyq/fBHnGOKpMRwcFvTnQCd3lAL0aAOEYYfdF1j3kuoInS cniNfS0FyiZKL+L065CoqAoA=
- In-reply-to: <LO2P302MB0282EE5D0951D0C8E4D890F2C6D79@LO2P302MB0282.GBRP302.PROD.OUTLOOK.COM>
- Mail-followup-to: pari-users@pari.math.u-bordeaux.fr
- References: <LO2P302MB0282EE5D0951D0C8E4D890F2C6D79@LO2P302MB0282.GBRP302.PROD.OUTLOOK.COM>
On Fri, Feb 03, 2023 at 11:48:05AM +0000, Brereton, Ashley wrote:
> Hi there,
>
> I was hoping to get some clarification for the precision argument used for pari variables. For the following code snippet:
>
> ***************
>
> use ISO_C_BINDING, only : C_PTR
> use PARI
>
> type(C_PTR) :: u
> integer(kind=C_LONG) :: prec = 5
>
> CALL pari_init(10000000_8, 2_8)
> u = Pi2n(1_8, prec)
>
> *****************
>
> As far as I can tell, 'prec = 5' gives me precision to higher than quad
> precision. Now I'm limited to prec >=3 , and I'd like to use lower precision
> in some cases.
>
> Is there a workaround to do lower precision calculations? The reason is there
> are some handy functions that I'd like to use that don't require high
> precision calculations, e.g. complementary error function with complex input.
Alas, the minimal precision is 3, which correspond to 64bit.
Cheers,
Bill