


Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -5.18597177717162e+27Initial program 35.2
Simplified35.2
Taylor expanded around -inf 6.4
if -5.18597177717162e+27 < b_2 < -3.4927265156327194e-305Initial program 10.7
Simplified10.7
rmApplied div-inv10.8
if -3.4927265156327194e-305 < b_2 < 2.65285615146505e+68Initial program 30.3
Simplified30.3
rmApplied div-inv30.4
rmApplied flip--30.5
Applied associate-*l/30.5
Simplified16.0
Taylor expanded around 0 9.3
Simplified9.3
if 2.65285615146505e+68 < b_2 Initial program 56.9
Simplified56.9
rmApplied div-inv56.9
Taylor expanded around inf 3.6
Final simplification7.5
herbie shell --seed 1137770271
(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: 28.2s)Debug log
| 39× | 45.0ms | (body real 2560) |
| 189× | 38.0ms | (body real 80) |
| 28× | 37.0ms | (body real 1280) |
| 93× | 17.0ms | (body nan 80) |
| 366× | 13.0ms | (pre true 80) |
| 9× | 7.0ms | (body real 640) |
| 3× | 4.0ms | (body real 5120) |
| 4× | 2.0ms | (body real 320) |
| 1× | 0.0ms | (body real 160) |
| 38.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 33.7b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.6b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.4b | (- (* 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 |
| 50.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 26.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 5.0ms | (- (* b_2 b_2) (* a c)) |
| 103.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 85.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 72.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 42.0ms | (- (* b_2 b_2) (* a c)) |
| 403.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 393.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 387.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 188.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 134.0ms | (- (pow (sqrt (- (* b_2 b_2) (* a c))) 3) (pow b_2 3)) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 6.7b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.6b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.4b | (- (* b_2 b_2) (* a c)) |
| 0.2b | (* (- (sqrt (- (* b_2 b_2) (* a c))) b_2) (/ 1 a)) |
| 8× | add-sqr-sqrt |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 6× | add-log-exp |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | pow1 |
| 5× | flip-- |
| 5× | flip3-- |
| 4× | associate-*r* |
| 3× | sqrt-prod |
| 3× | associate-*l* |
| 2× | associate-*l/ |
| 2× | sqrt-div |
| 2× | sub-neg |
| 2× | frac-times |
| 1× | div-inv |
| 1× | cbrt-unprod |
| 1× | *-commutative |
| 1× | un-div-inv |
| 1× | associate-*r/ |
| 1× | prod-exp |
| 1× | pow-prod-down |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | rem-sqrt-square |
| 137.0ms | (* (- (sqrt (- (* b_2 b_2) (* a c))) b_2) (/ 1 a)) |
| 58.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 10.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 9.0ms | (- (* b_2 b_2) (* a c)) |
| 160.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 153.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 137.0ms | (* (- (sqrt (- (* b_2 b_2) (* a c))) b_2) (/ 1 a)) |
| 58.0ms | (- (* b_2 b_2) (* a c)) |
| 997.0ms | (* (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) 1) |
| 692.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))) a) |
| 516.0ms | (* (* (* (- (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)) (* (* (/ 1 a) (/ 1 a)) (/ 1 a))) |
| 490.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 453.0ms | (* (- (pow (sqrt (- (* b_2 b_2) (* a c))) 3) (pow b_2 3)) (/ 1 a)) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b_2 b_2) (* a c))) |
| 15.3b | (/ (+ (- (* a c)) 0) a) |
| 0.6b | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.4b | (- (* 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× | associate-/l* |
| 3× | sqrt-prod |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | flip3-+ |
| 1× | div-inv |
| 1× | associate-+l- |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | neg-sub0 |
| 22.0ms | (/ (+ (- (* a c)) 0) a) |
| 19.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 8.0ms | (- (* b_2 b_2) (* a c)) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 80.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 72.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 47.0ms | (- (* b_2 b_2) (* a c)) |
| 13.0ms | (/ (+ (- (* a c)) 0) a) |
| 798.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 744.0ms | (+ (pow (- (* a c)) 3) (pow 0 3)) |
| 464.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 307.0ms | (- (* (- (* a c)) (- (* a c))) (* 0 0)) |
| 216.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
6 alts after pruning (5 fresh and 1 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b_2 b_2) (* a c))) |
| 0.6b | (+ (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)) |
| 15× | add-sqr-sqrt |
| 14× | add-cube-cbrt |
| 14× | *-un-lft-identity |
| 12× | times-frac |
| 6× | add-log-exp |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 4× | associate-/l* |
| 4× | neg-mul-1 |
| 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× | distribute-frac-neg |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 1× | neg-sub0 |
| 33.0ms | (/ (- c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 19.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 5.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 5.0ms | (- (* b_2 b_2) (* a c)) |
| 206.0ms | (/ (- c) (+ (sqrt (- (* b_2 b_2) (* a c))) b_2)) |
| 87.0ms | (+ (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 79.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 44.0ms | (- (* b_2 b_2) (* a c)) |
| 433.0ms | (/ (* (* (- c) (- c)) (- 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))) |
| 402.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 346.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 287.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))) |
| 191.0ms | (- (* 2 b_2) (* 1/2 (/ (* a c) b_2))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 6.6b
92.7% (1.7b remaining)
Error of 7.5b against oracle of 5.8b and baseline of 28.9b
| 1068× | 1.7s | (body real 2560) |
| 5866× | 1.3s | (body real 80) |
| 910× | 1.0s | (body real 1280) |
| 11702× | 870.0ms | (pre true 80) |
| 2943× | 625.0ms | (body nan 80) |
| 483× | 428.0ms | (body real 640) |
| 111× | 233.0ms | (body real 5120) |
| 220× | 137.0ms | (body real 320) |
| 101× | 40.0ms | (body real 160) |