x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -3.85925241219433013420020968270754857507 \cdot 10^{-145}:\\
\;\;\;\;x + \left(y - x\right) \cdot \frac{z - t}{a - t}\\
\mathbf{elif}\;a \le 4.817774934231587927443345859010053721242 \cdot 10^{-234}:\\
\;\;\;\;\left(y + \frac{x \cdot z}{t}\right) - \frac{z \cdot y}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y - x}{\frac{a}{z - t} - \frac{t}{z - t}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r528319 = x;
double r528320 = y;
double r528321 = r528320 - r528319;
double r528322 = z;
double r528323 = t;
double r528324 = r528322 - r528323;
double r528325 = r528321 * r528324;
double r528326 = a;
double r528327 = r528326 - r528323;
double r528328 = r528325 / r528327;
double r528329 = r528319 + r528328;
return r528329;
}
double f(double x, double y, double z, double t, double a) {
double r528330 = a;
double r528331 = -3.85925241219433e-145;
bool r528332 = r528330 <= r528331;
double r528333 = x;
double r528334 = y;
double r528335 = r528334 - r528333;
double r528336 = z;
double r528337 = t;
double r528338 = r528336 - r528337;
double r528339 = r528330 - r528337;
double r528340 = r528338 / r528339;
double r528341 = r528335 * r528340;
double r528342 = r528333 + r528341;
double r528343 = 4.817774934231588e-234;
bool r528344 = r528330 <= r528343;
double r528345 = r528333 * r528336;
double r528346 = r528345 / r528337;
double r528347 = r528334 + r528346;
double r528348 = r528336 * r528334;
double r528349 = r528348 / r528337;
double r528350 = r528347 - r528349;
double r528351 = r528330 / r528338;
double r528352 = r528337 / r528338;
double r528353 = r528351 - r528352;
double r528354 = r528335 / r528353;
double r528355 = r528333 + r528354;
double r528356 = r528344 ? r528350 : r528355;
double r528357 = r528332 ? r528342 : r528356;
return r528357;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 25.2 |
|---|---|
| Target | 9.6 |
| Herbie | 10.5 |
if a < -3.85925241219433e-145Initial program 23.9
rmApplied *-un-lft-identity23.9
Applied times-frac9.4
Simplified9.4
if -3.85925241219433e-145 < a < 4.817774934231588e-234Initial program 31.1
Taylor expanded around inf 11.5
if 4.817774934231588e-234 < a Initial program 24.0
rmApplied associate-/l*11.1
rmApplied div-sub11.1
Final simplification10.5
herbie shell --seed 2019323
(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))))