


Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -6.952882728129189e+102Initial program 44.6
Simplified44.6
Taylor expanded around -inf 3.6
if -6.952882728129189e+102 < b_2 < -3.033451162286551e-263Initial program 8.3
Simplified8.3
rmApplied div-inv8.5
if -3.033451162286551e-263 < b_2 < 2.1436459426001753e+111Initial program 30.3
Simplified30.3
rmApplied flip--30.4
Applied associate-/l/34.9
Simplified19.4
rmApplied distribute-frac-neg19.4
Simplified8.5
if 2.1436459426001753e+111 < b_2 Initial program 59.5
Simplified59.5
Taylor expanded around inf 2.8
Final simplification6.4
herbie shell --seed 1042149663
(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.3s)Debug log
| 369× | (pre true 80) |
| 181× | (body real 80) |
| 88× | (body nan 80) |
| 31× | (body real 1280) |
| 28× | (body real 2560) |
| 22× | (body real 640) |
| 13× | (body real 320) |
| 3× | (body real 160) |
| 3× | (body real 5120) |
| 41.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.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.8b | (- (* 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 |
| 46.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 24.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
| 101.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 79.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 77.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 41.0ms | (- (* b_2 b_2) (* a c)) |
| 412.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 406.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 386.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 190.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 163.0ms | (- (pow (sqrt (- (* b_2 b_2) (* a c))) 3) (pow b_2 3)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 5.5b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b_2 b_2) (* a c))) |
| 17.9b | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.4b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 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 |
| 21.0ms | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 18.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
| 210.0ms | (/ (- (* a c)) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 85.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 74.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 52.0ms | (- (* b_2 b_2) (* a c)) |
| 639.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)))) |
| 467.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 418.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 403.0ms | (- (log (- (* a c))) (log (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)))) |
| 323.0ms | (/ (* a c) (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 5.4b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.4b | (+ b_2 (sqrt (- (* 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 |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
| 3.0ms | (+ b_2 (sqrt (- (* b_2 b_2) (* a c)))) |
| 3.0ms | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
| 200.0ms | (/ c (+ b_2 (sqrt (- (* b_2 b_2) (* a c))))) |
| 82.0ms | (+ b_2 (sqrt (- (* b_2 b_2) (* a c)))) |
| 70.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 41.0ms | (- (* b_2 b_2) (* a c)) |
| 413.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 356.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 277.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)))))) |
| 198.0ms | (- (* 2 b_2) (* 1/2 (/ (* a c) b_2))) |
| 192.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
8 alts after pruning (6 fresh and 2 done)
Merged error: 5.4b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b_2 b_2) (* a c))) |
| 17.9b | (/ (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a c)) |
| 0.8b | (- (* b_2 b_2) (* a c)) |
| 0.5b | (/ -1 (/ (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a c))) |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | *-un-lft-identity |
| 6× | associate-/r* |
| 4× | pow1 |
| 3× | div-inv |
| 3× | sqrt-prod |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | associate-*r/ |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | times-frac |
| 2× | frac-2neg |
| 2× | clear-num |
| 1× | associate-/l* |
| 1× | flip-+ |
| 1× | associate-/r/ |
| 1× | div-exp |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | flip3-+ |
| 1× | sub-neg |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 29.0ms | (/ (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a c)) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 4.0ms | (- (* b_2 b_2) (* a c)) |
| 1.0ms | (/ -1 (/ (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a c))) |
| 189.0ms | (/ -1 (/ (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a c))) |
| 115.0ms | (/ (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a c)) |
| 70.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 42.0ms | (- (* b_2 b_2) (* a c)) |
| 710.0ms | (- (log (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2))) (log (* a c))) |
| 427.0ms | (* a (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2))) |
| 418.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 378.0ms | (/ (* (* (* 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))) (* (* (* a c) (* a c)) (* a c))) |
| 375.0ms | (/ (* a (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) (* a c)) |
9 alts after pruning (6 fresh and 3 done)
Merged error: 5.4b
95.7% (1.0b remaining)
Error of 6.4b against oracle of 5.4b and baseline of 29.1b
| 11580× | (pre true 80) |
| 5842× | (body real 80) |
| 2894× | (body nan 80) |
| 1079× | (body real 2560) |
| 861× | (body real 1280) |
| 446× | (body real 640) |
| 230× | (body real 320) |
| 132× | (body real 160) |
| 96× | (body real 5120) |