


Bits error versus a



Bits error versus b_2



Bits error versus c
if b_2 < -4.50441748916982e-54Initial program 26.0
Simplified26.0
Taylor expanded around -inf 10.3
Simplified10.3
if -4.50441748916982e-54 < b_2 < 5.148536202069843e+105Initial program 24.7
Simplified24.7
rmApplied flip--28.5
Applied associate-/l/33.4
Simplified22.8
rmApplied distribute-frac-neg22.8
Simplified13.6
if 5.148536202069843e+105 < b_2 Initial program 58.9
Simplified58.9
Taylor expanded around inf 2.5
Final simplification10.2
herbie shell --seed 2019002 +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: 20.7s)Debug log
1 calls. Slowest were:
| 49.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 37.0b
Found 4 expressions with local error:
| 23.6b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.1b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.6b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
4 calls. Slowest were:
| 32.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:
| 73.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 58.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 48.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 29.0ms | (- (* b_2 b_2) (* a c)) |
31 calls. Slowest were:
| 548.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 454.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 426.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 4.5b
Found 4 expressions with local error:
| 23.6b | (sqrt (- (* b_2 b_2) (* a c))) |
| 14.2b | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 1.4b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.6b | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 19.0ms | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 17.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 117.0ms | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 43.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 41.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 33.0ms | (- (* b_2 b_2) (* a c)) |
48 calls. Slowest were:
| 516.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 453.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))) |
| 429.0ms | (/ (* (* (- (* a c)) (- (* a c))) (- (* a c))) (* (* (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 4.5b
Found 4 expressions with local error:
| 23.6b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.4b | (+ b_2 (sqrt (- (* b_2 b_2) (* a c)))) |
| 0.6b | (- (* b_2 b_2) (* a c)) |
| 0.0b | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
4 calls. Slowest were:
| 3.0ms | (- (* b_2 b_2) (* a c)) |
| 3.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 3.0ms | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
4 calls. Slowest were:
| 111.0ms | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
| 56.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 40.0ms | (+ b_2 (sqrt (- (* b_2 b_2) (* a c)))) |
| 20.0ms | (- (* b_2 b_2) (* a c)) |
26 calls. Slowest were:
| 497.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 329.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)))))) |
| 310.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 4.5b
Found 4 expressions with local error:
| 23.6b | (sqrt (- (* b_2 b_2) (* a c))) |
| 7.6b | (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 6.1b | (* (/ -1 a) (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 1.4b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
4 calls. Slowest were:
| 28.0ms | (* (/ -1 a) (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 18.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 12.0ms | (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
4 calls. Slowest were:
| 130.0ms | (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 111.0ms | (* (/ -1 a) (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 47.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 41.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
56 calls. Slowest were:
| 518.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 454.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))) |
| 415.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))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 4.5b
80.6% (4.6b remaining)
Error of 10.2b against oracle of 5.6b and baseline of 29.3b