x + y \cdot \frac{z - t}{z - a}\frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{z - a} \cdot \sqrt[3]{z - a}} \cdot \frac{y}{\frac{\sqrt[3]{z - a}}{\sqrt[3]{z - t}}} + xdouble f(double x, double y, double z, double t, double a) {
double r814448 = x;
double r814449 = y;
double r814450 = z;
double r814451 = t;
double r814452 = r814450 - r814451;
double r814453 = a;
double r814454 = r814450 - r814453;
double r814455 = r814452 / r814454;
double r814456 = r814449 * r814455;
double r814457 = r814448 + r814456;
return r814457;
}
double f(double x, double y, double z, double t, double a) {
double r814458 = z;
double r814459 = t;
double r814460 = r814458 - r814459;
double r814461 = cbrt(r814460);
double r814462 = r814461 * r814461;
double r814463 = a;
double r814464 = r814458 - r814463;
double r814465 = cbrt(r814464);
double r814466 = r814465 * r814465;
double r814467 = r814462 / r814466;
double r814468 = y;
double r814469 = r814465 / r814461;
double r814470 = r814468 / r814469;
double r814471 = r814467 * r814470;
double r814472 = x;
double r814473 = r814471 + r814472;
return r814473;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 1.5 |
|---|---|
| Target | 1.4 |
| Herbie | 0.6 |
Initial program 1.5
Simplified1.5
rmApplied clear-num1.6
rmApplied fma-udef1.6
Simplified1.4
rmApplied add-cube-cbrt1.9
Applied add-cube-cbrt1.8
Applied times-frac1.8
Applied *-un-lft-identity1.8
Applied times-frac0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020060 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisLine from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(+ x (/ y (/ (- z a) (- z t))))
(+ x (* y (/ (- z t) (- z a)))))