x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;a \le -1.2559578325599803 \cdot 10^{-148}:\\
\;\;\;\;x + \frac{1}{\frac{a - z}{y - z}} \cdot \left(t - x\right)\\
\mathbf{elif}\;a \le 9.426461174744231 \cdot 10^{-94}:\\
\;\;\;\;\left(\frac{x}{z} - \frac{t}{z}\right) \cdot y + t\\
\mathbf{else}:\\
\;\;\;\;x + \frac{1}{\frac{a - z}{y - z}} \cdot \left(t - x\right)\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r26988363 = x;
double r26988364 = y;
double r26988365 = z;
double r26988366 = r26988364 - r26988365;
double r26988367 = t;
double r26988368 = r26988367 - r26988363;
double r26988369 = r26988366 * r26988368;
double r26988370 = a;
double r26988371 = r26988370 - r26988365;
double r26988372 = r26988369 / r26988371;
double r26988373 = r26988363 + r26988372;
return r26988373;
}
double f(double x, double y, double z, double t, double a) {
double r26988374 = a;
double r26988375 = -1.2559578325599803e-148;
bool r26988376 = r26988374 <= r26988375;
double r26988377 = x;
double r26988378 = 1.0;
double r26988379 = z;
double r26988380 = r26988374 - r26988379;
double r26988381 = y;
double r26988382 = r26988381 - r26988379;
double r26988383 = r26988380 / r26988382;
double r26988384 = r26988378 / r26988383;
double r26988385 = t;
double r26988386 = r26988385 - r26988377;
double r26988387 = r26988384 * r26988386;
double r26988388 = r26988377 + r26988387;
double r26988389 = 9.426461174744231e-94;
bool r26988390 = r26988374 <= r26988389;
double r26988391 = r26988377 / r26988379;
double r26988392 = r26988385 / r26988379;
double r26988393 = r26988391 - r26988392;
double r26988394 = r26988393 * r26988381;
double r26988395 = r26988394 + r26988385;
double r26988396 = r26988390 ? r26988395 : r26988388;
double r26988397 = r26988376 ? r26988388 : r26988396;
return r26988397;
}




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.9 |
| Herbie | 9.3 |
if a < -1.2559578325599803e-148 or 9.426461174744231e-94 < a Initial program 21.6
Simplified11.0
rmApplied fma-udef11.0
rmApplied div-inv11.1
Applied associate-*l*8.7
Simplified8.6
rmApplied clear-num8.7
if -1.2559578325599803e-148 < a < 9.426461174744231e-94Initial program 29.8
Simplified24.7
rmApplied fma-udef24.7
rmApplied div-inv24.8
Applied associate-*l*20.5
Simplified20.4
Taylor expanded around inf 14.0
Simplified10.8
Final simplification9.3
herbie shell --seed 2019162 +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))))