| author | nsz <nsz@port70.net> | 2012年03月27日 22:17:36 +0200 |
|---|---|---|
| committer | nsz <nsz@port70.net> | 2012年03月27日 22:17:36 +0200 |
| commit | bbfbc7edaf992abe1d3d09868be07c4c1cc44db7 (patch) | |
| tree | dd101f1fad64e09dd18cc7ceb71bcb0804f9aae7 /src/math/tgamma.c | |
| parent | 1b229a2098a35795aa20bd09b8c81c5143e64277 (diff) | |
| download | musl-bbfbc7edaf992abe1d3d09868be07c4c1cc44db7.tar.gz | |
| -rw-r--r-- | src/math/tgamma.c | 16 |
diff --git a/src/math/tgamma.c b/src/math/tgamma.c new file mode 100644 index 00000000..f3bbe370 --- /dev/null +++ b/src/math/tgamma.c @@ -0,0 +1,16 @@ +#include <math.h> + +// FIXME: use lanczos approximation + +double __lgamma_r(double, int *); + +double tgamma(double x) +{ + int sign; + double y; + + y = exp(__lgamma_r(x, &sign)); + if (sign < 0) + y = -y; + return y; +} |