x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}\begin{array}{l}
\mathbf{if}\;x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t} \le -1.403564582318585211719020373880695846914 \cdot 10^{-264}:\\
\;\;\;\;\left(y - x\right) \cdot \frac{z - t}{a - t} + x\\
\mathbf{elif}\;x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t} \le 0.0:\\
\;\;\;\;\left(\frac{z \cdot x}{t} + y\right) - \frac{z \cdot y}{t}\\
\mathbf{else}:\\
\;\;\;\;\left(y - x\right) \cdot \frac{z - t}{a - t} + x\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r33053321 = x;
double r33053322 = y;
double r33053323 = r33053322 - r33053321;
double r33053324 = z;
double r33053325 = t;
double r33053326 = r33053324 - r33053325;
double r33053327 = r33053323 * r33053326;
double r33053328 = a;
double r33053329 = r33053328 - r33053325;
double r33053330 = r33053327 / r33053329;
double r33053331 = r33053321 + r33053330;
return r33053331;
}
double f(double x, double y, double z, double t, double a) {
double r33053332 = x;
double r33053333 = y;
double r33053334 = r33053333 - r33053332;
double r33053335 = z;
double r33053336 = t;
double r33053337 = r33053335 - r33053336;
double r33053338 = r33053334 * r33053337;
double r33053339 = a;
double r33053340 = r33053339 - r33053336;
double r33053341 = r33053338 / r33053340;
double r33053342 = r33053332 + r33053341;
double r33053343 = -1.4035645823185852e-264;
bool r33053344 = r33053342 <= r33053343;
double r33053345 = r33053337 / r33053340;
double r33053346 = r33053334 * r33053345;
double r33053347 = r33053346 + r33053332;
double r33053348 = 0.0;
bool r33053349 = r33053342 <= r33053348;
double r33053350 = r33053335 * r33053332;
double r33053351 = r33053350 / r33053336;
double r33053352 = r33053351 + r33053333;
double r33053353 = r33053335 * r33053333;
double r33053354 = r33053353 / r33053336;
double r33053355 = r33053352 - r33053354;
double r33053356 = r33053349 ? r33053355 : r33053347;
double r33053357 = r33053344 ? r33053347 : r33053356;
return r33053357;
}




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 | 9.5 |
| Herbie | 8.8 |
if (+ x (/ (* (- y x) (- z t)) (- a t))) < -1.4035645823185852e-264 or 0.0 < (+ x (/ (* (- y x) (- z t)) (- a t))) Initial program 21.1
rmApplied associate-/l*7.7
rmApplied div-inv7.9
Simplified7.8
if -1.4035645823185852e-264 < (+ x (/ (* (- y x) (- z t)) (- a t))) < 0.0Initial program 57.6
Taylor expanded around inf 19.6
Final simplification8.8
herbie shell --seed 2019172
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:linMap from Chart-1.5.3"
:herbie-target
(if (< a -1.6153062845442575e-142) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t)))) (if (< a 3.774403170083174e-182) (- y (* (/ z t) (- y x))) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t))))))
(+ x (/ (* (- y x) (- z t)) (- a t))))