Guido van Rossum | 57b1822 | 1996-08-29 18:10:41 +0000 | [diff] [blame] | 1 | /* hypot() replacement */ |
2 | |||||
3 | #include "config.h" | ||||
4 | #include "myproto.h" | ||||
5 | #include "mymath.h" | ||||
6 | |||||
Thomas Wouters | f70ef4f | 2000-07-22 18:47:25 +0000 | [diff] [blame^] | 7 | double hypot(double x, double y) |
Guido van Rossum | 57b1822 | 1996-08-29 18:10:41 +0000 | [diff] [blame] | 8 | { |
9 | double yx; | ||||
10 | |||||
11 | x = fabs(x); | ||||
12 | y = fabs(y); | ||||
13 | if (x < y) { | ||||
14 | double temp = x; | ||||
15 | x = y; | ||||
16 | y = temp; | ||||
17 | } | ||||
18 | if (x == 0.) | ||||
19 | return 0.; | ||||
20 | else { | ||||
21 | yx = y/x; | ||||
22 | return x*sqrt(1.+yx*yx); | ||||
23 | } | ||||
24 | } |