Charles Greathouse on Fri, 03 Feb 2023 20:57:11 +0100

 Re: Fortran Pari precision

• To: pari-users@pari.math.u-bordeaux.fr
• Subject: Re: Fortran Pari precision
• From: Charles Greathouse <crgreathouse@gmail.com>
• Date: Fri, 3 Feb 2023 14:55:37 -0500
• Delivery-date: Fri, 03 Feb 2023 20:57:11 +0100
• References: <LO2P302MB0282EE5D0951D0C8E4D890F2C6D79@LO2P302MB0282.GBRP302.PROD.OUTLOOK.COM> <Y91dfeaWEB7PZ7Ap@seventeen>

If 64-bit is too much, I guess that means you're looking for a 32-bit complex incomplete gamma implementation?

On Fri, Feb 3, 2023 at 2:16 PM Bill Allombert <Bill.Allombert@math.u-bordeaux.fr> wrote:
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