



Bits error versus x
Results
| Original | 52.2 |
|---|---|
| Target | 44.6 |
| Herbie | 0.1 |
if x < -1.0435353126017766Initial program 61.7
Simplified60.9
Taylor expanded around -inf 0.2
Simplified0.2
if -1.0435353126017766 < x < 0.007938820784157252Initial program 59.0
Simplified59.0
Taylor expanded around 0 0.1
if 0.007938820784157252 < x Initial program 30.1
Simplified0.0
Final simplification0.1
herbie shell --seed 2019002 +o rules:numerics
(FPCore (x)
:name "Hyperbolic arcsine"
:herbie-target
(if (< x 0) (log (/ -1 (- x (sqrt (+ (* x x) 1))))) (log (+ x (sqrt (+ (* x x) 1)))))
(log (+ x (sqrt (+ (* x x) 1)))))
Time bar (total: 19.1s)Debug log
1 calls. Slowest were:
| 4.0ms | (log (+ x (sqrt (+ (* x x) 1)))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 44.9b
Found 3 expressions with local error:
| 3.2b | (log (+ x (hypot 1 x))) |
| 0.3b | (+ x (hypot 1 x)) |
| 0.0b | (hypot 1 x) |
3 calls. Slowest were:
| 1.0ms | (log (+ x (hypot 1 x))) |
| 1.0ms | (+ x (hypot 1 x)) |
| 0.0ms | (hypot 1 x) |
3 calls. Slowest were:
| 137.0ms | (log (+ x (hypot 1 x))) |
| 31.0ms | (hypot 1 x) |
| 30.0ms | (+ x (hypot 1 x)) |
15 calls. Slowest were:
| 361.0ms | (- (+ (log 2) (* 1/4 (/ 1 (pow x 2)))) (+ (log (/ 1 x)) (* 3/32 (/ 1 (pow x 4))))) |
| 295.0ms | (- (* 1/8 (/ 1 (pow x 3))) (+ (* 1/16 (/ 1 (pow x 5))) (* 1/2 (/ 1 x)))) |
| 214.0ms | (- (* 1/8 (/ 1 (pow x 3))) (+ x (* 1/2 (/ 1 x)))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0.0b
Found 3 expressions with local error:
| 18.3b | (- (+ x (* 3/40 (pow x 5))) (* 1/6 (pow x 3))) |
| 0.1b | (* 1/6 (pow x 3)) |
| 0.0b | (* 3/40 (pow x 5)) |
3 calls. Slowest were:
| 10.0ms | (- (+ x (* 3/40 (pow x 5))) (* 1/6 (pow x 3))) |
| 2.0ms | (* 1/6 (pow x 3)) |
| 1.0ms | (* 3/40 (pow x 5)) |
3 calls. Slowest were:
| 24.0ms | (- (+ x (* 3/40 (pow x 5))) (* 1/6 (pow x 3))) |
| 21.0ms | (* 3/40 (pow x 5)) |
| 19.0ms | (* 1/6 (pow x 3)) |
25 calls. Slowest were:
| 194.0ms | (fma 1 (+ x (* 3/40 (pow x 5))) (- (* (pow x 3) 1/6))) |
| 128.0ms | (fma (sqrt (+ x (* 3/40 (pow x 5)))) (sqrt (+ x (* 3/40 (pow x 5)))) (- (* (pow x 3) 1/6))) |
| 125.0ms | (/ (exp (+ x (* 3/40 (pow x 5)))) (exp (* 1/6 (pow x 3)))) |
4 alts after pruning (2 fresh and 2 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 3.2b | (log (* (cbrt (+ x (hypot 1 x))) (cbrt (+ x (hypot 1 x))))) |
| 3.2b | (log (cbrt (+ x (hypot 1 x)))) |
| 0.3b | (+ x (hypot 1 x)) |
| 0.3b | (+ x (hypot 1 x)) |
4 calls. Slowest were:
| 8.0ms | (log (* (cbrt (+ x (hypot 1 x))) (cbrt (+ x (hypot 1 x))))) |
| 3.0ms | (log (cbrt (+ x (hypot 1 x)))) |
| 1.0ms | (+ x (hypot 1 x)) |
4 calls. Slowest were:
| 293.0ms | (log (* (cbrt (+ x (hypot 1 x))) (cbrt (+ x (hypot 1 x))))) |
| 228.0ms | (log (cbrt (+ x (hypot 1 x)))) |
| 34.0ms | (+ x (hypot 1 x)) |
| 30.0ms | (+ x (hypot 1 x)) |
35 calls. Slowest were:
| 336.0ms | (- (+ (* 1/3 (log (/ -1 x))) (+ (* 1/3 (log 1/2)) (* 1/32 (/ 1 (pow x 4))))) (* 1/12 (/ 1 (pow x 2)))) |
| 304.0ms | (- (* 1/8 (/ 1 (pow x 3))) (+ (* 1/16 (/ 1 (pow x 5))) (* 1/2 (/ 1 x)))) |
| 273.0ms | (- (* 1/8 (/ 1 (pow x 3))) (+ (* 1/16 (/ 1 (pow x 5))) (* 1/2 (/ 1 x)))) |
4 alts after pruning (2 fresh and 2 done)
Merged error: 0b
Found 4 expressions with local error:
| 3.3b | (log (sqrt (cbrt (+ x (hypot 1 x))))) |
| 3.3b | (log (sqrt (cbrt (+ x (hypot 1 x))))) |
| 3.2b | (log (* (cbrt (+ x (hypot 1 x))) (cbrt (+ x (hypot 1 x))))) |
| 0.3b | (+ x (hypot 1 x)) |
4 calls. Slowest were:
| 6.0ms | (log (* (cbrt (+ x (hypot 1 x))) (cbrt (+ x (hypot 1 x))))) |
| 2.0ms | (log (sqrt (cbrt (+ x (hypot 1 x))))) |
| 2.0ms | (log (sqrt (cbrt (+ x (hypot 1 x))))) |
4 calls. Slowest were:
| 285.0ms | (log (sqrt (cbrt (+ x (hypot 1 x))))) |
| 239.0ms | (log (sqrt (cbrt (+ x (hypot 1 x))))) |
| 209.0ms | (log (* (cbrt (+ x (hypot 1 x))) (cbrt (+ x (hypot 1 x))))) |
| 42.0ms | (+ x (hypot 1 x)) |
50 calls. Slowest were:
| 392.0ms | (- (+ (* 1/6 (log (/ -1 x))) (+ (* 1/6 (log 1/2)) (* 1/64 (/ 1 (pow x 4))))) (* 1/24 (/ 1 (pow x 2)))) |
| 338.0ms | (- (+ (* 1/6 (log (/ -1 x))) (+ (* 1/6 (log 1/2)) (* 1/64 (/ 1 (pow x 4))))) (* 1/24 (/ 1 (pow x 2)))) |
| 325.0ms | (- (+ (* 1/6 (log 2)) (* 1/24 (/ 1 (pow x 2)))) (+ (* 1/6 (log (/ 1 x))) (* 1/64 (/ 1 (pow x 4))))) |
4 alts after pruning (1 fresh and 3 done)
Merged error: 0b
99.7% (0.1b remaining)
Error of 0.1b against oracle of 0.0b and baseline of 31.1b