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 r31264033 = x;
double r31264034 = y;
double r31264035 = r31264034 - r31264033;
double r31264036 = z;
double r31264037 = t;
double r31264038 = r31264036 - r31264037;
double r31264039 = r31264035 * r31264038;
double r31264040 = a;
double r31264041 = r31264040 - r31264037;
double r31264042 = r31264039 / r31264041;
double r31264043 = r31264033 + r31264042;
return r31264043;
}
double f(double x, double y, double z, double t, double a) {
double r31264044 = x;
double r31264045 = y;
double r31264046 = r31264045 - r31264044;
double r31264047 = z;
double r31264048 = t;
double r31264049 = r31264047 - r31264048;
double r31264050 = r31264046 * r31264049;
double r31264051 = a;
double r31264052 = r31264051 - r31264048;
double r31264053 = r31264050 / r31264052;
double r31264054 = r31264044 + r31264053;
double r31264055 = -1.4035645823185852e-264;
bool r31264056 = r31264054 <= r31264055;
double r31264057 = r31264049 / r31264052;
double r31264058 = r31264046 * r31264057;
double r31264059 = r31264058 + r31264044;
double r31264060 = 0.0;
bool r31264061 = r31264054 <= r31264060;
double r31264062 = r31264047 * r31264044;
double r31264063 = r31264062 / r31264048;
double r31264064 = r31264063 + r31264045;
double r31264065 = r31264047 * r31264045;
double r31264066 = r31264065 / r31264048;
double r31264067 = r31264064 - r31264066;
double r31264068 = r31264061 ? r31264067 : r31264059;
double r31264069 = r31264056 ? r31264059 : r31264068;
return r31264069;
}




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))))