\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({2}^{\left(\frac{1 - k}{2}\right)} \cdot {\pi}^{\left(\frac{1 - k}{2}\right)}\right)\right) \cdot {n}^{\left(\frac{1 - k}{2}\right)}double f(double k, double n) {
double r65039 = 1.0;
double r65040 = k;
double r65041 = sqrt(r65040);
double r65042 = r65039 / r65041;
double r65043 = 2.0;
double r65044 = atan2(1.0, 0.0);
double r65045 = r65043 * r65044;
double r65046 = n;
double r65047 = r65045 * r65046;
double r65048 = r65039 - r65040;
double r65049 = r65048 / r65043;
double r65050 = pow(r65047, r65049);
double r65051 = r65042 * r65050;
return r65051;
}
double f(double k, double n) {
double r65052 = 1.0;
double r65053 = k;
double r65054 = sqrt(r65053);
double r65055 = r65052 / r65054;
double r65056 = 2.0;
double r65057 = r65052 - r65053;
double r65058 = r65057 / r65056;
double r65059 = pow(r65056, r65058);
double r65060 = atan2(1.0, 0.0);
double r65061 = pow(r65060, r65058);
double r65062 = r65059 * r65061;
double r65063 = r65055 * r65062;
double r65064 = n;
double r65065 = pow(r65064, r65058);
double r65066 = r65063 * r65065;
return r65066;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
rmApplied unpow-prod-down0.6
Applied associate-*r*0.6
rmApplied unpow-prod-down0.6
Final simplification0.6
herbie shell --seed 2019326
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))