\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)}\sqrt{1 \cdot 1 + \left(\left(3 \cdot \left(v \cdot v\right)\right) \cdot \left(3 \cdot \left(v \cdot v\right)\right) + \left(3 \cdot \left(v \cdot v\right)\right) \cdot 1\right)} \cdot \frac{\frac{\frac{1 - \left(v \cdot v\right) \cdot 5}{1 - v \cdot v}}{\pi \cdot \sqrt{\left(\left(1 \cdot 1\right) \cdot 1 - \left(\left(3 \cdot \left(v \cdot v\right)\right) \cdot \left(3 \cdot \left(v \cdot v\right)\right)\right) \cdot \left(3 \cdot \left(v \cdot v\right)\right)\right) \cdot 2}}}{t}double f(double v, double t) {
double r8062245 = 1.0;
double r8062246 = 5.0;
double r8062247 = v;
double r8062248 = r8062247 * r8062247;
double r8062249 = r8062246 * r8062248;
double r8062250 = r8062245 - r8062249;
double r8062251 = atan2(1.0, 0.0);
double r8062252 = t;
double r8062253 = r8062251 * r8062252;
double r8062254 = 2.0;
double r8062255 = 3.0;
double r8062256 = r8062255 * r8062248;
double r8062257 = r8062245 - r8062256;
double r8062258 = r8062254 * r8062257;
double r8062259 = sqrt(r8062258);
double r8062260 = r8062253 * r8062259;
double r8062261 = r8062245 - r8062248;
double r8062262 = r8062260 * r8062261;
double r8062263 = r8062250 / r8062262;
return r8062263;
}
double f(double v, double t) {
double r8062264 = 1.0;
double r8062265 = r8062264 * r8062264;
double r8062266 = 3.0;
double r8062267 = v;
double r8062268 = r8062267 * r8062267;
double r8062269 = r8062266 * r8062268;
double r8062270 = r8062269 * r8062269;
double r8062271 = r8062269 * r8062264;
double r8062272 = r8062270 + r8062271;
double r8062273 = r8062265 + r8062272;
double r8062274 = sqrt(r8062273);
double r8062275 = 5.0;
double r8062276 = r8062268 * r8062275;
double r8062277 = r8062264 - r8062276;
double r8062278 = r8062264 - r8062268;
double r8062279 = r8062277 / r8062278;
double r8062280 = atan2(1.0, 0.0);
double r8062281 = r8062265 * r8062264;
double r8062282 = r8062270 * r8062269;
double r8062283 = r8062281 - r8062282;
double r8062284 = 2.0;
double r8062285 = r8062283 * r8062284;
double r8062286 = sqrt(r8062285);
double r8062287 = r8062280 * r8062286;
double r8062288 = r8062279 / r8062287;
double r8062289 = t;
double r8062290 = r8062288 / r8062289;
double r8062291 = r8062274 * r8062290;
return r8062291;
}



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
Simplified0.3
rmApplied associate-/r*0.1
Final simplification0.1
herbie shell --seed 2019171 +o rules:numerics
(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)))))