The R Project SVN R

Rev

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

Rev 42303 Rev 42549
Line 79... Line 79...
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, "beta");
 83 
	ML_ERROR(ME_UNDERFLOW, "beta");
84 
	return ML_UNDERFLOW;
 84 
	/* return ML_UNDERFLOW; pointless giving incorrect value */
85 
    }
 85 
    }
86 
    return exp(val);
 86 
    return exp(val);
87 
}
 87 
}