hypot.c

Go to the documentation of this file.

/* public domain rewrite of hypot */

#include "ruby/missing.h"
#include <math.h>

double hypot(double x, double y)
{
    if (x < 0) x = -x;
    if (y < 0) y = -y;
    if (x < y) {
        double tmp = x;
        x = y; y = tmp;
    }
    if (y == 0.0) return x;
    y /= x;
    return x * sqrt(1.0+y*y);
}