The R Project SVN R

Rev

Rev 36820 | Rev 42302 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 36820 Rev 37461
Line 78... Line 78...
78 
	return gammafn(a) * gammafn(b) / gammafn(a+b);
 78 
	return gammafn(a) * gammafn(b) / gammafn(a+b);
79 
 
 79 
 
80 
    val = lbeta(a, b);
 80 
    val = lbeta(a, b);
81 
    if (val < lnsml) {
 81 
    if (val < lnsml) {
82 
	/* a and/or b so big that beta underflows */
 82 
	/* a and/or b so big that beta underflows */
83 
	ML_ERROR(ME_UNDERFLOW);
 83 
	ML_ERROR(ME_UNDERFLOW, "beta");
84 
	return ML_UNDERFLOW;
 84 
	return ML_UNDERFLOW;
85 
    }
 85 
    }
86 
    return exp(val);
 86 
    return exp(val);
87 
}
 87 
}