\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 r328030 = 1.0;
double r328031 = 5.0;
double r328032 = v;
double r328033 = r328032 * r328032;
double r328034 = r328031 * r328033;
double r328035 = r328030 - r328034;
double r328036 = atan2(1.0, 0.0);
double r328037 = t;
double r328038 = r328036 * r328037;
double r328039 = 2.0;
double r328040 = 3.0;
double r328041 = r328040 * r328033;
double r328042 = r328030 - r328041;
double r328043 = r328039 * r328042;
double r328044 = sqrt(r328043);
double r328045 = r328038 * r328044;
double r328046 = r328030 - r328033;
double r328047 = r328045 * r328046;
double r328048 = r328035 / r328047;
return r328048;
}
double f(double v, double t) {
double r328049 = 1.0;
double r328050 = 5.0;
double r328051 = v;
double r328052 = r328051 * r328051;
double r328053 = r328050 * r328052;
double r328054 = r328049 - r328053;
double r328055 = atan2(1.0, 0.0);
double r328056 = t;
double r328057 = r328055 * r328056;
double r328058 = 2.0;
double r328059 = 3.0;
double r328060 = pow(r328049, r328059);
double r328061 = 3.0;
double r328062 = r328061 * r328052;
double r328063 = pow(r328062, r328059);
double r328064 = r328060 - r328063;
double r328065 = r328058 * r328064;
double r328066 = sqrt(r328065);
double r328067 = r328057 * r328066;
double r328068 = r328049 - r328052;
double r328069 = r328067 * r328068;
double r328070 = r328054 / r328069;
double r328071 = r328049 * r328049;
double r328072 = r328062 * r328062;
double r328073 = r328049 * r328062;
double r328074 = r328072 + r328073;
double r328075 = r328071 + r328074;
double r328076 = sqrt(r328075);
double r328077 = r328070 * r328076;
return r328077;
}



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