source: vendor/python/2.5/Python/hypot.c@ 3397

/* hypot() replacement */

#include "pyconfig.h"
#include "pyport.h"

double hypot(double x, double y)
{
        double yx;

        x = fabs(x);
        y = fabs(y);
        if (x < y) {
                double temp = x;
                x = y;
                y = temp;
        }
        if (x == 0.)
                return 0.;
        else {
                yx = y/x;
                return x*sqrt(1.+yx*yx);
        }
}