x + y \cdot \frac{z - t}{a - t}x + \left(y \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}\right) \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{a - t}}double f(double x, double y, double z, double t, double a) {
double r357622 = x;
double r357623 = y;
double r357624 = z;
double r357625 = t;
double r357626 = r357624 - r357625;
double r357627 = a;
double r357628 = r357627 - r357625;
double r357629 = r357626 / r357628;
double r357630 = r357623 * r357629;
double r357631 = r357622 + r357630;
return r357631;
}
double f(double x, double y, double z, double t, double a) {
double r357632 = x;
double r357633 = y;
double r357634 = z;
double r357635 = t;
double r357636 = r357634 - r357635;
double r357637 = cbrt(r357636);
double r357638 = r357637 * r357637;
double r357639 = a;
double r357640 = r357639 - r357635;
double r357641 = cbrt(r357640);
double r357642 = r357641 * r357641;
double r357643 = r357638 / r357642;
double r357644 = r357633 * r357643;
double r357645 = r357637 / r357641;
double r357646 = r357644 * r357645;
double r357647 = r357632 + r357646;
return r357647;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 1.3 |
|---|---|
| Target | 0.4 |
| Herbie | 0.5 |
Initial program 1.3
rmApplied add-cube-cbrt1.8
Applied add-cube-cbrt1.6
Applied times-frac1.6
Applied associate-*r*0.5
Final simplification0.5
herbie shell --seed 2019347 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisLine from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(if (< y -8.508084860551241e-17) (+ x (* y (/ (- z t) (- a t)))) (if (< y 2.894426862792089e-49) (+ x (* (* y (- z t)) (/ 1 (- a t)))) (+ x (* y (/ (- z t) (- a t))))))
(+ x (* y (/ (- z t) (- a t)))))