


Bits error versus a



Bits error versus b_2



Bits error versus c
if b_2 < -2.927037124445895e+148Initial program 58.7
Simplified58.7
Taylor expanded around -inf 2.1
Simplified2.1
if -2.927037124445895e+148 < b_2 < 3.587545867948427e-80Initial program 12.4
Simplified12.4
if 3.587545867948427e-80 < b_2 Initial program 52.5
Simplified52.5
Taylor expanded around inf 8.4
Final simplification9.7
herbie shell --seed 2019004 +o rules:numerics
(FPCore (a b_2 c)
:name "quad2p (problem 3.2.1, positive)"
(/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))
Time bar (total: 25.0s)Debug log
1 calls. Slowest were:
| 24.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 32.6b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 0.0b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
4 calls. Slowest were:
| 23.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 14.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 4.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 56.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 44.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 39.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 20.0ms | (- (* b_2 b_2) (* a c)) |
31 calls. Slowest were:
| 553.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 351.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 334.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.2b | (/ 1 (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 0.1b | (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
4 calls. Slowest were:
| 13.0ms | (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 3.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 2.0ms | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 137.0ms | (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 86.0ms | (/ 1 (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 60.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 34.0ms | (- (* b_2 b_2) (* a c)) |
34 calls. Slowest were:
| 484.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 385.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 329.0ms | (+ (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (+ (* b_2 b_2) (* (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
8 alts after pruning (6 fresh and 2 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.3b | (/ (/ 1 a) (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 0.1b | (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
4 calls. Slowest were:
| 7.0ms | (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 6.0ms | (/ (/ 1 a) (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 6.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 108.0ms | (/ (/ 1 a) (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 64.0ms | (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 37.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 20.0ms | (- (* b_2 b_2) (* a c)) |
87 calls. Slowest were:
| 1.6s | (/ (* (* (/ 1 a) (/ 1 a)) (/ 1 a)) (* (* (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) (/ 1 (- (sqrt (- (* b_2 b_2) (* a c))) b_2)))) |
| 570.0ms | (/ 1 (/ 1 (- (pow (sqrt (- (* b_2 b_2) (* a c))) 3) (pow b_2 3)))) |
| 564.0ms | (/ (sqrt (/ 1 a)) (/ 1 (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)))) |
9 alts after pruning (6 fresh and 3 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 22.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 3.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 3.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 2.0ms | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 54.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 44.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 26.0ms | (- (* b_2 b_2) (* a c)) |
| 22.0ms | (- (* b_2 b_2) (* a c)) |
30 calls. Slowest were:
| 434.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 349.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 346.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
9 alts after pruning (5 fresh and 4 done)
Merged error: 4.8b
84.7% (4.3b remaining)
Error of 9.7b against oracle of 5.4b and baseline of 33.5b