Ed,
gfortran uses libquadmath for support of its REAL(16)
intrinsic subprograms. I checked the list of intrinsics
and tgamma is current not in the list. I don't have
time to try Fortran's C interop feature to see if an
ordinary user can access __float128. While your patch
is probably useful for
Here is a patch for tammaq for negative argument pr68686.
I know about depending on ports from upstream but this was done recently
and this (tgammaq) was left out.
This patch is basically a one-liner.
I have test cases but libquadmath doesn't have a testsuite.
One test just shows alternating
