\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 r248375 = 1.0;
double r248376 = 5.0;
double r248377 = v;
double r248378 = r248377 * r248377;
double r248379 = r248376 * r248378;
double r248380 = r248375 - r248379;
double r248381 = atan2(1.0, 0.0);
double r248382 = t;
double r248383 = r248381 * r248382;
double r248384 = 2.0;
double r248385 = 3.0;
double r248386 = r248385 * r248378;
double r248387 = r248375 - r248386;
double r248388 = r248384 * r248387;
double r248389 = sqrt(r248388);
double r248390 = r248383 * r248389;
double r248391 = r248375 - r248378;
double r248392 = r248390 * r248391;
double r248393 = r248380 / r248392;
return r248393;
}
double f(double v, double t) {
double r248394 = 1.0;
double r248395 = 5.0;
double r248396 = v;
double r248397 = r248396 * r248396;
double r248398 = r248395 * r248397;
double r248399 = r248394 - r248398;
double r248400 = atan2(1.0, 0.0);
double r248401 = t;
double r248402 = r248400 * r248401;
double r248403 = 2.0;
double r248404 = 3.0;
double r248405 = pow(r248394, r248404);
double r248406 = 3.0;
double r248407 = r248406 * r248397;
double r248408 = pow(r248407, r248404);
double r248409 = r248405 - r248408;
double r248410 = r248403 * r248409;
double r248411 = sqrt(r248410);
double r248412 = r248402 * r248411;
double r248413 = r248394 - r248397;
double r248414 = r248412 * r248413;
double r248415 = r248399 / r248414;
double r248416 = r248394 * r248394;
double r248417 = r248407 * r248407;
double r248418 = r248394 * r248407;
double r248419 = r248417 + r248418;
double r248420 = r248416 + r248419;
double r248421 = sqrt(r248420);
double r248422 = r248415 * r248421;
return r248422;
}



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