math.functions: adding lgamma

diff --git a/basis/math/functions/functions.factor b/basis/math/functions/functions.factor
index 0c1ce94b031b12822dbe012aaa5435dc262fcdf9..5ac743b76c98ed27eba6a3359274b2cf7461b48a 100644 (file)
--- a/basis/math/functions/functions.factor
+++ b/basis/math/functions/functions.factor
@@ -198,6 +198,12 @@ M: complex log >polar [ flog ] dip rect> ; inline
 
 : logn ( x n -- y ) [ log ] bi@ / ;
 
+GENERIC: lgamma ( x -- y )
+
+M: float lgamma flgamma ;
+
+M: real lgamma >float lgamma ;
+
 <PRIVATE
 
 : most-negative-finite-float ( -- x )

diff --git a/basis/math/libm/libm.factor b/basis/math/libm/libm.factor
index 1149695414c65cd78894038a24e9178ccfc5d680..04ee64543f089c7591eb4705d05095fbd2c091d9 100644 (file)
--- a/basis/math/libm/libm.factor
+++ b/basis/math/libm/libm.factor
@@ -73,3 +73,6 @@ FUNCTION-ALIAS: fasinh
 
 FUNCTION-ALIAS: fatanh
     double atanh ( double x )
+
+FUNCTION-ALIAS: flgamma
+    double lgamma ( double x )