\frac{x - y \cdot z}{t - a \cdot z}\frac{x - y \cdot z}{t - a \cdot z}double f(double x, double y, double z, double t, double a) {
double r708772 = x;
double r708773 = y;
double r708774 = z;
double r708775 = r708773 * r708774;
double r708776 = r708772 - r708775;
double r708777 = t;
double r708778 = a;
double r708779 = r708778 * r708774;
double r708780 = r708777 - r708779;
double r708781 = r708776 / r708780;
return r708781;
}
double f(double x, double y, double z, double t, double a) {
double r708782 = x;
double r708783 = y;
double r708784 = z;
double r708785 = r708783 * r708784;
double r708786 = r708782 - r708785;
double r708787 = t;
double r708788 = a;
double r708789 = r708788 * r708784;
double r708790 = r708787 - r708789;
double r708791 = r708786 / r708790;
return r708791;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.5 |
|---|---|
| Target | 1.9 |
| Herbie | 10.5 |
Initial program 10.5
rmApplied clear-num10.8
rmApplied *-un-lft-identity10.8
Applied *-un-lft-identity10.8
Applied times-frac10.8
Applied add-cube-cbrt10.8
Applied times-frac10.8
Simplified10.8
Simplified10.5
Final simplification10.5
herbie shell --seed 2020018
(FPCore (x y z t a)
:name "Diagrams.Solve.Tridiagonal:solveTriDiagonal from diagrams-solve-0.1, A"
:precision binary64
:herbie-target
(if (< z -32113435955957344) (- (/ x (- t (* a z))) (/ y (- (/ t z) a))) (if (< z 3.5139522372978296e-86) (* (- x (* y z)) (/ 1 (- t (* a z)))) (- (/ x (- t (* a z))) (/ y (- (/ t z) a)))))
(/ (- x (* y z)) (- t (* a z))))