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