\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 r202413 = 1.0;
double r202414 = 5.0;
double r202415 = v;
double r202416 = r202415 * r202415;
double r202417 = r202414 * r202416;
double r202418 = r202413 - r202417;
double r202419 = atan2(1.0, 0.0);
double r202420 = t;
double r202421 = r202419 * r202420;
double r202422 = 2.0;
double r202423 = 3.0;
double r202424 = r202423 * r202416;
double r202425 = r202413 - r202424;
double r202426 = r202422 * r202425;
double r202427 = sqrt(r202426);
double r202428 = r202421 * r202427;
double r202429 = r202413 - r202416;
double r202430 = r202428 * r202429;
double r202431 = r202418 / r202430;
return r202431;
}
double f(double v, double t) {
double r202432 = 1.0;
double r202433 = 5.0;
double r202434 = v;
double r202435 = r202434 * r202434;
double r202436 = r202433 * r202435;
double r202437 = r202432 - r202436;
double r202438 = atan2(1.0, 0.0);
double r202439 = t;
double r202440 = r202438 * r202439;
double r202441 = 2.0;
double r202442 = 3.0;
double r202443 = pow(r202432, r202442);
double r202444 = 3.0;
double r202445 = r202444 * r202435;
double r202446 = pow(r202445, r202442);
double r202447 = r202443 - r202446;
double r202448 = r202441 * r202447;
double r202449 = sqrt(r202448);
double r202450 = r202440 * r202449;
double r202451 = r202432 - r202435;
double r202452 = r202450 * r202451;
double r202453 = r202437 / r202452;
double r202454 = r202432 * r202432;
double r202455 = r202445 * r202445;
double r202456 = r202432 * r202445;
double r202457 = r202455 + r202456;
double r202458 = r202454 + r202457;
double r202459 = sqrt(r202458);
double r202460 = r202453 * r202459;
return r202460;
}



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