Re: Any good package for mathematical function in Ada?

["Dmitry A. Kazakov" <mailbox@dmitry-kazakov.de>]

On 01/06/2020 12:48, Nasser M. Abbasi wrote:
> On 6/1/2020 5:19 AM, Dmitry A. Kazakov wrote:
> 
>> The test produces:
>>
>> J0(1)= 7.65197686557967E-01
> 
> That is good. In Mathematica
> 
> N[BesselJ[0, 1], 100]
> 
> 0.76519768655796655144971752610266322090927428975532524186154754911927\
> 89122152724401671806000989156340
> 
> In Maple
> 
> Digits:=100;
>                           Digits := 100
> evalf(BesselJ(0,1))
> 0.76519768655796655144971752610266322090927428975532524186154754911927\
> 89122152724401671806000989156340
> 
> Want 1,000 digits? 2,000 digits? all can be done.

In how many microseconds?

> I think these systems both link to GMP "GNU Multiple Precision 
> Arithmetic Library"
> for this. "There are no practical limits to the precision "
> https://en.wikipedia.org/wiki/GNU_Multiple_Precision_Arithmetic_Library

There is no best numerical library as requirements greatly vary. One 
would probably never want to use multiple precision arithmetic in 
production code. Less improbable, but still, in laboratory computations 
when other factors control precision.

As for GMP specifically, I think that arbitrary precision numeric types 
must be an integral part of Ada. Unfortunately, this would introduce 
same mess Unbounded_String did. So, for now, I would not push for them 
until the language type system mature to accommodate them smoothly.

-- 
Regards,
Dmitry A. Kazakov
http://www.dmitry-kazakov.de