


Bits error versus a



Bits error versus b_2



Bits error versus c
if b_2 < -4.0614473548971374e+148Initial program 58.9
Simplified58.9
rmApplied *-un-lft-identity58.9
Applied associate-/l*58.9
rmApplied associate-/r/58.9
rmApplied associate-*l/58.9
Simplified58.9
Taylor expanded around -inf 2.2
Simplified2.2
if -4.0614473548971374e+148 < b_2 < 1.2612246360442721e-22Initial program 14.7
Simplified14.7
rmApplied *-un-lft-identity14.7
Applied associate-/l*14.9
rmApplied associate-/r/14.9
rmApplied associate-*l/14.7
Simplified14.7
if 1.2612246360442721e-22 < b_2 Initial program 54.5
Simplified54.5
Taylor expanded around inf 6.4
Final simplification10.3
herbie shell --seed 2019007 +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: 18.8s)Debug log
1 calls. Slowest were:
| 48.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 33.9b
Found 4 expressions with local error:
| 23.9b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.5b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.0b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
4 calls. Slowest were:
| 49.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 25.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 8.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 69.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 65.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 42.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 34.0ms | (- (* b_2 b_2) (* a c)) |
31 calls. Slowest were:
| 414.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 395.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 313.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 23.9b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.5b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.0b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ 1 (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
4 calls. Slowest were:
| 21.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 6.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 2.0ms | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 109.0ms | (/ 1 (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 57.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 44.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 23.0ms | (- (* b_2 b_2) (* a c)) |
34 calls. Slowest were:
| 444.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 340.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 248.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))) |
11 alts after pruning (10 fresh and 1 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 23.9b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.5b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.0b | (- (* b_2 b_2) (* a c)) |
| 0.2b | (* (/ 1 a) (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
4 calls. Slowest were:
| 25.0ms | (* (/ 1 a) (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 13.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 3.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 119.0ms | (* (/ 1 a) (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 68.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 59.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 37.0ms | (- (* b_2 b_2) (* a c)) |
45 calls. Slowest were:
| 781.0ms | (* (* (* (/ 1 a) (/ 1 a)) (/ 1 a)) (* (* (- (sqrt (- (* b_2 b_2) (* a c))) b_2) (- (sqrt (- (* b_2 b_2) (* a c))) b_2)) (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 546.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 491.0ms | (* 1 (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2))) |
11 alts after pruning (9 fresh and 2 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 23.9b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.5b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.0b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
4 calls. Slowest were:
| 28.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 15.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 3.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 67.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 55.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 38.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 25.0ms | (- (* b_2 b_2) (* a c)) |
31 calls. Slowest were:
| 504.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 386.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 338.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
11 alts after pruning (8 fresh and 3 done)
Merged error: 6.1b
83.3% (4.8b remaining)
Error of 10.3b against oracle of 5.5b and baseline of 34.1b