This stuff is mostly untested as of yet, but should not cause any crashes.
Results of lgamma are miserable caused by rounding.

Still unimplemented, though in <math.h>:

copysign, drem, infnan, isinf, isnan, rint, cbrt, erf, fabs, fmod,
frexp, j0, j1, jn, ldexp, modf, pow.

Also we will need something like test-math.c (name says it all)
to test availability and accuracy of these functions.

Go for it, test it. Eddie.
<ecd@dressler.de>