\frac{x}{y - z \cdot t}\frac{x}{y - z \cdot t}double f(double x, double y, double z, double t) {
double r1257760 = x;
double r1257761 = y;
double r1257762 = z;
double r1257763 = t;
double r1257764 = r1257762 * r1257763;
double r1257765 = r1257761 - r1257764;
double r1257766 = r1257760 / r1257765;
return r1257766;
}
double f(double x, double y, double z, double t) {
double r1257767 = x;
double r1257768 = y;
double r1257769 = z;
double r1257770 = t;
double r1257771 = r1257769 * r1257770;
double r1257772 = r1257768 - r1257771;
double r1257773 = r1257767 / r1257772;
return r1257773;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.9 |
|---|---|
| Target | 1.8 |
| Herbie | 2.9 |
Initial program 2.9
Final simplification2.9
herbie shell --seed 2019318
(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))))