


Bits error versus a



Bits error versus k



Bits error versus m
if k < 2.9740929793241473e+109Initial program 0.0
Initial simplification0.0
rmApplied add-sqr-sqrt0.0
Applied associate-*l*0.0
if 2.9740929793241473e+109 < k Initial program 8.4
Initial simplification8.4
Taylor expanded around -inf 63.0
Simplified0.7
Final simplification0.2
herbie shell --seed 2018362 +o rules:numerics
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
(/ (* a (pow k m)) (+ (+ 1 (* 10 k)) (* k k))))
Time bar (total: 21.7s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 2.3b
Found 2 expressions with local error:
| 2.3b | (/ (* (pow k m) a) (fma k (+ k 10) 1)) |
| 0.0b | (* (pow k m) a) |
2 calls. Slowest were:
| 10.0ms | (/ (* (pow k m) a) (fma k (+ k 10) 1)) |
| 4.0ms | (* (pow k m) a) |
2 calls. Slowest were:
| 132.0ms | (/ (* (pow k m) a) (fma k (+ k 10) 1)) |
| 110.0ms | (* (pow k m) a) |
20 calls. Slowest were:
| 922.0ms | (/ (* (* (* (pow k m) a) (* (pow k m) a)) (* (pow k m) a)) (* (* (fma k (+ k 10) 1) (fma k (+ k 10) 1)) (fma k (+ k 10) 1))) |
| 594.0ms | (- (+ (* 99 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 4))) (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 2))) (* 10 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 3)))) |
| 293.0ms | (- (+ a (* (log k) (* m a))) (* 10 (* a k))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.0b
Found 3 expressions with local error:
| 2.3b | (/ (* (sqrt (pow k m)) (* (sqrt (pow k m)) a)) (fma k (+ k 10) 1)) |
| 0.0b | (* (sqrt (pow k m)) (* (sqrt (pow k m)) a)) |
| 0.0b | (* (sqrt (pow k m)) a) |
3 calls. Slowest were:
| 11.0ms | (* (sqrt (pow k m)) (* (sqrt (pow k m)) a)) |
| 11.0ms | (/ (* (sqrt (pow k m)) (* (sqrt (pow k m)) a)) (fma k (+ k 10) 1)) |
| 6.0ms | (* (sqrt (pow k m)) a) |
3 calls. Slowest were:
| 177.0ms | (* (sqrt (pow k m)) a) |
| 113.0ms | (/ (* (sqrt (pow k m)) (* (sqrt (pow k m)) a)) (fma k (+ k 10) 1)) |
| 113.0ms | (* (sqrt (pow k m)) (* (sqrt (pow k m)) a)) |
36 calls. Slowest were:
| 612.0ms | (- (+ (* 99 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 4))) (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 2))) (* 10 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 3)))) |
| 486.0ms | (* (* (* (sqrt (pow k m)) (sqrt (pow k m))) (sqrt (pow k m))) (* (* (* (sqrt (pow k m)) a) (* (sqrt (pow k m)) a)) (* (sqrt (pow k m)) a))) |
| 338.0ms | (- (+ (* 99 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 4))) (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 2))) (* 10 (/ (* a (exp (* m (- (log -1) (log (/ -1 k)))))) (pow k 3)))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 14.0b | (sqrt (fma k (+ k 10) 1)) |
| 14.0b | (sqrt (fma k (+ k 10) 1)) |
| 0.0b | (* (/ (pow k m) (sqrt (fma k (+ k 10) 1))) (/ a (sqrt (fma k (+ k 10) 1)))) |
| 0.0b | (/ (pow k m) (sqrt (fma k (+ k 10) 1))) |
4 calls. Slowest were:
| 7.0ms | (* (/ (pow k m) (sqrt (fma k (+ k 10) 1))) (/ a (sqrt (fma k (+ k 10) 1)))) |
| 4.0ms | (/ (pow k m) (sqrt (fma k (+ k 10) 1))) |
| 1.0ms | (sqrt (fma k (+ k 10) 1)) |
4 calls. Slowest were:
| 105.0ms | (* (/ (pow k m) (sqrt (fma k (+ k 10) 1))) (/ a (sqrt (fma k (+ k 10) 1)))) |
| 64.0ms | (/ (pow k m) (sqrt (fma k (+ k 10) 1))) |
| 25.0ms | (sqrt (fma k (+ k 10) 1)) |
| 24.0ms | (sqrt (fma k (+ k 10) 1)) |
102 calls. Slowest were:
| 660.0ms | (+ (- (log (pow k m)) (log (sqrt (fma k (+ k 10) 1)))) (log (/ a (sqrt (fma k (+ k 10) 1))))) |
| 611.0ms | (- (+ (* 99 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 4))) (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 2))) (* 10 (/ (* (exp (* -1 (* (log (/ 1 k)) m))) a) (pow k 3)))) |
| 473.0ms | (* (* (* (/ (pow k m) (sqrt (fma k (+ k 10) 1))) (/ (pow k m) (sqrt (fma k (+ k 10) 1)))) (/ (pow k m) (sqrt (fma k (+ k 10) 1)))) (* (* (/ a (sqrt (fma k (+ k 10) 1))) (/ a (sqrt (fma k (+ k 10) 1)))) (/ a (sqrt (fma k (+ k 10) 1))))) |
3 alts after pruning (2 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 14.4b | (cbrt (fma k (+ k 10) 1)) |
| 14.4b | (cbrt (fma (+ 10 k) k 1)) |
| 14.0b | (sqrt (fma k (+ k 10) 1)) |
| 0.2b | (* (fabs (cbrt (fma (+ 10 k) k 1))) (sqrt (cbrt (fma k (+ k 10) 1)))) |
4 calls. Slowest were:
| 6.0ms | (* (fabs (cbrt (fma (+ 10 k) k 1))) (sqrt (cbrt (fma k (+ k 10) 1)))) |
| 1.0ms | (sqrt (fma k (+ k 10) 1)) |
| 1.0ms | (cbrt (fma k (+ k 10) 1)) |
4 calls. Slowest were:
| 51.0ms | (cbrt (fma k (+ k 10) 1)) |
| 46.0ms | (cbrt (fma (+ 10 k) k 1)) |
| 40.0ms | (* (fabs (cbrt (fma (+ 10 k) k 1))) (sqrt (cbrt (fma k (+ k 10) 1)))) |
| 22.0ms | (sqrt (fma k (+ k 10) 1)) |
36 calls. Slowest were:
| 420.0ms | (- (+ (fabs (pow (fma (+ 10 k) k 1) 1/3)) (* 5/3 (* (fabs (pow (fma (+ 10 k) k 1) 1/3)) k))) (* 61/9 (* (fabs (pow (fma (+ 10 k) k 1) 1/3)) (pow k 2)))) |
| 304.0ms | (- (+ (* (fabs (pow (fma (+ 10 k) k 1) 1/3)) (pow k 1/3)) (* 5/3 (* (fabs (pow (fma (+ 10 k) k 1) 1/3)) (pow (/ 1 (pow k 2)) 1/3)))) (* 61/9 (* (fabs (pow (fma (+ 10 k) k 1) 1/3)) (pow (/ 1 (pow k 5)) 1/3)))) |
| 270.0ms | (- (+ (* (pow (* k -1) 1/3) (fabs (pow (fma (+ 10 k) k 1) 1/3))) (* 5/3 (* (pow (/ -1 (pow k 2)) 1/3) (fabs (pow (fma (+ 10 k) k 1) 1/3))))) (* 61/9 (* (pow (/ -1 (pow k 5)) 1/3) (fabs (pow (fma (+ 10 k) k 1) 1/3))))) |
3 alts after pruning (1 fresh and 2 done)
Merged error: 0.0b
97.2% (0.1b remaining)
Error of 0.2b against oracle of 0.1b and baseline of 2.1b