\frac{x}{y - z \cdot t}\frac{x}{y - z \cdot t}double f(double x, double y, double z, double t) {
double r525948 = x;
double r525949 = y;
double r525950 = z;
double r525951 = t;
double r525952 = r525950 * r525951;
double r525953 = r525949 - r525952;
double r525954 = r525948 / r525953;
return r525954;
}
double f(double x, double y, double z, double t) {
double r525955 = x;
double r525956 = y;
double r525957 = z;
double r525958 = t;
double r525959 = r525957 * r525958;
double r525960 = r525956 - r525959;
double r525961 = r525955 / r525960;
return r525961;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 3.0 |
|---|---|
| Target | 1.7 |
| Herbie | 3.0 |
Initial program 3.0
Final simplification3.0
herbie shell --seed 2019303 +o rules:numerics
(FPCore (x y z t)
:name "Diagrams.Solve.Tridiagonal:solveTriDiagonal from diagrams-solve-0.1, B"
:precision binary64
:herbie-target
(if (< x -1.618195973607049e50) (/ 1 (- (/ y x) (* (/ z x) t))) (if (< x 2.13783064348764444e131) (/ x (- y (* z t))) (/ 1 (- (/ y x) (* (/ z x) t)))))
(/ x (- y (* z t))))