\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)}1.5 \cdot \frac{{v}^{2}}{t \cdot \left(\sqrt{2} \cdot \left(\sqrt{1} \cdot \pi\right)\right)} + \left(1 \cdot \frac{\sqrt{1}}{t \cdot \left(\sqrt{2} \cdot \pi\right)} - \left(\left(1.5 \cdot \frac{{v}^{4}}{t \cdot \left(\sqrt{2} \cdot \left(\sqrt{1} \cdot \pi\right)\right)} + 1.125 \cdot \frac{{v}^{4}}{t \cdot \left(\sqrt{2} \cdot \left({\left(\sqrt{1}\right)}^{3} \cdot \pi\right)\right)}\right) + 4 \cdot \left(\frac{{v}^{2} \cdot \sqrt{1}}{t \cdot \left(\sqrt{2} \cdot \pi\right)} + \frac{{v}^{4} \cdot \sqrt{1}}{t \cdot \left(\sqrt{2} \cdot \pi\right)}\right)\right)\right)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) (1.5 * ((double) (((double) pow(v, 2.0)) / ((double) (t * ((double) (((double) sqrt(2.0)) * ((double) (((double) sqrt(1.0)) * ((double) M_PI))))))))))) + ((double) (((double) (1.0 * ((double) (((double) sqrt(1.0)) / ((double) (t * ((double) (((double) sqrt(2.0)) * ((double) M_PI))))))))) - ((double) (((double) (((double) (1.5 * ((double) (((double) pow(v, 4.0)) / ((double) (t * ((double) (((double) sqrt(2.0)) * ((double) (((double) sqrt(1.0)) * ((double) M_PI))))))))))) + ((double) (1.125 * ((double) (((double) pow(v, 4.0)) / ((double) (t * ((double) (((double) sqrt(2.0)) * ((double) (((double) pow(((double) sqrt(1.0)), 3.0)) * ((double) M_PI))))))))))))) + ((double) (4.0 * ((double) (((double) (((double) (((double) pow(v, 2.0)) * ((double) sqrt(1.0)))) / ((double) (t * ((double) (((double) sqrt(2.0)) * ((double) M_PI))))))) + ((double) (((double) (((double) pow(v, 4.0)) * ((double) sqrt(1.0)))) / ((double) (t * ((double) (((double) sqrt(2.0)) * ((double) M_PI)))))))))))))))));
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
Taylor expanded around 0 0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020174
(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)))))