예제 #1
0
파일: choose.c 프로젝트: 6e441f9c/julia
/* mathematically the same:
   less stable typically, but useful if n-k+1 < 0 : */
static
double lfastchoose2(double n, double k, int *s_choose)
{
    double r;
    r = lgammafn_sign(n - k + 1., s_choose);
    return lgammafn(n + 1.) - lgammafn(k + 1.) - r;
}
예제 #2
0
파일: lgamma.c 프로젝트: Bgods/r-source
double lgammafn(double x)
{
    return lgammafn_sign(x, NULL);
}