[PATCH 20/25] math: Use tgamma from CORE-MATH

[Paul Zimmermann]

       Hi DJ,

> > +static inline double
> > +splt (double x, double *xl)
> > +{
> > +  static const double off = 0x1.8p27;
> > +  double xh = (x + off) - off;
> > +  *xl = x - xh;
> > +  return xh;
> > +}
> 
> These could really use comments ;-)

you are right. One could add:

/* split x into a high part xh (return value) and a low part xl,
   where ulp(xh) = 2^-25. If |x| <= 2^26, the splitting is exact:
   x = xh + xl */

Paul