x + y \cdot \frac{z - t}{z - a}x + \left(y \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{z - a} \cdot \sqrt[3]{z - a}}\right) \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{z - a}}double f(double x, double y, double z, double t, double a) {
double r508500 = x;
double r508501 = y;
double r508502 = z;
double r508503 = t;
double r508504 = r508502 - r508503;
double r508505 = a;
double r508506 = r508502 - r508505;
double r508507 = r508504 / r508506;
double r508508 = r508501 * r508507;
double r508509 = r508500 + r508508;
return r508509;
}
double f(double x, double y, double z, double t, double a) {
double r508510 = x;
double r508511 = y;
double r508512 = z;
double r508513 = t;
double r508514 = r508512 - r508513;
double r508515 = cbrt(r508514);
double r508516 = r508515 * r508515;
double r508517 = a;
double r508518 = r508512 - r508517;
double r508519 = cbrt(r508518);
double r508520 = r508519 * r508519;
double r508521 = r508516 / r508520;
double r508522 = r508511 * r508521;
double r508523 = r508515 / r508519;
double r508524 = r508522 * r508523;
double r508525 = r508510 + r508524;
return r508525;
}




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 | 1.2 |
| Herbie | 0.5 |
Initial program 1.3
rmApplied add-cube-cbrt1.9
Applied add-cube-cbrt1.7
Applied times-frac1.7
Applied associate-*r*0.5
Final simplification0.5
herbie shell --seed 2019208
(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)))))