\frac{x - y \cdot z}{t - a \cdot z}x \cdot \frac{1}{t - a \cdot z} - \frac{y}{\frac{t}{z} - a}double f(double x, double y, double z, double t, double a) {
double r477764 = x;
double r477765 = y;
double r477766 = z;
double r477767 = r477765 * r477766;
double r477768 = r477764 - r477767;
double r477769 = t;
double r477770 = a;
double r477771 = r477770 * r477766;
double r477772 = r477769 - r477771;
double r477773 = r477768 / r477772;
return r477773;
}
double f(double x, double y, double z, double t, double a) {
double r477774 = x;
double r477775 = 1.0;
double r477776 = t;
double r477777 = a;
double r477778 = z;
double r477779 = r477777 * r477778;
double r477780 = r477776 - r477779;
double r477781 = r477775 / r477780;
double r477782 = r477774 * r477781;
double r477783 = y;
double r477784 = r477776 / r477778;
double r477785 = r477784 - r477777;
double r477786 = r477783 / r477785;
double r477787 = r477782 - r477786;
return r477787;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.4 |
|---|---|
| Target | 1.6 |
| Herbie | 2.9 |
if z < -3.528978657649998e+26 or 4.507429170566249e-285 < z Initial program 14.4
rmApplied div-sub14.4
rmApplied associate-/l*9.7
rmApplied div-sub9.7
Simplified2.8
rmApplied div-inv2.9
if -3.528978657649998e+26 < z < 4.507429170566249e-285Initial program 0.4
rmApplied div-sub0.4
rmApplied clear-num0.7
Final simplification2.9
herbie shell --seed 2019298
(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.51395223729782958e-86) (* (- x (* y z)) (/ 1 (- t (* a z)))) (- (/ x (- t (* a z))) (/ y (- (/ t z) a)))))
(/ (- x (* y z)) (- t (* a z))))