


Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -1.5880262840008754e+117Initial program 48.6
Simplified48.6
Taylor expanded around -inf 3.1
if -1.5880262840008754e+117 < b_2 < -2.4925302910670146e-306Initial program 9.1
Simplified9.1
rmApplied add-sqr-sqrt9.1
Applied sqrt-prod9.3
if -2.4925302910670146e-306 < b_2 < 3.1234815487678856e+63Initial program 29.9
Simplified29.9
rmApplied flip--30.1
Applied associate-/l/34.7
Simplified21.5
rmApplied distribute-frac-neg21.5
Simplified9.2
if 3.1234815487678856e+63 < b_2 Initial program 56.6
Simplified56.6
Taylor expanded around inf 3.5
Final simplification6.8
herbie shell --seed 2019016
(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.7s)Debug log
| 378× | (pre true 80) |
| 186× | (body real 80) |
| 93× | (body nan 80) |
| 36× | (body real 2560) |
| 30× | (body real 1280) |
| 17× | (body real 640) |
| 9× | (body real 160) |
| 5× | (body real 320) |
| 2× | (body real 5120) |
| 386.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.8b
Found 4 expressions with local error:
| 23.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 2.2b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.2b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 4× | flip-- |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | associate-/l* |
| 3× | sqrt-prod |
| 2× | associate-/l/ |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 79.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 41.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 28.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
| 61.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 55.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 45.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 21.0ms | (- (* b_2 b_2) (* a c)) |
| 325.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 186.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 163.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 139.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 88.0ms | (- (pow (sqrt (- (* b_2 b_2) (* a c))) 3) (pow b_2 3)) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 5.0b
Found 4 expressions with local error:
| 23.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 16.1b | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 1.2b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.2b | (- (* b_2 b_2) (* a c)) |
| 8× | add-sqr-sqrt |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 6× | add-log-exp |
| 6× | associate-/l* |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | times-frac |
| 4× | pow1 |
| 3× | sqrt-prod |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | associate-/r/ |
| 2× | associate-*r/ |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | flip3-+ |
| 2× | distribute-rgt-neg-in |
| 2× | neg-mul-1 |
| 2× | distribute-lft-neg-in |
| 1× | div-inv |
| 1× | distribute-frac-neg |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 1× | neg-sub0 |
| 10.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 8.0ms | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 3.0ms | (- (* b_2 b_2) (* a c)) |
| 2.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 111.0ms | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 48.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 42.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 20.0ms | (- (* b_2 b_2) (* a c)) |
| 433.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)))) |
| 272.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 246.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 228.0ms | (/ (* a c) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 168.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))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 5.0b
Found 4 expressions with local error:
| 23.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.2b | (+ b_2 (sqrt (- (* b_2 b_2) (* a c)))) |
| 1.2b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | sqrt-prod |
| 3× | associate-/r* |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | associate-/r/ |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | flip3-+ |
| 1× | div-inv |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 2.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 2.0ms | (- (* b_2 b_2) (* a c)) |
| 1.0ms | (+ b_2 (sqrt (- (* b_2 b_2) (* a c)))) |
| 1.0ms | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
| 137.0ms | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
| 52.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 51.0ms | (+ b_2 (sqrt (- (* b_2 b_2) (* a c)))) |
| 31.0ms | (- (* b_2 b_2) (* a c)) |
| 300.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 245.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 198.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)))))) |
| 147.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 135.0ms | (- (* 2 b_2) (* 1/2 (/ (* a c) b_2))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 5.0b
Found 4 expressions with local error:
| 23.0b | (sqrt (- (* b_2 b_2) (* a c))) |
| 10.2b | (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 3.5b | (* (/ -1 a) (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 1.2b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 11× | add-exp-log |
| 11× | add-cbrt-cube |
| 11× | add-sqr-sqrt |
| 10× | add-cube-cbrt |
| 10× | *-un-lft-identity |
| 9× | associate-*r* |
| 6× | times-frac |
| 6× | pow1 |
| 4× | add-log-exp |
| 4× | associate-/r/ |
| 4× | associate-*l* |
| 3× | div-inv |
| 3× | flip-+ |
| 3× | sqrt-prod |
| 3× | flip3-+ |
| 3× | associate-/r* |
| 2× | cbrt-unprod |
| 2× | prod-exp |
| 2× | div-exp |
| 2× | sqrt-div |
| 2× | cbrt-undiv |
| 1× | associate-/l* |
| 1× | flip-- |
| 1× | *-commutative |
| 1× | associate-*r/ |
| 1× | associate-*l/ |
| 1× | pow-prod-down |
| 1× | pow1/2 |
| 1× | flip3-- |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | frac-times |
| 27.0ms | (* (/ -1 a) (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 19.0ms | (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 14.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 2.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 146.0ms | (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 102.0ms | (* (/ -1 a) (/ (* a c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 66.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 35.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 217.0ms | (* (* (* (/ -1 a) (/ -1 a)) (/ -1 a)) (/ (* (* (* a c) (* a c)) (* a c)) (* (* (+ (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)))) |
| 217.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 206.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 203.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))) |
| 194.0ms | (+ (log (/ -1 a)) (- (log (* a c)) (log (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 5.0b
95% (1.2b remaining)
Error of 6.8b against oracle of 5.7b and baseline of 29.6b
| 11497× | (pre true 80) |
| 5928× | (body real 80) |
| 2743× | (body nan 80) |
| 1015× | (body real 2560) |
| 886× | (body real 1280) |
| 456× | (body real 640) |
| 225× | (body real 320) |
| 126× | (body real 160) |
| 118× | (body real 5120) |