The R Project SVN R

Rev

Rev 79203 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 79203 Rev 88401
Line 98... Line 98...
98 
	    return(M_LN_SQRT_2PI + (x - 0.5) * log(x) - x);
 98 
	    return(M_LN_SQRT_2PI + (x - 0.5) * log(x) - x);
99 
	else
 99 
	else
100 
#endif
 100 
#endif
101 
	    return M_LN_SQRT_2PI + (x - 0.5) * log(x) - x + lgammacor(x);
 101 
	    return M_LN_SQRT_2PI + (x - 0.5) * log(x) - x + lgammacor(x);
102 
    }
 102 
    }
103 
    /* else: x < -10; y = -x */
 103 
    /* else: x < -10; y = -x > 10 */
104 
    sinpiy = fabs(sinpi(y));
 104 
    sinpiy = fabs(sinpi(y));
105 
 
 105 
 
106 
    if (sinpiy == 0) { /* Negative integer argument ===
 106 
    if (sinpiy == 0) { /* Negative integer argument ===
107 
			  Now UNNECESSARY: caught above */
 107 
			  Now UNNECESSARY: caught above */
108 
	MATHLIB_WARNING(" ** should NEVER happen! *** [lgamma.c: Neg.int, y=%g]\n",y);
 108 
	MATHLIB_WARNING(" ** should NEVER happen! *** [lgamma.c: Neg.int, y=%g]\n",y);