



Bits error versus a




Bits error versus b




Bits error versus eps
Results
| Original | 58.6 |
|---|---|
| Target | 14.4 |
| Herbie | 3.0 |
if eps < -1.0077585081869481e-76Initial program 50.9
Simplified7.5
rmApplied add-cube-cbrt8.1
Applied add-cube-cbrt8.0
Applied times-frac8.0
Applied associate-*l*7.7
if -1.0077585081869481e-76 < eps Initial program 59.5
Simplified29.9
Taylor expanded around 0 2.5
Final simplification3.0
herbie shell --seed 2019007 +o rules:numerics
(FPCore (a b eps)
:name "expq3 (problem 3.4.2)"
:pre (and (< -1 eps) (< eps 1))
:herbie-target
(/ (+ a b) (* a b))
(/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1))))
Time bar (total: 28.0s)Debug log
1 calls. Slowest were:
| 446.0ms | (/ (* eps (- (exp (* (+ a b) eps)) 1)) (* (- (exp (* a eps)) 1) (- (exp (* b eps)) 1))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 27.3b
Found 4 expressions with local error:
| 17.1b | (/ eps (expm1 (* eps a))) |
| 10.3b | (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) |
| 4.8b | (* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ eps (expm1 (* eps a)))) |
| 0.1b | (expm1 (* (+ a b) eps)) |
4 calls. Slowest were:
| 6.0ms | (* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ eps (expm1 (* eps a)))) |
| 3.0ms | (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) |
| 1.0ms | (/ eps (expm1 (* eps a))) |
4 calls. Slowest were:
| 407.0ms | (* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ eps (expm1 (* eps a)))) |
| 88.0ms | (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) |
| 55.0ms | (expm1 (* (+ a b) eps)) |
| 37.0ms | (/ eps (expm1 (* eps a))) |
67 calls. Slowest were:
| 348.0ms | (* (* (* (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b))) (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))) (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))) (* (* (/ eps (expm1 (* eps a))) (/ eps (expm1 (* eps a)))) (/ eps (expm1 (* eps a))))) |
| 300.0ms | (+ (- (log (expm1 (* (+ a b) eps))) (log (expm1 (* eps b)))) (log (/ eps (expm1 (* eps a))))) |
| 190.0ms | (+ (log (/ (expm1 (* (+ a b) eps)) (expm1 (* eps b)))) (log (/ eps (expm1 (* eps a))))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.0b
Found 1 expressions with local error:
| 0.0b | (+ (/ 1 a) (/ 1 b)) |
1 calls. Slowest were:
| 7.0ms | (+ (/ 1 a) (/ 1 b)) |
1 calls. Slowest were:
| 12.0ms | (+ (/ 1 a) (/ 1 b)) |
12 calls. Slowest were:
| 6.0ms | (* (exp (/ 1 a)) (exp (/ 1 b))) |
| 4.0ms | (+ (* 1 b) (* a 1)) |
| 4.0ms | (* (cbrt (/ 1 a)) (cbrt (/ 1 a))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 17.1b | (/ eps (expm1 (* eps a))) |
| 15.0b | (cbrt (expm1 (* eps b))) |
| 15.0b | (cbrt (expm1 (* eps b))) |
| 15.0b | (cbrt (expm1 (* eps b))) |
4 calls. Slowest were:
| 1.0ms | (cbrt (expm1 (* eps b))) |
| 1.0ms | (cbrt (expm1 (* eps b))) |
| 1.0ms | (cbrt (expm1 (* eps b))) |
4 calls. Slowest were:
| 111.0ms | (cbrt (expm1 (* eps b))) |
| 99.0ms | (cbrt (expm1 (* eps b))) |
| 92.0ms | (cbrt (expm1 (* eps b))) |
| 38.0ms | (/ eps (expm1 (* eps a))) |
24 calls. Slowest were:
| 227.0ms | (+ (exp (* 1/3 (+ (log eps) (log b)))) (+ (* 1/36 (* (exp (* 1/3 (+ (log eps) (log b)))) (* (pow eps 2) (pow b 2)))) (* 1/6 (* (exp (* 1/3 (+ (log eps) (log b)))) (* eps b))))) |
| 199.0ms | (+ (exp (* 1/3 (+ (log eps) (log b)))) (+ (* 1/36 (* (exp (* 1/3 (+ (log eps) (log b)))) (* (pow eps 2) (pow b 2)))) (* 1/6 (* (exp (* 1/3 (+ (log eps) (log b)))) (* eps b))))) |
| 196.0ms | (+ (exp (* 1/3 (+ (log eps) (log b)))) (+ (* 1/36 (* (exp (* 1/3 (+ (log eps) (log b)))) (* (pow eps 2) (pow b 2)))) (* 1/6 (* (exp (* 1/3 (+ (log eps) (log b)))) (* eps b))))) |
4 alts after pruning (2 fresh and 2 done)
Merged error: 0.0b
Found 3 expressions with local error:
| 0.1b | (fma (sqrt (/ 1 a)) (sqrt (/ 1 a)) (/ 1 b)) |
| 0.0b | (sqrt (/ 1 a)) |
| 0.0b | (sqrt (/ 1 a)) |
3 calls. Slowest were:
| 1.0ms | (sqrt (/ 1 a)) |
| 1.0ms | (sqrt (/ 1 a)) |
| 0.0ms | (fma (sqrt (/ 1 a)) (sqrt (/ 1 a)) (/ 1 b)) |
3 calls. Slowest were:
| 44.0ms | (fma (sqrt (/ 1 a)) (sqrt (/ 1 a)) (/ 1 b)) |
| 11.0ms | (sqrt (/ 1 a)) |
| 10.0ms | (sqrt (/ 1 a)) |
21 calls. Slowest were:
| 2.0ms | (sqrt (* (cbrt (/ 1 a)) (cbrt (/ 1 a)))) |
| 2.0ms | (sqrt (* (cbrt (/ 1 a)) (cbrt (/ 1 a)))) |
| 2.0ms | (sqrt (sqrt (/ 1 a))) |
4 alts after pruning (1 fresh and 3 done)
Merged error: 0.0b
9.3% (3.0b remaining)
Error of 3.0b against oracle of 0.0b and baseline of 3.3b