\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\frac{1 \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1}{2}\right)}}{\sqrt{k} \cdot \left({\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{\frac{k}{2}}{2}\right)} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{\frac{k}{2}}{2}\right)}\right)}double f(double k, double n) {
double r114839 = 1.0;
double r114840 = k;
double r114841 = sqrt(r114840);
double r114842 = r114839 / r114841;
double r114843 = 2.0;
double r114844 = atan2(1.0, 0.0);
double r114845 = r114843 * r114844;
double r114846 = n;
double r114847 = r114845 * r114846;
double r114848 = r114839 - r114840;
double r114849 = r114848 / r114843;
double r114850 = pow(r114847, r114849);
double r114851 = r114842 * r114850;
return r114851;
}
double f(double k, double n) {
double r114852 = 1.0;
double r114853 = 2.0;
double r114854 = atan2(1.0, 0.0);
double r114855 = r114853 * r114854;
double r114856 = n;
double r114857 = r114855 * r114856;
double r114858 = r114852 / r114853;
double r114859 = pow(r114857, r114858);
double r114860 = r114852 * r114859;
double r114861 = k;
double r114862 = sqrt(r114861);
double r114863 = r114861 / r114853;
double r114864 = 2.0;
double r114865 = r114863 / r114864;
double r114866 = pow(r114857, r114865);
double r114867 = r114866 * r114866;
double r114868 = r114862 * r114867;
double r114869 = r114860 / r114868;
return r114869;
}



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 sqr-pow0.4
Final simplification0.4
herbie shell --seed 2020057
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))