


Bits error versus a



Bits error versus b_2



Bits error versus c
if b_2 < -2.8328889513248502e-121Initial program 50.8
Taylor expanded around -inf 11.6
if -2.8328889513248502e-121 < b_2 < 3.0964691378674393e+108Initial program 11.9
rmApplied *-un-lft-identity11.9
Applied associate-/l*12.0
if 3.0964691378674393e+108 < b_2 Initial program 46.3
Taylor expanded around inf 3.5
Simplified3.5
Final simplification10.6
herbie shell --seed 2018365 +o rules:numerics
(FPCore (a b_2 c)
:name "NMSE problem 3.2.1"
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))
Time bar (total: 31.9s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.5b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.3b | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
4 calls. Slowest were:
| 20.0ms | (/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
| 13.0ms | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 63.0ms | (/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
| 60.0ms | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 53.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 28.0ms | (- (* b_2 b_2) (* a c)) |
84 calls. Slowest were:
| 482.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 291.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 242.0ms | (- (* (- b_2) (- b_2)) (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c))))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 5.9b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (pow b_2 2) (* a c))) |
| 1.3b | (- (- b_2) (sqrt (- (pow b_2 2) (* a c)))) |
| 0.8b | (- (pow b_2 2) (* a c)) |
| 0.1b | (/ (- (- b_2) (sqrt (- (pow b_2 2) (* a c)))) a) |
4 calls. Slowest were:
| 22.0ms | (/ (- (- b_2) (sqrt (- (pow b_2 2) (* a c)))) a) |
| 13.0ms | (- (- b_2) (sqrt (- (pow b_2 2) (* a c)))) |
| 4.0ms | (sqrt (- (pow b_2 2) (* a c))) |
4 calls. Slowest were:
| 68.0ms | (/ (- (- b_2) (sqrt (- (pow b_2 2) (* a c)))) a) |
| 63.0ms | (- (- b_2) (sqrt (- (pow b_2 2) (* a c)))) |
| 42.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 27.0ms | (- (pow b_2 2) (* a c)) |
92 calls. Slowest were:
| 571.0ms | (sqrt (- (* (pow b_2 2) (pow b_2 2)) (* (* a c) (* a c)))) |
| 424.0ms | (sqrt (- (pow (pow b_2 2) 3) (pow (* a c) 3))) |
| 186.0ms | (- (* (- b_2) (- b_2)) (* (sqrt (- (pow b_2 2) (* a c))) (sqrt (- (pow b_2 2) (* a c))))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 5.9b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.3b | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.2b | (/ 1 (/ a (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))) |
4 calls. Slowest were:
| 12.0ms | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 105.0ms | (/ 1 (/ a (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))) |
| 63.0ms | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 40.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 34.0ms | (- (* b_2 b_2) (* a c)) |
85 calls. Slowest were:
| 546.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 354.0ms | (+ (* (- b_2) (- b_2)) (+ (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))) |
| 321.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 5.9b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.3b | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.3b | (/ (/ 1 a) (/ 1 (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))) |
4 calls. Slowest were:
| 7.0ms | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 4.0ms | (/ (/ 1 a) (/ 1 (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))) |
| 3.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 91.0ms | (/ (/ 1 a) (/ 1 (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))) |
| 53.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 51.0ms | (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) |
| 29.0ms | (- (* b_2 b_2) (* a c)) |
137 calls. Slowest were:
| 1.5s | (/ (* (* (/ 1 a) (/ 1 a)) (/ 1 a)) (* (* (/ 1 (- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))) (/ 1 (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))) (/ 1 (- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))) |
| 771.0ms | (/ (cbrt (/ 1 a)) (+ (* (- b_2) (- b_2)) (+ (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))) |
| 731.0ms | (/ (sqrt (/ 1 a)) (+ (* (- b_2) (- b_2)) (+ (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 5.9b
82.3% (5.0b remaining)
Error of 10.6b against oracle of 5.6b and baseline of 33.6b