



Bits error versus i




Bits error versus n
Results
| Original | 42.0 |
|---|---|
| Target | 41.8 |
| Herbie | 20.0 |
if n < -8.764753048170504e+49Initial program 45.8
rmApplied add-exp-log53.6
Applied pow-exp53.6
Applied expm1-def53.6
Simplified27.4
Taylor expanded around inf 42.5
Simplified27.0
if -8.764753048170504e+49 < n < -3.869064973251871e-99Initial program 25.7
rmApplied associate-*r/25.6
if -3.869064973251871e-99 < n < -2.549728235978189e-174Initial program 18.1
rmApplied add-exp-log18.1
Applied pow-exp18.1
Applied expm1-def2.0
Simplified23.3
rmApplied add-sqr-sqrt25.4
Applied *-un-lft-identity25.4
Applied times-frac25.4
if -2.549728235978189e-174 < n < 5.831751074444269e-195Initial program 22.6
rmApplied add-exp-log22.6
Applied pow-exp22.6
Applied expm1-def15.2
Simplified15.0
Taylor expanded around 0 12.1
if 5.831751074444269e-195 < n Initial program 57.2
rmApplied add-exp-log57.2
Applied pow-exp57.2
Applied expm1-def50.7
Simplified15.5
rmApplied add-sqr-sqrt15.8
Applied associate-*r*15.8
Final simplification20.0
herbie shell --seed 2018362 +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: 23.0s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 39.9b
Found 4 expressions with local error:
| 10.2b | (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n)) |
| 9.2b | (pow (+ 1 (/ i n)) n) |
| 1.2b | (- (pow (+ 1 (/ i n)) n) 1) |
| 0.1b | (* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))) |
4 calls. Slowest were:
| 37.0ms | (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n)) |
| 25.0ms | (- (pow (+ 1 (/ i n)) n) 1) |
| 19.0ms | (* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))) |
4 calls. Slowest were:
| 253.0ms | (* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))) |
| 224.0ms | (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n)) |
| 168.0ms | (- (pow (+ 1 (/ i n)) n) 1) |
| 123.0ms | (pow (+ 1 (/ i n)) n) |
105 calls. Slowest were:
| 420.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))) |
| 326.0ms | (* 100 (/ (* n (- (exp (* n (- (log (/ -1 n)) (log (/ -1 i))))) 1)) i)) |
| 291.0ms | (/ (* (- (exp (* (- (log (/ 1 n)) (log (/ 1 i))) n)) 1) n) i) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 9.7b
Found 4 expressions with local error:
| 10.1b | (* n (log1p (/ i n))) |
| 10.1b | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 9.0b | (log1p (/ i n)) |
| 0.1b | (* 100 (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
4 calls. Slowest were:
| 6.0ms | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 3.0ms | (* 100 (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 2.0ms | (* n (log1p (/ i n))) |
4 calls. Slowest were:
| 311.0ms | (/ (expm1 (* n (log1p (/ i n)))) (/ i n)) |
| 288.0ms | (* 100 (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 63.0ms | (* n (log1p (/ i n))) |
| 42.0ms | (log1p (/ i n)) |
66 calls. Slowest were:
| 235.0ms | (/ (* (* (expm1 (* n (log1p (/ i n)))) (expm1 (* n (log1p (/ i n))))) (expm1 (* n (log1p (/ i n))))) (* (* (/ i n) (/ i n)) (/ i n))) |
| 80.0ms | (/ (* (- (* 100 (exp i)) 100) n) i) |
| 79.0ms | (- (log (/ -1 n)) (log (/ -1 i))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 9.6b
Found 4 expressions with local error:
| 19.3b | (cbrt (/ i n)) |
| 19.3b | (cbrt (/ i n)) |
| 19.3b | (cbrt (/ i n)) |
| 10.1b | (* n (log1p (/ i n))) |
4 calls. Slowest were:
| 1.0ms | (* n (log1p (/ i n))) |
| 1.0ms | (cbrt (/ i n)) |
| 0.0ms | (cbrt (/ i n)) |
4 calls. Slowest were:
| 91.0ms | (cbrt (/ i n)) |
| 70.0ms | (cbrt (/ i n)) |
| 62.0ms | (* n (log1p (/ i n))) |
| 61.0ms | (cbrt (/ i n)) |
30 calls. Slowest were:
| 232.0ms | (exp (* 1/3 (- (log (/ 1 n)) (log (/ 1 i))))) |
| 215.0ms | (exp (* 1/3 (- (log (/ 1 n)) (log (/ 1 i))))) |
| 178.0ms | (exp (* 1/3 (- (log (/ 1 n)) (log (/ 1 i))))) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 9.6b
Found 4 expressions with local error:
| 16.2b | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 16.2b | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 16.2b | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 10.1b | (* n (log1p (/ i n))) |
4 calls. Slowest were:
| 2.0ms | (* n (log1p (/ i n))) |
| 2.0ms | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 1.0ms | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
4 calls. Slowest were:
| 1.2s | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 1.1s | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 1.1s | (cbrt (/ (expm1 (* n (log1p (/ i n)))) (/ i n))) |
| 62.0ms | (* n (log1p (/ i n))) |
69 calls. Slowest were:
| 399.0ms | (exp (* 1/3 (- (+ (log (/ 1 i)) (log (- (exp i) 1))) (log (/ 1 n))))) |
| 338.0ms | (exp (* 1/3 (- (+ (log (/ 1 i)) (log (- (exp i) 1))) (log (/ 1 n))))) |
| 307.0ms | (exp (* 1/3 (- (+ (log (/ 1 i)) (log (- (exp i) 1))) (log (/ 1 n))))) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 9.6b
8.2% (10.5b remaining)
Error of 20.0b against oracle of 9.5b and baseline of 20.9b