\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 r329360 = 1.0;
double r329361 = 5.0;
double r329362 = v;
double r329363 = r329362 * r329362;
double r329364 = r329361 * r329363;
double r329365 = r329360 - r329364;
double r329366 = atan2(1.0, 0.0);
double r329367 = t;
double r329368 = r329366 * r329367;
double r329369 = 2.0;
double r329370 = 3.0;
double r329371 = r329370 * r329363;
double r329372 = r329360 - r329371;
double r329373 = r329369 * r329372;
double r329374 = sqrt(r329373);
double r329375 = r329368 * r329374;
double r329376 = r329360 - r329363;
double r329377 = r329375 * r329376;
double r329378 = r329365 / r329377;
return r329378;
}
double f(double v, double t) {
double r329379 = 1.0;
double r329380 = 5.0;
double r329381 = v;
double r329382 = r329381 * r329381;
double r329383 = r329380 * r329382;
double r329384 = r329379 - r329383;
double r329385 = atan2(1.0, 0.0);
double r329386 = t;
double r329387 = r329385 * r329386;
double r329388 = 2.0;
double r329389 = 3.0;
double r329390 = pow(r329379, r329389);
double r329391 = 3.0;
double r329392 = r329391 * r329382;
double r329393 = pow(r329392, r329389);
double r329394 = r329390 - r329393;
double r329395 = r329388 * r329394;
double r329396 = sqrt(r329395);
double r329397 = r329387 * r329396;
double r329398 = r329379 - r329382;
double r329399 = r329397 * r329398;
double r329400 = r329384 / r329399;
double r329401 = r329379 * r329379;
double r329402 = r329392 * r329392;
double r329403 = r329379 * r329392;
double r329404 = r329402 + r329403;
double r329405 = r329401 + r329404;
double r329406 = sqrt(r329405);
double r329407 = r329400 * r329406;
return r329407;
}



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