x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\left(\left(\sqrt[3]{\sqrt[3]{\frac{y - z}{a - z}}} \cdot \sqrt[3]{\sqrt[3]{\frac{y - z}{a - z}} \cdot \sqrt[3]{\frac{y - z}{a - z}}}\right) \cdot t\right) \cdot \left(\sqrt[3]{\frac{y - z}{a - z}} \cdot \sqrt[3]{\frac{y - z}{a - z}}\right) + \left(x - x \cdot \frac{y - z}{a - z}\right)double f(double x, double y, double z, double t, double a) {
double r24996285 = x;
double r24996286 = y;
double r24996287 = z;
double r24996288 = r24996286 - r24996287;
double r24996289 = t;
double r24996290 = r24996289 - r24996285;
double r24996291 = r24996288 * r24996290;
double r24996292 = a;
double r24996293 = r24996292 - r24996287;
double r24996294 = r24996291 / r24996293;
double r24996295 = r24996285 + r24996294;
return r24996295;
}
double f(double x, double y, double z, double t, double a) {
double r24996296 = y;
double r24996297 = z;
double r24996298 = r24996296 - r24996297;
double r24996299 = a;
double r24996300 = r24996299 - r24996297;
double r24996301 = r24996298 / r24996300;
double r24996302 = cbrt(r24996301);
double r24996303 = cbrt(r24996302);
double r24996304 = r24996302 * r24996302;
double r24996305 = cbrt(r24996304);
double r24996306 = r24996303 * r24996305;
double r24996307 = t;
double r24996308 = r24996306 * r24996307;
double r24996309 = r24996308 * r24996304;
double r24996310 = x;
double r24996311 = r24996310 * r24996301;
double r24996312 = r24996310 - r24996311;
double r24996313 = r24996309 + r24996312;
return r24996313;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.1 |
|---|---|
| Target | 12.1 |
| Herbie | 7.6 |
Initial program 24.1
Simplified11.6
rmApplied fma-udef11.6
rmApplied sub-neg11.6
Applied distribute-lft-in11.6
Applied associate-+l+7.2
Simplified7.2
rmApplied add-cube-cbrt7.5
Applied associate-*l*7.5
rmApplied add-cube-cbrt7.5
Applied cbrt-prod7.6
Final simplification7.6
herbie shell --seed 2019172 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:herbie-target
(if (< z -1.2536131056095036e+188) (- t (* (/ y z) (- t x))) (if (< z 4.446702369113811e+64) (+ x (/ (- y z) (/ (- a z) (- t x)))) (- t (* (/ y z) (- t x)))))
(+ x (/ (* (- y z) (- t x)) (- a z))))