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.
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
