https://gcc.gnu.org/bugzilla/show_bug.cgi?id=123409
--- Comment #6 from Infinity Integral <infinityintegral333 at gmail dot com> ---
also replacing line 3246 in (your install
location)/include/c++/13/bits/random.h which is originally
{ _M_param = __param; }
with
{ _M_param = __param;
_M_gd_x.param(std::gamma_distribution<double>::param_type(__param.m() / 2));
_M_gd_y.param(std::gamma_distribution<double>::param_type(__param.n() / 2));}
fixes std::fisher_f_distribution, the example programme (attachment 64876) now
produces the output
first average: 1.23719
second average: 1.43268
second average (new distribution): 1.28694
note that the expected second average is 1.4
