


Bits error versus a



Bits error versus b_2



Bits error versus c
if b_2 < -2.074336923720548e+100Initial program 45.9
Initial simplification45.9
rmApplied div-inv45.9
Taylor expanded around -inf 3.6
Simplified3.6
if -2.074336923720548e+100 < b_2 < -4.1996115297456014e-247Initial program 7.9
Initial simplification7.9
rmApplied div-inv8.0
if -4.1996115297456014e-247 < b_2 < 1.8812121325960807e+75Initial program 27.8
Initial simplification27.8
rmApplied div-inv27.8
rmApplied flip--28.0
Applied associate-*l/28.0
Simplified15.5
Taylor expanded around -inf 8.9
Simplified8.9
if 1.8812121325960807e+75 < b_2 Initial program 57.7
Initial simplification57.7
Taylor expanded around inf 3.7
Final simplification6.5
herbie shell --seed 2018362 +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: 16.8s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 30.9b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.8b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.4b | (- (* b_2 b_2) (* a c)) |
| 0.4b | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
4 calls. Slowest were:
| 44.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 24.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 8.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
4 calls. Slowest were:
| 82.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 66.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 64.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 30.0ms | (- (* b_2 b_2) (* a c)) |
31 calls. Slowest were:
| 494.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 466.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 331.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 4.2b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.8b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.5b | (* (- (sqrt (- (* b_2 b_2) (* a c))) b_2) (/ 1 a)) |
| 0.4b | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 48.0ms | (* (- (sqrt (- (* b_2 b_2) (* a c))) b_2) (/ 1 a)) |
| 22.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 66.0ms | (* (- (sqrt (- (* b_2 b_2) (* a c))) b_2) (/ 1 a)) |
| 62.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 49.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 38.0ms | (- (* b_2 b_2) (* a c)) |
41 calls. Slowest were:
| 777.0ms | (* (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) 1) |
| 501.0ms | (* (- (pow (sqrt (- (* b_2 b_2) (* a c))) 3) (pow b_2 3)) (/ 1 a)) |
| 473.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 4.2b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b_2 b_2) (* a c))) |
| 10.6b | (/ (fma a (- c) 0) a) |
| 0.9b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.4b | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 17.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 4.0ms | (/ (fma a (- c) 0) a) |
4 calls. Slowest were:
| 49.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 44.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 43.0ms | (/ (fma a (- c) 0) a) |
| 36.0ms | (- (* b_2 b_2) (* a c)) |
28 calls. Slowest were:
| 504.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 387.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 127.0ms | (- (* 2 b_2) (* 1/2 (/ (* a c) b_2))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 4.2b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.9b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.4b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
4 calls. Slowest were:
| 14.0ms | (/ (- c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 11.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
4 calls. Slowest were:
| 157.0ms | (/ (- c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 57.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 48.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 31.0ms | (- (* b_2 b_2) (* a c)) |
61 calls. Slowest were:
| 512.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 451.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))) |
| 389.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
8 alts after pruning (6 fresh and 2 done)
Merged error: 4.2b
94.6% (1.3b remaining)
Error of 6.5b against oracle of 5.2b and baseline of 28.7b