id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
5433,beta function: when b<epsilon return tgamma(b)?,Eric Butter <egbutter@…>,John Maddock,"In the beta.hpp header for the beta function (i.e., *not* the beta distribution) around line 131, you have: 

   if((c == a) && (b < tools::epsilon<T>()))
      return boost::math::tgamma(b, pol);
   else if((c == b) && (a < tools::epsilon<T>()))
      return boost::math::tgamma(a, pol);

Don't you actually want if b<epsilon, tgamma(a)?  beta(a,b) reduces to (approximately) gamma(a)b^-a for a>>b, which seems like the case for which you are trying to account?

Cheers,  

Eric ",Bugs,closed,To Be Determined,math,Boost 1.47.0,Problem,invalid,,