x - \frac{y \cdot \left(z - t\right)}{a}x - \left(\sqrt[3]{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{a}} \cdot \left(\frac{\sqrt[3]{y}}{\sqrt[3]{a}} \cdot \frac{z - t}{\sqrt[3]{a}}\right)} \cdot \sqrt[3]{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{a}} \cdot \left(\frac{\sqrt[3]{y}}{\sqrt[3]{a}} \cdot \frac{z - t}{\sqrt[3]{a}}\right)}\right) \cdot \sqrt[3]{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{a}} \cdot \left(\frac{\sqrt[3]{y}}{\sqrt[3]{a}} \cdot \frac{z - t}{\sqrt[3]{a}}\right)}double f(double x, double y, double z, double t, double a) {
double r377474 = x;
double r377475 = y;
double r377476 = z;
double r377477 = t;
double r377478 = r377476 - r377477;
double r377479 = r377475 * r377478;
double r377480 = a;
double r377481 = r377479 / r377480;
double r377482 = r377474 - r377481;
return r377482;
}
double f(double x, double y, double z, double t, double a) {
double r377483 = x;
double r377484 = y;
double r377485 = cbrt(r377484);
double r377486 = r377485 * r377485;
double r377487 = a;
double r377488 = cbrt(r377487);
double r377489 = r377486 / r377488;
double r377490 = r377485 / r377488;
double r377491 = z;
double r377492 = t;
double r377493 = r377491 - r377492;
double r377494 = r377493 / r377488;
double r377495 = r377490 * r377494;
double r377496 = r377489 * r377495;
double r377497 = cbrt(r377496);
double r377498 = r377497 * r377497;
double r377499 = r377498 * r377497;
double r377500 = r377483 - r377499;
return r377500;
}




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.7 |
| Herbie | 2.3 |
Initial program 6.2
rmApplied add-cube-cbrt6.6
Applied times-frac3.1
rmApplied add-cube-cbrt3.2
Applied times-frac3.2
Applied associate-*l*2.2
rmApplied add-cube-cbrt2.3
Final simplification2.3
herbie shell --seed 2020047
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, F"
:precision binary64
:herbie-target
(if (< y -1.0761266216389975e-10) (- x (/ 1 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (- x (/ (* y (- z t)) a)) (- x (/ y (/ a (- z t))))))
(- x (/ (* y (- z t)) a)))