\frac{x - y \cdot z}{t - a \cdot z}\frac{x - y \cdot z}{t - z \cdot a}double f(double x, double y, double z, double t, double a) {
double r35198485 = x;
double r35198486 = y;
double r35198487 = z;
double r35198488 = r35198486 * r35198487;
double r35198489 = r35198485 - r35198488;
double r35198490 = t;
double r35198491 = a;
double r35198492 = r35198491 * r35198487;
double r35198493 = r35198490 - r35198492;
double r35198494 = r35198489 / r35198493;
return r35198494;
}
double f(double x, double y, double z, double t, double a) {
double r35198495 = x;
double r35198496 = y;
double r35198497 = z;
double r35198498 = r35198496 * r35198497;
double r35198499 = r35198495 - r35198498;
double r35198500 = t;
double r35198501 = a;
double r35198502 = r35198497 * r35198501;
double r35198503 = r35198500 - r35198502;
double r35198504 = r35198499 / r35198503;
return r35198504;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.7 |
|---|---|
| Target | 1.7 |
| Herbie | 10.7 |
Initial program 10.7
rmApplied clear-num11.1
rmApplied *-un-lft-identity11.1
Applied *-un-lft-identity11.1
Applied times-frac11.1
Applied add-cube-cbrt11.1
Applied times-frac11.1
Simplified11.1
Simplified10.7
Final simplification10.7
herbie shell --seed 2019168 +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.0 (- t (* a z)))) (- (/ x (- t (* a z))) (/ y (- (/ t z) a)))))
(/ (- x (* y z)) (- t (* a z))))