Dear lazyweb,
I'm working on optimizing for a non-x86 target processor and left no
choice but to start writing optimized math functions.
So far I've compared
glibc
Sun solaris open sourced code
NetBSD mathlib (Which originated from Sun circa 1993)
netlib
{REDACTED-vendor-name} modified version of netlib
Surprisingly, glibc does a pretty respectable job in terms of
accuracy, but alas it's certainly not the fastest.
------------
One of the faster versions tested sometimes incorrectly
rounds for medium size and larger arguments, i.e. the least
significant bit differs from the other versions.
---
sin(6e5) =
bfcb40318b8c1728 vs
bfcb40318b8c1727 =
-0.212896531236929975605676 (REDACTED-vendor-name) vs
-0.212896531236929947850101 (nearest) vs
-0.212896531236929960235352 (higher precision) ...
-------
Question: Is the above acceptable to your scientific codes in general?
_______________________________________________
Beowulf mailing list, Beowulf@beowulf.org sponsored by Penguin Computing
To change your subscription (digest mode or unsubscribe) visit
http://www.beowulf.org/mailman/listinfo/beowulf
