diff options
author | Taru Karttunen <taruti@taruti.net> | 2011-03-30 15:46:40 +0300 |
---|---|---|
committer | Taru Karttunen <taruti@taruti.net> | 2011-03-30 15:46:40 +0300 |
commit | e5888a1ffdae813d7575f5fb02275c6bb07e5199 (patch) | |
tree | d8d51eac403f07814b9e936eed0c9a79195e2450 /sys/src/ape/lib/ap/math/hypot.c |
Import sources from 2011-03-30 iso image
Diffstat (limited to 'sys/src/ape/lib/ap/math/hypot.c')
-rwxr-xr-x | sys/src/ape/lib/ap/math/hypot.c | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/sys/src/ape/lib/ap/math/hypot.c b/sys/src/ape/lib/ap/math/hypot.c new file mode 100755 index 000000000..25457d8a9 --- /dev/null +++ b/sys/src/ape/lib/ap/math/hypot.c @@ -0,0 +1,29 @@ +#include <math.h> +/* + * sqrt(a^2 + b^2) + * (but carefully) + */ + +double +hypot(double a, double b) +{ + double t; + + if(a < 0) + a = -a; + if(b < 0) + b = -b; + if(a > b) { + t = a; + a = b; + b = t; + } + if(b == 0) + return 0; + a /= b; + /* + * pathological overflow possible + * in the next line. + */ + return b * sqrt(1 + a*a); +} |