\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\sqrt{\frac{\frac{{\left(\left(\pi \cdot 2\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \cdot 1}{\sqrt{\sqrt{k}}}}{\sqrt{\sqrt{k}}}} \cdot \sqrt{\frac{{\left(\left(\pi \cdot 2\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \cdot 1}{\sqrt{k}}}double f(double k, double n) {
double r5198414 = 1.0;
double r5198415 = k;
double r5198416 = sqrt(r5198415);
double r5198417 = r5198414 / r5198416;
double r5198418 = 2.0;
double r5198419 = atan2(1.0, 0.0);
double r5198420 = r5198418 * r5198419;
double r5198421 = n;
double r5198422 = r5198420 * r5198421;
double r5198423 = r5198414 - r5198415;
double r5198424 = r5198423 / r5198418;
double r5198425 = pow(r5198422, r5198424);
double r5198426 = r5198417 * r5198425;
return r5198426;
}
double f(double k, double n) {
double r5198427 = atan2(1.0, 0.0);
double r5198428 = 2.0;
double r5198429 = r5198427 * r5198428;
double r5198430 = n;
double r5198431 = r5198429 * r5198430;
double r5198432 = 1.0;
double r5198433 = k;
double r5198434 = r5198432 - r5198433;
double r5198435 = r5198434 / r5198428;
double r5198436 = pow(r5198431, r5198435);
double r5198437 = r5198436 * r5198432;
double r5198438 = sqrt(r5198433);
double r5198439 = sqrt(r5198438);
double r5198440 = r5198437 / r5198439;
double r5198441 = r5198440 / r5198439;
double r5198442 = sqrt(r5198441);
double r5198443 = r5198437 / r5198438;
double r5198444 = sqrt(r5198443);
double r5198445 = r5198442 * r5198444;
return r5198445;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.3
rmApplied associate-*l/0.3
rmApplied add-sqr-sqrt0.4
rmApplied add-sqr-sqrt0.4
Applied sqrt-prod0.4
Applied associate-/r*0.4
Final simplification0.4
herbie shell --seed 2019171 +o rules:numerics
(FPCore (k n)
:name "Migdal et al, Equation (51)"
(* (/ 1.0 (sqrt k)) (pow (* (* 2.0 PI) n) (/ (- 1.0 k) 2.0))))