\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{1 - 5 \cdot \left(v \cdot v\right)}{\left(\left(\pi \cdot t\right) \cdot \sqrt{2 \cdot \left({1}^{3} - {\left(3 \cdot \left(v \cdot v\right)\right)}^{3}\right)}\right) \cdot \left(1 - v \cdot v\right)} \cdot \sqrt{1 \cdot 1 + \left(\left(3 \cdot \left(v \cdot v\right)\right) \cdot \left(3 \cdot \left(v \cdot v\right)\right) + 1 \cdot \left(3 \cdot \left(v \cdot v\right)\right)\right)}double f(double v, double t) {
double r294854 = 1.0;
double r294855 = 5.0;
double r294856 = v;
double r294857 = r294856 * r294856;
double r294858 = r294855 * r294857;
double r294859 = r294854 - r294858;
double r294860 = atan2(1.0, 0.0);
double r294861 = t;
double r294862 = r294860 * r294861;
double r294863 = 2.0;
double r294864 = 3.0;
double r294865 = r294864 * r294857;
double r294866 = r294854 - r294865;
double r294867 = r294863 * r294866;
double r294868 = sqrt(r294867);
double r294869 = r294862 * r294868;
double r294870 = r294854 - r294857;
double r294871 = r294869 * r294870;
double r294872 = r294859 / r294871;
return r294872;
}
double f(double v, double t) {
double r294873 = 1.0;
double r294874 = 5.0;
double r294875 = v;
double r294876 = r294875 * r294875;
double r294877 = r294874 * r294876;
double r294878 = r294873 - r294877;
double r294879 = atan2(1.0, 0.0);
double r294880 = t;
double r294881 = r294879 * r294880;
double r294882 = 2.0;
double r294883 = 3.0;
double r294884 = pow(r294873, r294883);
double r294885 = 3.0;
double r294886 = r294885 * r294876;
double r294887 = pow(r294886, r294883);
double r294888 = r294884 - r294887;
double r294889 = r294882 * r294888;
double r294890 = sqrt(r294889);
double r294891 = r294881 * r294890;
double r294892 = r294873 - r294876;
double r294893 = r294891 * r294892;
double r294894 = r294878 / r294893;
double r294895 = r294873 * r294873;
double r294896 = r294886 * r294886;
double r294897 = r294873 * r294886;
double r294898 = r294896 + r294897;
double r294899 = r294895 + r294898;
double r294900 = sqrt(r294899);
double r294901 = r294894 * r294900;
return r294901;
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
rmApplied flip3--0.4
Applied associate-*r/0.4
Applied sqrt-div0.4
Applied associate-*r/0.4
Applied associate-*l/0.4
Applied associate-/r/0.4
Final simplification0.4
herbie shell --seed 2020047 +o rules:numerics
(FPCore (v t)
:name "Falkner and Boettcher, Equation (20:1,3)"
:precision binary64
(/ (- 1 (* 5 (* v v))) (* (* (* PI t) (sqrt (* 2 (- 1 (* 3 (* v v)))))) (- 1 (* v v)))))