> The problem with picking a random open (or closed, reallyl) software
> product
> for numerical work and using it is that you need to figure out if it
> is any good.
>
This is why using minuit for the fitting library is good idea. It's been
the standard fitting library for High Energy Physics for over 20 years.
It's bugs and limitations are now very well known.
Cheers
Martin
> A quick look at numpy suggests that the numerical quality is questionable.
> Take for example
>
> static double atanh(double x)
> {
> return 0.5*log((1.0+x)/(1.0-x));
> }
>
> This clearly has 100% loss of precision for small x. (Luckily this
> function is
> only used if libc fails to supply atanh.)
>
> Morten
> _______________________________________________
> gnumeric-list mailing list
> [email protected]
> http://mail.gnome.org/mailman/listinfo/gnumeric-list
>
_______________________________________________
gnumeric-list mailing list
[email protected]
http://mail.gnome.org/mailman/listinfo/gnumeric-list
