x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;a \le -2.145869478761536680500210929593317634387 \cdot 10^{-21} \lor \neg \left(a \le 5.909250579364137682050418316946687031018 \cdot 10^{-179}\right):\\
\;\;\;\;\frac{t - x}{\frac{a - z}{y - z}} + x\\
\mathbf{else}:\\
\;\;\;\;t + \left(\frac{y \cdot x}{z} - \frac{t}{\frac{z}{y}}\right)\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r514406 = x;
double r514407 = y;
double r514408 = z;
double r514409 = r514407 - r514408;
double r514410 = t;
double r514411 = r514410 - r514406;
double r514412 = r514409 * r514411;
double r514413 = a;
double r514414 = r514413 - r514408;
double r514415 = r514412 / r514414;
double r514416 = r514406 + r514415;
return r514416;
}
double f(double x, double y, double z, double t, double a) {
double r514417 = a;
double r514418 = -2.1458694787615367e-21;
bool r514419 = r514417 <= r514418;
double r514420 = 5.909250579364138e-179;
bool r514421 = r514417 <= r514420;
double r514422 = !r514421;
bool r514423 = r514419 || r514422;
double r514424 = t;
double r514425 = x;
double r514426 = r514424 - r514425;
double r514427 = z;
double r514428 = r514417 - r514427;
double r514429 = y;
double r514430 = r514429 - r514427;
double r514431 = r514428 / r514430;
double r514432 = r514426 / r514431;
double r514433 = r514432 + r514425;
double r514434 = r514429 * r514425;
double r514435 = r514434 / r514427;
double r514436 = r514427 / r514429;
double r514437 = r514424 / r514436;
double r514438 = r514435 - r514437;
double r514439 = r514424 + r514438;
double r514440 = r514423 ? r514433 : r514439;
return r514440;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 23.8 |
|---|---|
| Target | 11.4 |
| Herbie | 10.7 |
if a < -2.1458694787615367e-21 or 5.909250579364138e-179 < a Initial program 21.8
Simplified8.1
rmApplied clear-num8.2
rmApplied *-un-lft-identity8.2
Applied add-cube-cbrt8.2
Applied times-frac8.2
Applied associate-*l*8.2
Simplified8.1
if -2.1458694787615367e-21 < a < 5.909250579364138e-179Initial program 28.3
Simplified18.7
Taylor expanded around inf 18.0
Simplified16.4
Final simplification10.7
herbie shell --seed 2019179
(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))))