


Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -3.9662470494954583e+148Initial program 58.7
Simplified58.7
Taylor expanded around 0 58.7
Taylor expanded around -inf 2.0
if -3.9662470494954583e+148 < b_2 < 1.9841336792382835e-77Initial program 11.9
Simplified11.9
Taylor expanded around 0 11.9
if 1.9841336792382835e-77 < b_2 Initial program 52.4
Simplified52.4
Taylor expanded around 0 52.4
Taylor expanded around inf 8.8
Final simplification9.5
herbie shell --seed 323378853
(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: 21.9s)Debug log
| 38.0ms | 199× | body | 80 | valid |
| 35.0ms | 30× | body | 2560 | valid |
| 24.0ms | 25× | body | 1280 | valid |
| 14.0ms | 84× | body | 80 | nan |
| 10.0ms | 13× | body | 640 | valid |
| 5.0ms | 8× | body | 320 | valid |
| 3.0ms | 2× | body | 5120 | valid |
| 2.0ms | 5× | body | 160 | valid |
| 39.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.1b
Found 4 expressions with local error:
| 24.4b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.8b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.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 |
| 48.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 25.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)) |
| 103.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 81.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 78.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 42.0ms | (- (* b_2 b_2) (* a c)) |
| 421.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 385.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 347.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 185.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 135.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: 4.7b
Found 4 expressions with local error:
| 24.4b | (sqrt (- (pow b_2 2) (* a c))) |
| 1.8b | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 1.4b | (- (pow b_2 2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (pow b_2 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 |
| 73.0ms | (/ (- (sqrt (- (pow b_2 2) (* a c))) b_2) a) |
| 36.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 7.0ms | (- (pow b_2 2) (* a c)) |
| 7.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 105.0ms | (/ (- (sqrt (- (pow b_2 2) (* a c))) b_2) a) |
| 79.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 71.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 47.0ms | (- (pow b_2 2) (* a c)) |
| 562.0ms | (- (* (sqrt (- (pow b_2 2) (* a c))) (sqrt (- (pow b_2 2) (* a c)))) (* b_2 b_2)) |
| 504.0ms | (sqrt (- (pow (pow b_2 2) 3) (pow (* a c) 3))) |
| 346.0ms | (sqrt (- (* (pow b_2 2) (pow b_2 2)) (* (* a c) (* a c)))) |
| 181.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 142.0ms | (- (pow (sqrt (- (pow b_2 2) (* a c))) 3) (pow b_2 3)) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 4.7b
Found 4 expressions with local error:
| 24.4b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.8b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 1.4b | (- (* b_2 b_2) (* a c)) |
| 0.2b | (/ 1 (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 6× | associate-/r* |
| 5× | add-exp-log |
| 5× | pow1 |
| 4× | flip-- |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 3× | sqrt-prod |
| 3× | associate-/r/ |
| 2× | div-inv |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | inv-pow |
| 1× | pow-flip |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rec-exp |
| 1× | rem-sqrt-square |
| 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)) |
| 4.0ms | (/ 1 (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 140.0ms | (/ 1 (/ a (- (sqrt (- (* b_2 b_2) (* a c))) b_2))) |
| 80.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 74.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 41.0ms | (- (* b_2 b_2) (* a c)) |
| 402.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 349.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 201.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))) |
| 184.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 97.0ms | (- (* 1/2 (/ (* a c) b_2)) (* 2 b_2)) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 4.7b
Found 4 expressions with local error:
| 25.2b | (cbrt (- (* b_2 b_2) (* a c))) |
| 25.2b | (cbrt (- (* b_2 b_2) (* a c))) |
| 10.8b | (- (* (fabs (cbrt (- (* b_2 b_2) (* a c)))) (sqrt (cbrt (- (* b_2 b_2) (* a c))))) b_2) |
| 1.4b | (- (* b_2 b_2) (* a c)) |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 6× | add-sqr-sqrt |
| 6× | cbrt-prod |
| 4× | flip-- |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 4× | cbrt-div |
| 4× | pow1 |
| 2× | pow1/3 |
| 2× | sub-neg |
| 1× | diff-log |
| 20.0ms | (- (* (fabs (cbrt (- (* b_2 b_2) (* a c)))) (sqrt (cbrt (- (* b_2 b_2) (* a c))))) b_2) |
| 5.0ms | (- (* b_2 b_2) (* a c)) |
| 3.0ms | (cbrt (- (* b_2 b_2) (* a c))) |
| 3.0ms | (cbrt (- (* b_2 b_2) (* a c))) |
| 154.0ms | (cbrt (- (* b_2 b_2) (* a c))) |
| 154.0ms | (cbrt (- (* b_2 b_2) (* a c))) |
| 127.0ms | (- (* (fabs (cbrt (- (* b_2 b_2) (* a c)))) (sqrt (cbrt (- (* b_2 b_2) (* a c))))) b_2) |
| 42.0ms | (- (* b_2 b_2) (* a c)) |
| 405.0ms | (cbrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 401.0ms | (cbrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 380.0ms | (- (* (exp (* 1/6 (+ (log a) (+ (log -1) (log c))))) (fabs (pow (- (pow b_2 2) (* a c)) 1/3))) (+ b_2 (* 1/6 (* (exp (* 1/6 (- (log -1) (+ (* 5 (log a)) (* 5 (log c)))))) (* (pow b_2 2) (fabs (pow (- (pow b_2 2) (* a c)) 1/3))))))) |
| 345.0ms | (cbrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 339.0ms | (cbrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
9 alts after pruning (6 fresh and 3 done)
Merged error: 4.7b
84.6% (4.3b remaining)
Error of 9.5b against oracle of 5.2b and baseline of 33.4b
| 1.8s | 5856× | body | 80 | valid |
| 1.4s | 1063× | body | 2560 | valid |
| 948.0ms | 909× | body | 1280 | valid |
| 596.0ms | 2975× | body | 80 | nan |
| 384.0ms | 451× | body | 640 | valid |
| 184.0ms | 262× | body | 320 | valid |
| 132.0ms | 94× | body | 5120 | valid |
| 46.0ms | 119× | body | 160 | valid |