x + \frac{y \cdot \left(z - t\right)}{a}x + \frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\frac{\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}}{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{a} \cdot \sqrt[3]{a}}}} \cdot \frac{\sqrt[3]{1}}{\frac{\frac{\sqrt[3]{1}}{\sqrt[3]{z - t}}}{\frac{\sqrt[3]{y}}{\sqrt[3]{a}}}}double f(double x, double y, double z, double t, double a) {
double r254456 = x;
double r254457 = y;
double r254458 = z;
double r254459 = t;
double r254460 = r254458 - r254459;
double r254461 = r254457 * r254460;
double r254462 = a;
double r254463 = r254461 / r254462;
double r254464 = r254456 + r254463;
return r254464;
}
double f(double x, double y, double z, double t, double a) {
double r254465 = x;
double r254466 = 1.0;
double r254467 = cbrt(r254466);
double r254468 = r254467 * r254467;
double r254469 = z;
double r254470 = t;
double r254471 = r254469 - r254470;
double r254472 = cbrt(r254471);
double r254473 = r254472 * r254472;
double r254474 = r254468 / r254473;
double r254475 = y;
double r254476 = cbrt(r254475);
double r254477 = r254476 * r254476;
double r254478 = a;
double r254479 = cbrt(r254478);
double r254480 = r254479 * r254479;
double r254481 = r254477 / r254480;
double r254482 = r254474 / r254481;
double r254483 = r254468 / r254482;
double r254484 = r254467 / r254472;
double r254485 = r254476 / r254479;
double r254486 = r254484 / r254485;
double r254487 = r254467 / r254486;
double r254488 = r254483 * r254487;
double r254489 = r254465 + r254488;
return r254489;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.2 |
|---|---|
| Target | 0.5 |
| Herbie | 0.9 |
Initial program 6.2
rmApplied associate-/l*5.7
rmApplied div-inv5.7
Applied associate-/r*2.4
rmApplied *-un-lft-identity2.4
Applied *-un-lft-identity2.4
Applied times-frac2.4
Applied associate-/l*2.4
rmApplied add-cube-cbrt2.8
Applied add-cube-cbrt3.0
Applied times-frac3.0
Applied add-cube-cbrt3.0
Applied add-cube-cbrt3.0
Applied times-frac3.1
Applied times-frac0.9
Applied add-cube-cbrt0.9
Applied times-frac0.9
Simplified0.9
Simplified0.9
Final simplification0.9
herbie shell --seed 2019294
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, E"
:precision binary64
:herbie-target
(if (< y -1.07612662163899753e-10) (+ x (/ 1 (/ (/ a (- z t)) y))) (if (< y 2.8944268627920891e-49) (+ x (/ (* y (- z t)) a)) (+ x (/ y (/ a (- z t))))))
(+ x (/ (* y (- z t)) a)))