\frac{x + \frac{y \cdot z - x}{t \cdot z - x}}{x + 1}\begin{array}{l}
\mathbf{if}\;z \le 6.8077680963010561 \cdot 10^{135}:\\
\;\;\;\;\frac{\left(\frac{z}{t \cdot z - x} \cdot y + x\right) - \frac{x}{t \cdot z - x}}{x + 1}\\
\mathbf{else}:\\
\;\;\;\;\frac{x + \frac{y}{t}}{x + 1}\\
\end{array}double code(double x, double y, double z, double t) {
return ((x + (((y * z) - x) / ((t * z) - x))) / (x + 1.0));
}
double code(double x, double y, double z, double t) {
double VAR;
if ((z <= 6.807768096301056e+135)) {
VAR = (((((z / ((t * z) - x)) * y) + x) - (x / ((t * z) - x))) / (x + 1.0));
} else {
VAR = ((x + (y / t)) / (x + 1.0));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.3 |
|---|---|
| Target | 0.4 |
| Herbie | 2.2 |
if z < 6.807768096301056e+135Initial program 4.9
rmApplied div-sub4.9
Applied associate-+r-4.9
Simplified4.0
rmApplied clear-num4.1
rmApplied fma-udef4.1
Simplified1.5
if 6.807768096301056e+135 < z Initial program 23.2
Taylor expanded around inf 6.9
Final simplification2.2
herbie shell --seed 2020105 +o rules:numerics
(FPCore (x y z t)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, A"
:precision binary64
:herbie-target
(/ (+ x (- (/ y (- t (/ x z))) (/ x (- (* t z) x)))) (+ x 1))
(/ (+ x (/ (- (* y z) x) (- (* t z) x))) (+ x 1)))