



Bits error versus x




Bits error versus y
Results
| Original | 20.2 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
Initial program 20.2
Initial simplification20.2
rmApplied add-sqr-sqrt20.2
Applied times-frac20.2
Simplified20.2
Simplified0.0
rmApplied log1p-expm1-u0.0
Final simplification0.0
herbie shell --seed 2018362 +o rules:numerics
(FPCore (x y)
:name "Kahan p9 Example"
:pre (and (< 0 x 1) (< y 1))
:herbie-target
(if (< 0.5 (fabs (/ x y)) 2) (/ (* (- x y) (+ x y)) (+ (* x x) (* y y))) (- 1 (/ 2 (+ 1 (* (/ x y) (/ x y))))))
(/ (* (- x y) (+ x y)) (+ (* x x) (* y y))))
Time bar (total: 18.9s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 20.5b
Found 1 expressions with local error:
| 20.5b | (/ (* (- x y) (+ y x)) (fma x x (* y y))) |
1 calls. Slowest were:
| 13.0ms | (/ (* (- x y) (+ y x)) (fma x x (* y y))) |
1 calls. Slowest were:
| 59.0ms | (/ (* (- x y) (+ y x)) (fma x x (* y y))) |
25 calls. Slowest were:
| 713.0ms | (* (- (* x x) (* y y)) (- (* y y) (* x x))) |
| 418.0ms | (* (- (* x x) (* y y)) (+ y x)) |
| 390.0ms | (* (- (* x x) (* y y)) (+ (pow y 3) (pow x 3))) |
3 alts after pruning (3 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.0b | (/ (+ x y) (hypot x y)) |
| 0.0b | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 0.0b | (/ (- x y) (hypot x y)) |
| 0.0b | (hypot x y) |
4 calls. Slowest were:
| 16.0ms | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 6.0ms | (/ (+ x y) (hypot x y)) |
| 3.0ms | (/ (- x y) (hypot x y)) |
4 calls. Slowest were:
| 77.0ms | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 63.0ms | (/ (- x y) (hypot x y)) |
| 42.0ms | (/ (+ x y) (hypot x y)) |
| 28.0ms | (hypot x y) |
105 calls. Slowest were:
| 336.0ms | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 322.0ms | (+ (- (log (- x y)) (log (hypot x y))) (- (log (+ x y)) (log (hypot x y)))) |
| 301.0ms | (+ (- (log (- x y)) (log (hypot x y))) (log (/ (+ x y) (hypot x y)))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.3b | (expm1 (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y)))) |
| 0.0b | (/ (+ x y) (hypot x y)) |
| 0.0b | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 0.0b | (/ (- x y) (hypot x y)) |
4 calls. Slowest were:
| 9.0ms | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 3.0ms | (/ (+ x y) (hypot x y)) |
| 3.0ms | (/ (- x y) (hypot x y)) |
4 calls. Slowest were:
| 94.0ms | (expm1 (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y)))) |
| 52.0ms | (/ (- x y) (hypot x y)) |
| 50.0ms | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 49.0ms | (/ (+ x y) (hypot x y)) |
106 calls. Slowest were:
| 396.0ms | (+ (- (log (- x y)) (log (hypot x y))) (- (log (+ x y)) (log (hypot x y)))) |
| 341.0ms | (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) |
| 341.0ms | (+ (- (log (- x y)) (log (hypot x y))) (log (/ (+ x y) (hypot x y)))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.3b | (expm1 (cbrt (* (* (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y)))) (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y)))))) |
| 0.0b | (/ (+ x y) (hypot x y)) |
| 0.0b | (/ (+ x y) (hypot x y)) |
| 0.0b | (/ (+ x y) (hypot x y)) |
4 calls. Slowest were:
| 6.0ms | (/ (+ x y) (hypot x y)) |
| 6.0ms | (/ (+ x y) (hypot x y)) |
| 6.0ms | (/ (+ x y) (hypot x y)) |
4 calls. Slowest were:
| 145.0ms | (expm1 (cbrt (* (* (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y)))) (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y)))))) |
| 73.0ms | (/ (+ x y) (hypot x y)) |
| 42.0ms | (/ (+ x y) (hypot x y)) |
| 42.0ms | (/ (+ x y) (hypot x y)) |
94 calls. Slowest were:
| 1.2s | (cbrt (* (* (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))) (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y)))) (* (/ (- x y) (hypot x y)) (/ (+ x y) (hypot x y))))) |
| 314.0ms | (/ (* (* (+ x y) (+ x y)) (+ x y)) (* (* (hypot x y) (hypot x y)) (hypot x y))) |
| 275.0ms | (/ (* (* (+ x y) (+ x y)) (+ x y)) (* (* (hypot x y) (hypot x y)) (hypot x y))) |
4 alts after pruning (2 fresh and 2 done)
Merged error: 0b
0% (0.0b remaining)
Error of 0.0b against oracle of 0.0b and baseline of 0.0b