\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\left(\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{\frac{1 - k}{2}}{2}\right)}\right) \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{\frac{1 - k}{2}}{2}\right)}double f(double k, double n) {
double r86078 = 1.0;
double r86079 = k;
double r86080 = sqrt(r86079);
double r86081 = r86078 / r86080;
double r86082 = 2.0;
double r86083 = atan2(1.0, 0.0);
double r86084 = r86082 * r86083;
double r86085 = n;
double r86086 = r86084 * r86085;
double r86087 = r86078 - r86079;
double r86088 = r86087 / r86082;
double r86089 = pow(r86086, r86088);
double r86090 = r86081 * r86089;
return r86090;
}
double f(double k, double n) {
double r86091 = 1.0;
double r86092 = k;
double r86093 = sqrt(r86092);
double r86094 = r86091 / r86093;
double r86095 = 2.0;
double r86096 = atan2(1.0, 0.0);
double r86097 = r86095 * r86096;
double r86098 = n;
double r86099 = r86097 * r86098;
double r86100 = r86091 - r86092;
double r86101 = r86100 / r86095;
double r86102 = 2.0;
double r86103 = r86101 / r86102;
double r86104 = pow(r86099, r86103);
double r86105 = r86094 * r86104;
double r86106 = r86105 * r86104;
return r86106;
}



Bits error versus k



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