http://gcc.gnu.org/bugzilla/show_bug.cgi?id=48114
--- Comment #12 from Andrey Zholos <aaz at althenia dot net> 2011-03-14 18:07:22 UTC --- The double __p also fixes an infinite (maybe) loop in both binomial_distribution(100, .75) and binomial_distribution(100, .25) when _GLIBCXX_USE_C99_MATH_TR1 is defined. I guess the proper way to test the algorithms is to generate a lot of values and do a Pearson's chi-square test to check that the distribution is correct. Then that can be used to make a faster testcase for just a few values.