\frac{1 - 5 \cdot \left(v \cdot v\right)}{\left(\left(\pi \cdot t\right) \cdot \sqrt{2 \cdot \left(1 - 3 \cdot \left(v \cdot v\right)\right)}\right) \cdot \left(1 - v \cdot v\right)}\left(1 - 5 \cdot \left(v \cdot v\right)\right) \cdot \left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{1 - v \cdot v} \cdot \frac{\frac{\sqrt[3]{1}}{\sqrt{2 - \left(v \cdot v\right) \cdot 6}}}{t \cdot \pi}\right)(FPCore (v t) :precision binary64 (/ (- 1.0 (* 5.0 (* v v))) (* (* (* PI t) (sqrt (* 2.0 (- 1.0 (* 3.0 (* v v)))))) (- 1.0 (* v v)))))
(FPCore (v t) :precision binary64 (* (- 1.0 (* 5.0 (* v v))) (* (/ (* (cbrt 1.0) (cbrt 1.0)) (- 1.0 (* v v))) (/ (/ (cbrt 1.0) (sqrt (- 2.0 (* (* v v) 6.0)))) (* t PI)))))
double code(double v, double t) {
return (1.0 - (5.0 * (v * v))) / (((((double) M_PI) * t) * sqrt(2.0 * (1.0 - (3.0 * (v * v))))) * (1.0 - (v * v)));
}
double code(double v, double t) {
return (1.0 - (5.0 * (v * v))) * (((cbrt(1.0) * cbrt(1.0)) / (1.0 - (v * v))) * ((cbrt(1.0) / sqrt(2.0 - ((v * v) * 6.0))) / (t * ((double) M_PI))));
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
Simplified0.4
rmApplied div-inv_binary64_21210.4
Simplified0.4
rmApplied *-un-lft-identity_binary64_21240.4
Applied sqrt-prod_binary64_21400.4
Applied add-cube-cbrt_binary64_21590.4
Applied times-frac_binary64_21300.4
Applied times-frac_binary64_21300.4
Final simplification0.4
herbie shell --seed 2021015
(FPCore (v t)
:name "Falkner and Boettcher, Equation (20:1,3)"
:precision binary64
(/ (- 1.0 (* 5.0 (* v v))) (* (* (* PI t) (sqrt (* 2.0 (- 1.0 (* 3.0 (* v v)))))) (- 1.0 (* v v)))))