\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 \cdot 1 - \left(3 \cdot \left(v \cdot v\right)\right) \cdot \left(3 \cdot \left(v \cdot v\right)\right)\right)}\right) \cdot \left(1 \cdot 1 - \left(v \cdot v\right) \cdot \left(v \cdot v\right)\right)} \cdot \left(\sqrt{1 + 3 \cdot \left(v \cdot v\right)} \cdot \left(1 + v \cdot v\right)\right)double f(double v, double t) {
double r10312070 = 1.0;
double r10312071 = 5.0;
double r10312072 = v;
double r10312073 = r10312072 * r10312072;
double r10312074 = r10312071 * r10312073;
double r10312075 = r10312070 - r10312074;
double r10312076 = atan2(1.0, 0.0);
double r10312077 = t;
double r10312078 = r10312076 * r10312077;
double r10312079 = 2.0;
double r10312080 = 3.0;
double r10312081 = r10312080 * r10312073;
double r10312082 = r10312070 - r10312081;
double r10312083 = r10312079 * r10312082;
double r10312084 = sqrt(r10312083);
double r10312085 = r10312078 * r10312084;
double r10312086 = r10312070 - r10312073;
double r10312087 = r10312085 * r10312086;
double r10312088 = r10312075 / r10312087;
return r10312088;
}
double f(double v, double t) {
double r10312089 = 1.0;
double r10312090 = 5.0;
double r10312091 = v;
double r10312092 = r10312091 * r10312091;
double r10312093 = r10312090 * r10312092;
double r10312094 = r10312089 - r10312093;
double r10312095 = atan2(1.0, 0.0);
double r10312096 = t;
double r10312097 = r10312095 * r10312096;
double r10312098 = 2.0;
double r10312099 = r10312089 * r10312089;
double r10312100 = 3.0;
double r10312101 = r10312100 * r10312092;
double r10312102 = r10312101 * r10312101;
double r10312103 = r10312099 - r10312102;
double r10312104 = r10312098 * r10312103;
double r10312105 = sqrt(r10312104);
double r10312106 = r10312097 * r10312105;
double r10312107 = r10312092 * r10312092;
double r10312108 = r10312099 - r10312107;
double r10312109 = r10312106 * r10312108;
double r10312110 = r10312094 / r10312109;
double r10312111 = r10312089 + r10312101;
double r10312112 = sqrt(r10312111);
double r10312113 = r10312089 + r10312092;
double r10312114 = r10312112 * r10312113;
double r10312115 = r10312110 * r10312114;
return r10312115;
}



Bits error versus v



Bits error versus t
Results
Initial program 0.4
rmApplied flip--0.4
Applied flip--0.4
Applied associate-*r/0.4
Applied sqrt-div0.5
Applied associate-*r/0.5
Applied frac-times0.5
Applied associate-/r/0.4
Final simplification0.4
herbie shell --seed 2019174
(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)))))