Add tgammaf() as a simple wrapper around tgamma().

5 files changed, 63 insertions, 9 deletions

diff --git a/lib/msun/Makefile b/lib/msun/Makefile
index 64d5898..c870352 100644
--- a/lib/msun/Makefile
+++ b/lib/msun/Makefile
@@ -63,7 +63,7 @@ COMMON_SRCS= b_exp.c b_log.c b_tgamma.c \
 	s_rint.c s_rintf.c s_round.c s_roundf.c s_roundl.c \
 	s_scalbln.c s_scalbn.c s_scalbnf.c s_signbit.c \
 	s_signgam.c s_significand.c s_significandf.c s_sin.c s_sinf.c \
-	s_tan.c s_tanf.c s_tanh.c s_tanhf.c s_trunc.c s_truncf.c \
+	s_tan.c s_tanf.c s_tanh.c s_tanhf.c s_tgammaf.c s_trunc.c s_truncf.c \
 	w_cabs.c w_cabsf.c w_drem.c w_dremf.c
 
 # Location of fpmath.h and _fpmath.h
@@ -151,7 +151,8 @@ MLINKS+=ilogb.3 ilogbf.3 ilogb.3 ilogbl.3 \
 	ilogb.3 logb.3 ilogb.3 logbf.3 ilogb.3 logbl.3
 MLINKS+=j0.3 j1.3 j0.3 jn.3 j0.3 y0.3 j0.3 y1.3 j0.3 y1f.3 j0.3 yn.3
 MLINKS+=j0.3 j0f.3 j0.3 j1f.3 j0.3 jnf.3 j0.3 y0f.3 j0.3 ynf.3
-MLINKS+=lgamma.3 gamma.3 lgamma.3 gammaf.3 lgamma.3 lgammaf.3 lgamma.3 tgamma.3
+MLINKS+=lgamma.3 gamma.3 lgamma.3 gammaf.3 lgamma.3 lgammaf.3 \
+	lgamma.3 tgamma.3 lgamma.3 tgammaf.3
 MLINKS+=log.3 log10.3 log.3 log10f.3 log.3 log1p.3 log.3 log1pf.3 log.3 logf.3
 MLINKS+=lrint.3 llrint.3 lrint.3 llrintf.3 lrint.3 llrintl.3 \
 	lrint.3 lrintf.3 lrint.3 lrintl.3
diff --git a/lib/msun/Symbol.map b/lib/msun/Symbol.map
index 2065b22..1c47130 100644
--- a/lib/msun/Symbol.map
+++ b/lib/msun/Symbol.map
@@ -201,4 +201,5 @@ FBSD_1.1 {
 	sinl;
 	cosl;
 	tanl;
+	tgammaf;
 };
diff --git a/lib/msun/man/lgamma.3 b/lib/msun/man/lgamma.3
index aaa4d7c..ea6eb6f 100644
--- a/lib/msun/man/lgamma.3
+++ b/lib/msun/man/lgamma.3
@@ -40,7 +40,8 @@
 .Nm gamma_r ,
 .Nm gammaf ,
 .Nm gammaf_r ,
-.Nm tgamma
+.Nm tgamma ,
+.Nm tgammaf
 .Nd log gamma functions, gamma function
 .Sh LIBRARY
 .Lb libm
@@ -67,6 +68,8 @@
 .Fn gammaf_r "float x" "int *signgamp"
 .Ft double
 .Fn tgamma "double x"
+.Ft float
+.Fn tgammaf "float x"
 .Sh DESCRIPTION
 .Fn lgamma x
 and
@@ -80,7 +83,6 @@ return ln\||\(*G(x)| where
 .\}
 .if n \
 return ln\||\(*G(x)|.
-.Pp
 The external integer
 .Fa signgam
 returns the sign of \(*G(x).
@@ -94,6 +96,13 @@ and
 .Fn lgammaf x
 but the caller must provide an integer to store the sign of \(*G(x).
 .Pp
+The
+.Fn tgamma x
+and
+.Fn tgammaf x
+functions return \(*G(x), with no effect on
+.Fa signgam .
+.Pp
 .Fn gamma ,
 .Fn gammaf ,
 .Fn gamma_r ,
@@ -106,9 +115,6 @@ are deprecated aliases for
 and
 .Fn lgammaf_r ,
 respectively.
-.Fn tgamma x
-returns \(*G(x), with no effect on
-.Fa signgam .
 .Sh IDIOSYNCRASIES
 Do not use the expression
 .Dq Li signgam\(**exp(lgamma(x))
@@ -151,9 +157,11 @@ will underflow.
 .Xr math 3
 .Sh STANDARDS
 The
-.Fn lgamma
+.Fn lgamma ,
+.Fn lgammaf ,
+.Fn tgamma ,
 and
-.Fn tgamma
+.Fn tgammaf
 functions are expected to conform to
 .St -isoC-99 .
 .Sh HISTORY
diff --git a/lib/msun/src/math.h b/lib/msun/src/math.h
index 94bd919..85b8483 100644
--- a/lib/msun/src/math.h
+++ b/lib/msun/src/math.h
@@ -331,6 +331,7 @@ float	erff(float);
 float	erfcf(float);
 float	hypotf(float, float);
 float	lgammaf(float);
+float	tgammaf(float);
 
 float	acoshf(float);
 float	asinhf(float);
diff --git a/lib/msun/src/s_tgammaf.c b/lib/msun/src/s_tgammaf.c
new file mode 100644
index 0000000..9993d91
--- /dev/null
+++ b/lib/msun/src/s_tgammaf.c
@@ -0,0 +1,43 @@
+/*-
+ * Copyright (c) 2008 David Schultz <das@FreeBSD.ORG>
+ * All rights reserved.
+ *
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ *    notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ *    notice, this list of conditions and the following disclaimer in the
+ *    documentation and/or other materials provided with the distribution.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
+ * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
+ * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+ * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+ * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+ * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+ * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+ * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+ * SUCH DAMAGE.
+ */
+
+#include <sys/cdefs.h>
+__FBSDID("$FreeBSD$");
+
+#include <math.h>
+
+/*
+ * We simply call tgamma() rather than bloating the math library with
+ * a float-optimized version of it. The reason is that tgammaf() is
+ * essentially useless, since the function is superexponential and
+ * floats have very limited range.
+ */
+float
+tgammaf(float x)
+{
+
+	return (tgamma(x));
+}