From b9a133e5ad6038ca7008a611966897be732f1416 Mon Sep 17 00:00:00 2001
From: Damien George <damien.p.george@gmail.com>
Date: Tue, 4 Sep 2018 17:00:35 +1000
Subject: [PATCH] lib/libm/wf_tgamma: Fix tgammaf handling of -inf, should
 return nan.

---
 lib/libm/wf_tgamma.c | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/lib/libm/wf_tgamma.c b/lib/libm/wf_tgamma.c
index 64b2488d1d..3ff05f331d 100644
--- 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