1/* public domain rewrite of hypot */ 2 3#include "ruby/missing.h" 4#include <math.h> 5 6double hypot(double x, double y) 7{ 8 if (x < 0) x = -x; 9 if (y < 0) y = -y; 10 if (x < y) { 11 double tmp = x; 12 x = y; y = tmp; 13 } 14 if (y == 0.0) return x; 15 y /= x; 16 return x * sqrt(1.0+y*y); 17} 18