[Beowulf] Optimized math routines/transcendentals

Prentice Bisbal pbisbal at pppl.gov
Tue May 31 12:36:15 PDT 2016


On 04/29/2016 02:23 PM, C Bergström wrote:
> 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?
>
Is this the sort of problem John Gustafson's Unums was meant to eliminate?

http://www.johngustafson.net/unums.html

Unfortunately, that seems to be the best (only?) page on the Internet 
about Unums, other than his Wikipedia entry:

https://en.wikipedia.org/wiki/John_Gustafson_(scientist)#Unums

--
Prentice




More information about the Beowulf mailing list