\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}{\left(2 \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \cdot \left(\sqrt{{\pi}^{\left(\frac{1 - k}{2}\right)}} \cdot \frac{1 \cdot \sqrt{{\pi}^{\left(\frac{1 - k}{2}\right)}}}{\sqrt{k}}\right)double f(double k, double n) {
double r132312 = 1.0;
double r132313 = k;
double r132314 = sqrt(r132313);
double r132315 = r132312 / r132314;
double r132316 = 2.0;
double r132317 = atan2(1.0, 0.0);
double r132318 = r132316 * r132317;
double r132319 = n;
double r132320 = r132318 * r132319;
double r132321 = r132312 - r132313;
double r132322 = r132321 / r132316;
double r132323 = pow(r132320, r132322);
double r132324 = r132315 * r132323;
return r132324;
}
double f(double k, double n) {
double r132325 = 2.0;
double r132326 = n;
double r132327 = r132325 * r132326;
double r132328 = 1.0;
double r132329 = k;
double r132330 = r132328 - r132329;
double r132331 = r132330 / r132325;
double r132332 = pow(r132327, r132331);
double r132333 = atan2(1.0, 0.0);
double r132334 = pow(r132333, r132331);
double r132335 = sqrt(r132334);
double r132336 = r132328 * r132335;
double r132337 = sqrt(r132329);
double r132338 = r132336 / r132337;
double r132339 = r132335 * r132338;
double r132340 = r132332 * r132339;
return r132340;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
Simplified0.4
rmApplied unpow-prod-down0.5
Applied associate-*r*0.5
rmApplied add-sqr-sqrt0.4
Applied associate-*r*0.5
Simplified0.4
Final simplification0.4
herbie shell --seed 2019174
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1.0 (sqrt k)) (pow (* (* 2.0 PI) n) (/ (- 1.0 k) 2.0))))