\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 r287398 = 1.0;
double r287399 = 5.0;
double r287400 = v;
double r287401 = r287400 * r287400;
double r287402 = r287399 * r287401;
double r287403 = r287398 - r287402;
double r287404 = atan2(1.0, 0.0);
double r287405 = t;
double r287406 = r287404 * r287405;
double r287407 = 2.0;
double r287408 = 3.0;
double r287409 = r287408 * r287401;
double r287410 = r287398 - r287409;
double r287411 = r287407 * r287410;
double r287412 = sqrt(r287411);
double r287413 = r287406 * r287412;
double r287414 = r287398 - r287401;
double r287415 = r287413 * r287414;
double r287416 = r287403 / r287415;
return r287416;
}
double f(double v, double t) {
double r287417 = 1.0;
double r287418 = 5.0;
double r287419 = v;
double r287420 = r287419 * r287419;
double r287421 = r287418 * r287420;
double r287422 = r287417 - r287421;
double r287423 = atan2(1.0, 0.0);
double r287424 = t;
double r287425 = r287423 * r287424;
double r287426 = 2.0;
double r287427 = 3.0;
double r287428 = pow(r287417, r287427);
double r287429 = 3.0;
double r287430 = r287429 * r287420;
double r287431 = pow(r287430, r287427);
double r287432 = r287428 - r287431;
double r287433 = r287426 * r287432;
double r287434 = sqrt(r287433);
double r287435 = r287425 * r287434;
double r287436 = r287417 - r287420;
double r287437 = r287435 * r287436;
double r287438 = r287422 / r287437;
double r287439 = r287417 * r287417;
double r287440 = r287430 * r287430;
double r287441 = r287417 * r287430;
double r287442 = r287440 + r287441;
double r287443 = r287439 + r287442;
double r287444 = sqrt(r287443);
double r287445 = r287438 * r287444;
return r287445;
}



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