



Bits error versus i




Bits error versus n
Results
| Original | 42.1 |
|---|---|
| Target | 41.9 |
| Herbie | 19.3 |
if n < -1.5274599045984505e-107Initial program 39.4
rmApplied pow-to-exp45.4
Applied expm1-def41.1
Simplified26.1
rmApplied div-inv26.1
Applied *-un-lft-identity26.1
Applied times-frac25.4
Applied associate-*r*25.4
Simplified25.4
if -1.5274599045984505e-107 < n < 2.195458028809774e-180Initial program 22.8
rmApplied pow-to-exp22.8
Applied expm1-def13.4
Simplified17.5
Taylor expanded around 0 14.4
if 2.195458028809774e-180 < n Initial program 57.5
rmApplied pow-to-exp57.5
Applied expm1-def51.1
Simplified15.6
rmApplied *-commutative15.6
Final simplification19.3
herbie shell --seed 2019004 +o rules:numerics
(FPCore (i n)
:name "Compound Interest"
:herbie-target
(* 100 (/ (- (exp (* n (if (== (+ 1 (/ i n)) 1) (/ i n) (/ (* (/ i n) (log (+ 1 (/ i n)))) (- (+ (/ i n) 1) 1))))) 1) (/ i n)))
(* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))))
Time bar (total: 28.8s)Debug log
1 calls. Slowest were:
| 264.0ms | (* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 41.4b
Found 4 expressions with local error:
| 13.3b | (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n)) |
| 7.5b | (pow (+ 1 (/ i n)) n) |
| 1.8b | (- (pow (+ 1 (/ i n)) n) 1) |
| 0.0b | (* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))) |
4 calls. Slowest were:
| 28.0ms | (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n)) |
| 20.0ms | (- (pow (+ 1 (/ i n)) n) 1) |
| 18.0ms | (* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))) |
4 calls. Slowest were:
| 220.0ms | (* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))) |
| 215.0ms | (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n)) |
| 169.0ms | (- (pow (+ 1 (/ i n)) n) 1) |
| 108.0ms | (pow (+ 1 (/ i n)) n) |
105 calls. Slowest were:
| 317.0ms | (/ (* (- (exp (* (- (log (/ 1 n)) (log (/ 1 i))) n)) 1) n) i) |
| 303.0ms | (/ (* (* (- (pow (+ 1 (/ i n)) n) 1) (- (pow (+ 1 (/ i n)) n) 1)) (- (pow (+ 1 (/ i n)) n) 1)) (* (* (/ i n) (/ i n)) (/ i n))) |
| 223.0ms | (* 100 (/ (* n (- (exp (* n (- (log (/ -1 n)) (log (/ -1 i))))) 1)) i)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 12.7b
Found 4 expressions with local error:
| 13.1b | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 13.1b | (* n (log1p (/ i n))) |
| 6.6b | (log1p (/ i n)) |
| 0.1b | (* 100 (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
4 calls. Slowest were:
| 3.0ms | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 1.0ms | (* 100 (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 1.0ms | (* n (log1p (/ i n))) |
4 calls. Slowest were:
| 309.0ms | (* 100 (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 279.0ms | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 39.0ms | (log1p (/ i n)) |
| 36.0ms | (* n (log1p (/ i n))) |
66 calls. Slowest were:
| 189.0ms | (/ (* (* (expm1 (* n (log1p (/ i n)))) (expm1 (* n (log1p (/ i n))))) (expm1 (* n (log1p (/ i n))))) (* (* (/ i n) (/ i n)) (/ i n))) |
| 81.0ms | (- (log (/ -1 n)) (log (/ -1 i))) |
| 57.0ms | (/ (* (- (exp i) 1) n) i) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 12.6b
Found 4 expressions with local error:
| 13.1b | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 13.1b | (* n (log1p (/ i n))) |
| 6.6b | (log1p (/ i n)) |
| 0.1b | (* (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) 100) |
4 calls. Slowest were:
| 6.0ms | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 6.0ms | (* (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) 100) |
| 2.0ms | (* n (log1p (/ i n))) |
4 calls. Slowest were:
| 304.0ms | (* (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) 100) |
| 285.0ms | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 48.0ms | (* n (log1p (/ i n))) |
| 35.0ms | (log1p (/ i n)) |
66 calls. Slowest were:
| 181.0ms | (/ (* (* (expm1 (* n (log1p (/ i n)))) (expm1 (* n (log1p (/ i n))))) (expm1 (* n (log1p (/ i n))))) (* (* (/ i n) (/ i n)) (/ i n))) |
| 58.0ms | (/ (* (- (* 100 (exp i)) 100) n) i) |
| 56.0ms | (/ (* (- (exp i) 1) n) i) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 12.6b
Found 4 expressions with local error:
| 13.1b | (* (log1p (/ i n)) n) |
| 13.1b | (* (log1p (/ i n)) n) |
| 13.1b | (* n (log1p (/ i n))) |
| 10.2b | (* (* (cbrt (expm1 (* (log1p (/ i n)) n))) (cbrt (expm1 (* (log1p (/ i n)) n)))) (/ (cbrt (expm1 (* n (log1p (/ i n))))) (/ i n))) |
4 calls. Slowest were:
| 13.0ms | (* (* (cbrt (expm1 (* (log1p (/ i n)) n))) (cbrt (expm1 (* (log1p (/ i n)) n)))) (/ (cbrt (expm1 (* n (log1p (/ i n))))) (/ i n))) |
| 2.0ms | (* (log1p (/ i n)) n) |
| 2.0ms | (* (log1p (/ i n)) n) |
4 calls. Slowest were:
| 218.0ms | (* (* (cbrt (expm1 (* (log1p (/ i n)) n))) (cbrt (expm1 (* (log1p (/ i n)) n)))) (/ (cbrt (expm1 (* n (log1p (/ i n))))) (/ i n))) |
| 68.0ms | (* (log1p (/ i n)) n) |
| 55.0ms | (* (log1p (/ i n)) n) |
| 53.0ms | (* n (log1p (/ i n))) |
59 calls. Slowest were:
| 676.0ms | (+ (log (* (cbrt (expm1 (* (log1p (/ i n)) n))) (cbrt (expm1 (* (log1p (/ i n)) n))))) (log (/ (cbrt (expm1 (* n (log1p (/ i n))))) (/ i n)))) |
| 511.0ms | (+ (log (* (cbrt (expm1 (* (log1p (/ i n)) n))) (cbrt (expm1 (* (log1p (/ i n)) n))))) (- (log (cbrt (expm1 (* n (log1p (/ i n)))))) (log (/ i n)))) |
| 491.0ms | (* (* (expm1 (* (log1p (/ i n)) n)) (expm1 (* (log1p (/ i n)) n))) (* (* (/ (cbrt (expm1 (* n (log1p (/ i n))))) (/ i n)) (/ (cbrt (expm1 (* n (log1p (/ i n))))) (/ i n))) (/ (cbrt (expm1 (* n (log1p (/ i n))))) (/ i n)))) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 12.6b
9.3% (9.7b remaining)
Error of 19.3b against oracle of 9.5b and baseline of 20.3b