x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -5.26347676420508287 \cdot 10^{-166}:\\
\;\;\;\;\left(\left(y - x\right) \cdot \frac{1}{a - t}\right) \cdot \left(z - t\right) + x\\
\mathbf{elif}\;a \le 5.08907771502049331 \cdot 10^{-194}:\\
\;\;\;\;\mathsf{fma}\left(\frac{x}{t}, z, y - \frac{z \cdot y}{t}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y - x, \frac{z - t}{a - t}, x\right)\\
\end{array}double code(double x, double y, double z, double t, double a) {
return (x + (((y - x) * (z - t)) / (a - t)));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((a <= -5.263476764205083e-166)) {
VAR = ((((y - x) * (1.0 / (a - t))) * (z - t)) + x);
} else {
double VAR_1;
if ((a <= 5.0890777150204933e-194)) {
VAR_1 = fma((x / t), z, (y - ((z * y) / t)));
} else {
VAR_1 = fma((y - x), ((z - t) / (a - t)), x);
}
VAR = VAR_1;
}
return VAR;
}




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.3 |
| Herbie | 11.1 |
if a < -5.263476764205083e-166Initial program 23.7
Simplified12.0
rmApplied fma-udef12.1
rmApplied div-inv12.1
if -5.263476764205083e-166 < a < 5.0890777150204933e-194Initial program 29.6
Simplified25.4
rmApplied fma-udef25.5
rmApplied div-inv25.5
Applied associate-*l*20.4
Simplified20.3
Taylor expanded around inf 11.4
Simplified11.5
if 5.0890777150204933e-194 < a Initial program 23.2
Simplified12.7
rmApplied fma-udef12.8
rmApplied div-inv12.9
Applied associate-*l*10.0
Simplified9.9
rmApplied fma-def9.9
Final simplification11.1
herbie shell --seed 2020078 +o rules:numerics
(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))))