\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{\frac{{\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(2 \cdot \frac{\frac{1}{2}}{2}\right)}}{\frac{\sqrt{k}}{1}}}{{\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{k}{2}\right)}}double f(double k, double n) {
double r148738 = 1.0;
double r148739 = k;
double r148740 = sqrt(r148739);
double r148741 = r148738 / r148740;
double r148742 = 2.0;
double r148743 = atan2(1.0, 0.0);
double r148744 = r148742 * r148743;
double r148745 = n;
double r148746 = r148744 * r148745;
double r148747 = r148738 - r148739;
double r148748 = r148747 / r148742;
double r148749 = pow(r148746, r148748);
double r148750 = r148741 * r148749;
return r148750;
}
double f(double k, double n) {
double r148751 = 2.0;
double r148752 = atan2(1.0, 0.0);
double r148753 = r148751 * r148752;
double r148754 = n;
double r148755 = r148753 * r148754;
double r148756 = 2.0;
double r148757 = 1.0;
double r148758 = r148757 / r148751;
double r148759 = r148758 / r148756;
double r148760 = r148756 * r148759;
double r148761 = pow(r148755, r148760);
double r148762 = k;
double r148763 = sqrt(r148762);
double r148764 = r148763 / r148757;
double r148765 = r148761 / r148764;
double r148766 = r148762 / r148751;
double r148767 = pow(r148755, r148766);
double r148768 = r148765 / r148767;
return r148768;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.5
rmApplied div-sub0.5
Applied pow-sub0.4
Applied frac-times0.4
rmApplied add-sqr-sqrt0.4
Applied associate-*r*0.4
rmApplied associate-/r*0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020062
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))