\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\left(\left(\frac{1}{\sqrt{k}} \cdot {2}^{\left(\frac{1 - k}{2}\right)}\right) \cdot {\pi}^{\left(\frac{1 - k}{2}\right)}\right) \cdot {n}^{\left(\frac{1 - k}{2}\right)}double f(double k, double n) {
double r8018614 = 1.0;
double r8018615 = k;
double r8018616 = sqrt(r8018615);
double r8018617 = r8018614 / r8018616;
double r8018618 = 2.0;
double r8018619 = atan2(1.0, 0.0);
double r8018620 = r8018618 * r8018619;
double r8018621 = n;
double r8018622 = r8018620 * r8018621;
double r8018623 = r8018614 - r8018615;
double r8018624 = r8018623 / r8018618;
double r8018625 = pow(r8018622, r8018624);
double r8018626 = r8018617 * r8018625;
return r8018626;
}
double f(double k, double n) {
double r8018627 = 1.0;
double r8018628 = k;
double r8018629 = sqrt(r8018628);
double r8018630 = r8018627 / r8018629;
double r8018631 = 2.0;
double r8018632 = r8018627 - r8018628;
double r8018633 = r8018632 / r8018631;
double r8018634 = pow(r8018631, r8018633);
double r8018635 = r8018630 * r8018634;
double r8018636 = atan2(1.0, 0.0);
double r8018637 = pow(r8018636, r8018633);
double r8018638 = r8018635 * r8018637;
double r8018639 = n;
double r8018640 = pow(r8018639, r8018633);
double r8018641 = r8018638 * r8018640;
return r8018641;
}



Bits error versus k



Bits error versus n
Results
Initial program 0.4
rmApplied unpow-prod-down0.5
Applied associate-*r*0.5
rmApplied unpow-prod-down0.4
Applied associate-*r*0.4
Final simplification0.4
herbie shell --seed 2019174 +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))))