\frac{x + \frac{y \cdot z}{t}}{\left(a + 1\right) + \frac{y \cdot b}{t}}\frac{z}{b + \frac{t}{y} \cdot \left(a + 1\right)} + \frac{x}{a + \left(1 + \frac{1}{\frac{t}{b \cdot y}}\right)}(FPCore (x y z t a b) :precision binary64 (/ (+ x (/ (* y z) t)) (+ (+ a 1.0) (/ (* y b) t))))
(FPCore (x y z t a b) :precision binary64 (+ (/ z (+ b (* (/ t y) (+ a 1.0)))) (/ x (+ a (+ 1.0 (/ 1.0 (/ t (* b y))))))))
double code(double x, double y, double z, double t, double a, double b) {
return (x + ((y * z) / t)) / ((a + 1.0) + ((y * b) / t));
}
double code(double x, double y, double z, double t, double a, double b) {
return (z / (b + ((t / y) * (a + 1.0)))) + (x / (a + (1.0 + (1.0 / (t / (b * y))))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 16.5 |
|---|---|
| Target | 13.0 |
| Herbie | 4.5 |
Initial program 16.5
Taylor expanded around 0 15.2
Simplified15.2
rmApplied associate-/l*_binary6413.2
Simplified14.3
Taylor expanded around 0 8.0
Simplified4.5
rmApplied clear-num_binary644.5
Final simplification4.5
herbie shell --seed 2021110
(FPCore (x y z t a b)
:name "Diagrams.Solve.Tridiagonal:solveCyclicTriDiagonal from diagrams-solve-0.1, B"
:precision binary64
:herbie-target
(if (< t -1.3659085366310088e-271) (* 1.0 (* (+ x (* (/ y t) z)) (/ 1.0 (+ (+ a 1.0) (* (/ y t) b))))) (if (< t 3.036967103737246e-130) (/ z b) (* 1.0 (* (+ x (* (/ y t) z)) (/ 1.0 (+ (+ a 1.0) (* (/ y t) b)))))))
(/ (+ x (/ (* y z) t)) (+ (+ a 1.0) (/ (* y b) t))))