lib/libm/wf_tgamma: Fix tgammaf handling of -inf, should return nan.

dpgeorge · dpgeorge · commit b9a133e5ad60 · 2018-09-04T17:00:35.000+10:00
diff --git a/lib/libm/wf_tgamma.c b/lib/libm/wf_tgamma.c
@@ -35,6 +35,10 @@
 {
         float y;
 	int local_signgam;
+	if (!isfinite(x)) {
+	  /* special cases: tgammaf(nan)=nan, tgammaf(inf)=inf, tgammaf(-inf)=nan */
+	  return x + INFINITY;
+	}
 	y = expf(__ieee754_lgammaf_r(x,&local_signgam));
 	if (local_signgam < 0) y = -y;
 #ifdef _IEEE_LIBM
