\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}x + \left(y - \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t}} \cdot \left(\frac{\sqrt[3]{\sqrt[3]{z - t}} \cdot \sqrt[3]{\sqrt[3]{z - t}}}{\sqrt[3]{\sqrt[3]{a - t}} \cdot \sqrt[3]{\sqrt[3]{a - t}}} \cdot \left(\frac{\sqrt[3]{\sqrt[3]{z - t}}}{\sqrt[3]{\sqrt[3]{a - t}}} \cdot \frac{y}{\sqrt[3]{a - t}}\right)\right)\right)double f(double x, double y, double z, double t, double a) {
double r600631 = x;
double r600632 = y;
double r600633 = r600631 + r600632;
double r600634 = z;
double r600635 = t;
double r600636 = r600634 - r600635;
double r600637 = r600636 * r600632;
double r600638 = a;
double r600639 = r600638 - r600635;
double r600640 = r600637 / r600639;
double r600641 = r600633 - r600640;
return r600641;
}
double f(double x, double y, double z, double t, double a) {
double r600642 = x;
double r600643 = y;
double r600644 = z;
double r600645 = t;
double r600646 = r600644 - r600645;
double r600647 = cbrt(r600646);
double r600648 = r600647 * r600647;
double r600649 = a;
double r600650 = r600649 - r600645;
double r600651 = cbrt(r600650);
double r600652 = r600648 / r600651;
double r600653 = cbrt(r600647);
double r600654 = r600653 * r600653;
double r600655 = cbrt(r600651);
double r600656 = r600655 * r600655;
double r600657 = r600654 / r600656;
double r600658 = r600653 / r600655;
double r600659 = r600643 / r600651;
double r600660 = r600658 * r600659;
double r600661 = r600657 * r600660;
double r600662 = r600652 * r600661;
double r600663 = r600643 - r600662;
double r600664 = r600642 + r600663;
return r600664;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 16.2 |
|---|---|
| Target | 8.7 |
| Herbie | 7.8 |
Initial program 16.2
rmApplied add-cube-cbrt16.3
Applied times-frac11.7
rmApplied add-cube-cbrt11.7
Applied times-frac11.7
Applied associate-*l*11.0
rmApplied associate--l+7.7
rmApplied add-cube-cbrt10.4
Applied add-cube-cbrt7.8
Applied times-frac7.8
Applied associate-*l*7.8
Final simplification7.8
herbie shell --seed 2020100
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(if (< (- (+ x y) (/ (* (- z t) y) (- a t))) -1.3664970889390727e-07) (- (+ y x) (* (* (- z t) (/ 1 (- a t))) y)) (if (< (- (+ x y) (/ (* (- z t) y) (- a t))) 1.4754293444577233e-239) (/ (- (* y (- a z)) (* x t)) (- a t)) (- (+ y x) (* (* (- z t) (/ 1 (- a t))) y))))
(- (+ x y) (/ (* (- z t) y) (- a t))))