\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)}\frac{\sqrt{1 - \left(v \cdot v\right) \cdot 5}}{1 - v \cdot v} \cdot \left(\frac{\frac{\frac{\sqrt{1 - \left(v \cdot v\right) \cdot 5}}{\sqrt{\left(1 + 3 \cdot \left(v \cdot v\right)\right) \cdot \left(2 \cdot \left(1 - 3 \cdot \left(v \cdot v\right)\right)\right)}}}{\pi}}{t} \cdot \sqrt{1 + 3 \cdot \left(v \cdot v\right)}\right)double f(double v, double t) {
double r7857779 = 1.0;
double r7857780 = 5.0;
double r7857781 = v;
double r7857782 = r7857781 * r7857781;
double r7857783 = r7857780 * r7857782;
double r7857784 = r7857779 - r7857783;
double r7857785 = atan2(1.0, 0.0);
double r7857786 = t;
double r7857787 = r7857785 * r7857786;
double r7857788 = 2.0;
double r7857789 = 3.0;
double r7857790 = r7857789 * r7857782;
double r7857791 = r7857779 - r7857790;
double r7857792 = r7857788 * r7857791;
double r7857793 = sqrt(r7857792);
double r7857794 = r7857787 * r7857793;
double r7857795 = r7857779 - r7857782;
double r7857796 = r7857794 * r7857795;
double r7857797 = r7857784 / r7857796;
return r7857797;
}
double f(double v, double t) {
double r7857798 = 1.0;
double r7857799 = v;
double r7857800 = r7857799 * r7857799;
double r7857801 = 5.0;
double r7857802 = r7857800 * r7857801;
double r7857803 = r7857798 - r7857802;
double r7857804 = sqrt(r7857803);
double r7857805 = r7857798 - r7857800;
double r7857806 = r7857804 / r7857805;
double r7857807 = 3.0;
double r7857808 = r7857807 * r7857800;
double r7857809 = r7857798 + r7857808;
double r7857810 = 2.0;
double r7857811 = r7857798 - r7857808;
double r7857812 = r7857810 * r7857811;
double r7857813 = r7857809 * r7857812;
double r7857814 = sqrt(r7857813);
double r7857815 = r7857804 / r7857814;
double r7857816 = atan2(1.0, 0.0);
double r7857817 = r7857815 / r7857816;
double r7857818 = t;
double r7857819 = r7857817 / r7857818;
double r7857820 = sqrt(r7857809);
double r7857821 = r7857819 * r7857820;
double r7857822 = r7857806 * r7857821;
return r7857822;
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
rmApplied add-sqr-sqrt0.4
Applied times-frac0.4
rmApplied *-un-lft-identity0.4
Applied sqrt-prod0.4
Applied times-frac0.5
Simplified0.5
rmApplied flip--0.5
Applied associate-*r/0.5
Applied sqrt-div0.5
Applied associate-/r/0.5
Applied associate-*r*0.5
Simplified0.1
Final simplification0.1
herbie shell --seed 2019200
(FPCore (v t)
:name "Falkner and Boettcher, Equation (20:1,3)"
(/ (- 1.0 (* 5.0 (* v v))) (* (* (* PI t) (sqrt (* 2.0 (- 1.0 (* 3.0 (* v v)))))) (- 1.0 (* v v)))))