x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -4.0098790375020589 \cdot 10^{-70}:\\
\;\;\;\;x + \frac{y - x}{a - t} \cdot \left(z - t\right)\\
\mathbf{elif}\;a \le 3.8034022935305039 \cdot 10^{-125}:\\
\;\;\;\;\left(y + \frac{x \cdot z}{t}\right) - \frac{z \cdot y}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y - x}{\frac{a - t}{z - t}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r1227058 = x;
double r1227059 = y;
double r1227060 = r1227059 - r1227058;
double r1227061 = z;
double r1227062 = t;
double r1227063 = r1227061 - r1227062;
double r1227064 = r1227060 * r1227063;
double r1227065 = a;
double r1227066 = r1227065 - r1227062;
double r1227067 = r1227064 / r1227066;
double r1227068 = r1227058 + r1227067;
return r1227068;
}
double f(double x, double y, double z, double t, double a) {
double r1227069 = a;
double r1227070 = -4.009879037502059e-70;
bool r1227071 = r1227069 <= r1227070;
double r1227072 = x;
double r1227073 = y;
double r1227074 = r1227073 - r1227072;
double r1227075 = t;
double r1227076 = r1227069 - r1227075;
double r1227077 = r1227074 / r1227076;
double r1227078 = z;
double r1227079 = r1227078 - r1227075;
double r1227080 = r1227077 * r1227079;
double r1227081 = r1227072 + r1227080;
double r1227082 = 3.803402293530504e-125;
bool r1227083 = r1227069 <= r1227082;
double r1227084 = r1227072 * r1227078;
double r1227085 = r1227084 / r1227075;
double r1227086 = r1227073 + r1227085;
double r1227087 = r1227078 * r1227073;
double r1227088 = r1227087 / r1227075;
double r1227089 = r1227086 - r1227088;
double r1227090 = r1227076 / r1227079;
double r1227091 = r1227074 / r1227090;
double r1227092 = r1227072 + r1227091;
double r1227093 = r1227083 ? r1227089 : r1227092;
double r1227094 = r1227071 ? r1227081 : r1227093;
return r1227094;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.6 |
|---|---|
| Target | 9.1 |
| Herbie | 11.6 |
if a < -4.009879037502059e-70Initial program 23.0
rmApplied associate-/l*8.1
rmApplied associate-/r/10.2
if -4.009879037502059e-70 < a < 3.803402293530504e-125Initial program 29.2
Taylor expanded around inf 16.9
if 3.803402293530504e-125 < a Initial program 22.2
rmApplied associate-/l*8.2
Final simplification11.6
herbie shell --seed 2020047
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:linMap from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< a -1.6153062845442575e-142) (+ x (* (/ (- y x) 1) (/ (- z t) (- a t)))) (if (< a 3.774403170083174e-182) (- y (* (/ z t) (- y x))) (+ x (* (/ (- y x) 1) (/ (- z t) (- a t))))))
(+ x (/ (* (- y x) (- z t)) (- a t))))