\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 r250879 = 1.0;
double r250880 = 5.0;
double r250881 = v;
double r250882 = r250881 * r250881;
double r250883 = r250880 * r250882;
double r250884 = r250879 - r250883;
double r250885 = atan2(1.0, 0.0);
double r250886 = t;
double r250887 = r250885 * r250886;
double r250888 = 2.0;
double r250889 = 3.0;
double r250890 = r250889 * r250882;
double r250891 = r250879 - r250890;
double r250892 = r250888 * r250891;
double r250893 = sqrt(r250892);
double r250894 = r250887 * r250893;
double r250895 = r250879 - r250882;
double r250896 = r250894 * r250895;
double r250897 = r250884 / r250896;
return r250897;
}
double f(double v, double t) {
double r250898 = 1.0;
double r250899 = 5.0;
double r250900 = v;
double r250901 = r250900 * r250900;
double r250902 = r250899 * r250901;
double r250903 = r250898 - r250902;
double r250904 = atan2(1.0, 0.0);
double r250905 = t;
double r250906 = r250904 * r250905;
double r250907 = 2.0;
double r250908 = 3.0;
double r250909 = pow(r250898, r250908);
double r250910 = 3.0;
double r250911 = r250910 * r250901;
double r250912 = pow(r250911, r250908);
double r250913 = r250909 - r250912;
double r250914 = r250907 * r250913;
double r250915 = sqrt(r250914);
double r250916 = r250906 * r250915;
double r250917 = r250898 - r250901;
double r250918 = r250916 * r250917;
double r250919 = r250903 / r250918;
double r250920 = r250898 * r250898;
double r250921 = r250911 * r250911;
double r250922 = r250898 * r250911;
double r250923 = r250921 + r250922;
double r250924 = r250920 + r250923;
double r250925 = sqrt(r250924);
double r250926 = r250919 * r250925;
return r250926;
}



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 2019353 +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)))))