\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\sqrt{n} \cdot \frac{\frac{\sqrt{2} \cdot \sqrt{\pi}}{{\left(n \cdot \left(2 \cdot \pi\right)\right)}^{\left(\frac{k}{2}\right)}}}{\sqrt{k}}(FPCore (k n) :precision binary64 (* (/ 1.0 (sqrt k)) (pow (* (* 2.0 PI) n) (/ (- 1.0 k) 2.0))))
(FPCore (k n) :precision binary64 (* (sqrt n) (/ (/ (* (sqrt 2.0) (sqrt PI)) (pow (* n (* 2.0 PI)) (/ k 2.0))) (sqrt k))))
double code(double k, double n) {
return ((double) ((1.0 / ((double) sqrt(k))) * ((double) pow(((double) (((double) (2.0 * ((double) M_PI))) * n)), (((double) (1.0 - k)) / 2.0)))));
}
double code(double k, double n) {
return ((double) (((double) sqrt(n)) * ((((double) (((double) sqrt(2.0)) * ((double) sqrt(((double) M_PI))))) / ((double) pow(((double) (n * ((double) (2.0 * ((double) M_PI))))), (k / 2.0)))) / ((double) sqrt(k)))));
}



Bits error versus k



Bits error versus n
Results
Initial program 0.5
Simplified0.5
rmApplied div-sub_binary640.5
Applied pow-sub_binary640.4
Simplified0.4
Simplified0.4
rmApplied *-un-lft-identity_binary640.4
Applied sqrt-prod_binary640.4
Applied *-un-lft-identity_binary640.4
Applied sqrt-prod_binary640.6
Applied times-frac_binary640.6
Applied times-frac_binary640.6
Simplified0.6
rmApplied sqrt-prod_binary640.5
Final simplification0.5
herbie shell --seed 2020219
(FPCore (k n)
:name "Migdal et al, Equation (51)"
:precision binary64
(* (/ 1.0 (sqrt k)) (pow (* (* 2.0 PI) n) (/ (- 1.0 k) 2.0))))