\frac{x - y \cdot z}{t - a \cdot z}\frac{x}{t - a \cdot z} - \frac{y}{\frac{t}{z} - a}double f(double x, double y, double z, double t, double a) {
double r12946560 = x;
double r12946561 = y;
double r12946562 = z;
double r12946563 = r12946561 * r12946562;
double r12946564 = r12946560 - r12946563;
double r12946565 = t;
double r12946566 = a;
double r12946567 = r12946566 * r12946562;
double r12946568 = r12946565 - r12946567;
double r12946569 = r12946564 / r12946568;
return r12946569;
}
double f(double x, double y, double z, double t, double a) {
double r12946570 = x;
double r12946571 = t;
double r12946572 = a;
double r12946573 = z;
double r12946574 = r12946572 * r12946573;
double r12946575 = r12946571 - r12946574;
double r12946576 = r12946570 / r12946575;
double r12946577 = y;
double r12946578 = r12946571 / r12946573;
double r12946579 = r12946578 - r12946572;
double r12946580 = r12946577 / r12946579;
double r12946581 = r12946576 - r12946580;
return r12946581;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.0 |
|---|---|
| Target | 1.7 |
| Herbie | 2.9 |
Initial program 10.0
rmApplied div-sub10.0
rmApplied associate-/l*7.7
Taylor expanded around 0 2.9
Final simplification2.9
herbie shell --seed 2019156 +o rules:numerics
(FPCore (x y z t a)
:name "Diagrams.Solve.Tridiagonal:solveTriDiagonal from diagrams-solve-0.1, A"
:herbie-target
(if (< z -32113435955957344.0) (- (/ 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))))