x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;a \le -2.669212786079050965464533172334442113648 \cdot 10^{-109}:\\
\;\;\;\;x + \frac{y - z}{a - z} \cdot \left(t - x\right)\\
\mathbf{elif}\;a \le 2.177902234001483733637034366006302694435 \cdot 10^{-87}:\\
\;\;\;\;\left(t + \frac{x}{\frac{z}{y}}\right) - \frac{t}{\frac{z}{y}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y - z}{a - z} \cdot \left(t - x\right)\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r27110431 = x;
double r27110432 = y;
double r27110433 = z;
double r27110434 = r27110432 - r27110433;
double r27110435 = t;
double r27110436 = r27110435 - r27110431;
double r27110437 = r27110434 * r27110436;
double r27110438 = a;
double r27110439 = r27110438 - r27110433;
double r27110440 = r27110437 / r27110439;
double r27110441 = r27110431 + r27110440;
return r27110441;
}
double f(double x, double y, double z, double t, double a) {
double r27110442 = a;
double r27110443 = -2.669212786079051e-109;
bool r27110444 = r27110442 <= r27110443;
double r27110445 = x;
double r27110446 = y;
double r27110447 = z;
double r27110448 = r27110446 - r27110447;
double r27110449 = r27110442 - r27110447;
double r27110450 = r27110448 / r27110449;
double r27110451 = t;
double r27110452 = r27110451 - r27110445;
double r27110453 = r27110450 * r27110452;
double r27110454 = r27110445 + r27110453;
double r27110455 = 2.1779022340014837e-87;
bool r27110456 = r27110442 <= r27110455;
double r27110457 = r27110447 / r27110446;
double r27110458 = r27110445 / r27110457;
double r27110459 = r27110451 + r27110458;
double r27110460 = r27110451 / r27110457;
double r27110461 = r27110459 - r27110460;
double r27110462 = r27110456 ? r27110461 : r27110454;
double r27110463 = r27110444 ? r27110454 : r27110462;
return r27110463;
}




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 | 9.0 |
if a < -2.669212786079051e-109 or 2.1779022340014837e-87 < a Initial program 22.0
rmApplied associate-/l*10.5
rmApplied associate-/r/8.0
if -2.669212786079051e-109 < a < 2.1779022340014837e-87Initial program 28.8
rmApplied add-cube-cbrt29.5
Applied times-frac21.2
Taylor expanded around inf 15.9
Simplified11.2
Final simplification9.0
herbie shell --seed 2019172
(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))))