\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(\frac{\frac{\left|\sqrt[3]{1 - 5 \cdot \left(v \cdot v\right)}\right|}{\pi}}{t} \cdot \frac{\sqrt{\sqrt[3]{1 - 5 \cdot \left(v \cdot v\right)}}}{\sqrt{2 \cdot \left(1 - 3 \cdot \left(v \cdot v\right)\right)}}\right) \cdot \frac{\sqrt{1 - 5 \cdot \left(v \cdot v\right)}}{1 - v \cdot v}double code(double v, double t) {
return ((double) (((double) (1.0 - ((double) (5.0 * ((double) (v * v)))))) / ((double) (((double) (((double) (((double) M_PI) * t)) * ((double) sqrt(((double) (2.0 * ((double) (1.0 - ((double) (3.0 * ((double) (v * v)))))))))))) * ((double) (1.0 - ((double) (v * v))))))));
}
double code(double v, double t) {
return ((double) (((double) (((double) (((double) (((double) fabs(((double) cbrt(((double) (1.0 - ((double) (5.0 * ((double) (v * v)))))))))) / ((double) M_PI))) / t)) * ((double) (((double) sqrt(((double) cbrt(((double) (1.0 - ((double) (5.0 * ((double) (v * v)))))))))) / ((double) sqrt(((double) (2.0 * ((double) (1.0 - ((double) (3.0 * ((double) (v * v)))))))))))))) * ((double) (((double) sqrt(((double) (1.0 - ((double) (5.0 * ((double) (v * v)))))))) / ((double) (1.0 - ((double) (v * v))))))));
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
rmApplied add-sqr-sqrt0.5
Applied times-frac0.5
rmApplied add-cube-cbrt0.5
Applied sqrt-prod0.5
Applied times-frac0.5
Simplified0.3
Final simplification0.3
herbie shell --seed 2020155
(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)))))