🚀 go-pugleaf

In article <9002220012.AA05757@csvax.caltech.edu>
andy@CSVAX.CALTECH.EDU (Andy Fyfe) suggests the following patch
>
>*** math-68881.h	Wed Feb 14 01:07:28 1990
>--- /usr/local/lib/gcc-include/math-68881.h	Wed Feb 21 12:01:32 1990
>***************
>*** 250,253 ****
>--- 250,258 ----
>  }
>  
>+ __inline static const double hypot (const double x, const double y)
>+ {
>+     return sqrt(x*x + y*y);
>+ }
>+ 
>  __inline static const double pow (const double x, const double y)
>  {

That may be fine if you're VERY concerned about speed, but it
should be worth a small sacrifice to assure more accurate results;
the version of hypot() above can experience unwarranted loss of
precision, overflow, or destructive underflow.

__inline static const double
hypot(const double x, const double y)
{
    register double u;
    double sqrt();
#define ABS(x) ((x) < 0.0 ? -(x) : (x))

    if (ABS(x) >= ABS(y)) {
	if (x == 0.0)
	    return (0.0);
	u = y / x;
	return (x * sqrt(1.0 + u * u));
    }
    else {
	u = x / y;
	return (y * sqrt(1.0 + u * u));
    }
#undef ABS
}
-- 
Jim Meyering      meyering@cs.utexas.edu