



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.4 |
|---|---|
| Target | 20.4 |
| Herbie | 7.9 |
if b < -3.9662470494954583e+148Initial program 58.7
Simplified58.7
Taylor expanded around -inf 2.0
if -3.9662470494954583e+148 < b < 2.2500462861079674e-93Initial program 11.4
Simplified11.4
if 2.2500462861079674e-93 < b < 1.6889965443539738e+108Initial program 42.7
Simplified42.7
rmApplied flip--42.7
Applied associate-/l/45.2
Simplified16.7
rmApplied associate-/l*9.2
if 1.6889965443539738e+108 < b Initial program 59.4
Simplified59.4
Taylor expanded around inf 2.3
Simplified2.3
Final simplification7.9
herbie shell --seed 323378853
(FPCore (a b c)
:name "quadp (p42, positive)"
:herbie-target
(if (< b 0) (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) (/ c (* a (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))))
(/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
Time bar (total: 30.5s)Debug log
| 46.0ms | 199× | body | 80 | valid |
| 42.0ms | 30× | body | 2560 | valid |
| 29.0ms | 25× | body | 1280 | valid |
| 18.0ms | 84× | body | 80 | nan |
| 12.0ms | 13× | body | 640 | valid |
| 5.0ms | 7× | body | 320 | valid |
| 3.0ms | 2× | body | 5120 | valid |
| 3.0ms | 6× | body | 160 | valid |
| 119.0ms | (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.1b
Found 4 expressions with local error:
| 24.4b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.8b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 1.4b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 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 |
| 4× | flip-- |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | associate-/l* |
| 3× | sqrt-prod |
| 3× | times-frac |
| 2× | associate-/l/ |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 49.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 25.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 5.0ms | (- (* b b) (* (* 4 a) c)) |
| 5.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 116.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 89.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 83.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 51.0ms | (- (* b b) (* (* 4 a) c)) |
| 549.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 486.0ms | (- (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (* b b)) |
| 276.0ms | (- (pow (sqrt (- (* b b) (* (* 4 a) c))) 3) (pow b 3)) |
| 220.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
| 180.0ms | (- b (* 2 (/ (* a c) b))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 4.7b
Found 4 expressions with local error:
| 24.4b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 24.4b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.8b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 1.8b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 8× | add-sqr-sqrt |
| 6× | sqrt-prod |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 4× | add-log-exp |
| 4× | flip-- |
| 4× | add-exp-log |
| 4× | sqrt-div |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 4× | pow1 |
| 2× | pow1/2 |
| 2× | sub-neg |
| 2× | rem-sqrt-square |
| 28.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 25.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 5.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 5.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 89.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 89.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 81.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 80.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 551.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 551.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 217.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
| 212.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
| 181.0ms | (- b (* 2 (/ (* a c) b))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 4.7b
Found 4 expressions with local error:
| 25.2b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 25.2b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 10.8b | (- (* (fabs (cbrt (- (* b b) (* (* 4 a) c)))) (sqrt (cbrt (- (* b b) (* (* 4 a) c))))) b) |
| 1.4b | (- (* b b) (* (* 4 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 |
| 18.0ms | (- (* (fabs (cbrt (- (* b b) (* (* 4 a) c)))) (sqrt (cbrt (- (* b b) (* (* 4 a) c))))) b) |
| 5.0ms | (- (* b b) (* (* 4 a) c)) |
| 3.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 3.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 164.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 157.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 131.0ms | (- (* (fabs (cbrt (- (* b b) (* (* 4 a) c)))) (sqrt (cbrt (- (* b b) (* (* 4 a) c))))) b) |
| 52.0ms | (- (* b b) (* (* 4 a) c)) |
| 949.0ms | (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 572.0ms | (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 439.0ms | (- (exp (* 1/3 (+ (log -4) (+ (log a) (log c))))) (* 1/12 (* (exp (* 1/3 (- (log -4) (+ (* 2 (log a)) (* 2 (log c)))))) (pow b 2)))) |
| 393.0ms | (- (exp (* 1/3 (+ (log -4) (+ (log a) (log c))))) (* 1/12 (* (exp (* 1/3 (- (log -4) (+ (* 2 (log a)) (* 2 (log c)))))) (pow b 2)))) |
| 323.0ms | (- (* (fabs (pow (- (pow b 2) (* 4 (* a c))) 1/3)) (exp (* 1/6 (+ (log -4) (+ (log a) (log c)))))) (+ b (* 1/24 (* (pow b 2) (* (fabs (pow (- (pow b 2) (* 4 (* a c))) 1/3)) (exp (* 1/6 (- (log -4) (+ (* 5 (log a)) (* 5 (log c))))))))))) |
7 alts after pruning (4 fresh and 3 done)
Merged error: 4.7b
Found 4 expressions with local error:
| 24.4b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 15.7b | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 1.4b | (- (* b b) (* (* 4 a) c)) |
| 0.6b | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 9× | add-exp-log |
| 9× | add-cbrt-cube |
| 6× | add-log-exp |
| 6× | add-sqr-sqrt |
| 5× | add-cube-cbrt |
| 5× | *-un-lft-identity |
| 4× | pow1 |
| 3× | sqrt-prod |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | associate-/r/ |
| 2× | associate-*r/ |
| 2× | div-exp |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | flip3-+ |
| 2× | cbrt-undiv |
| 1× | associate-/l* |
| 1× | div-inv |
| 1× | cbrt-unprod |
| 1× | prod-exp |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | times-frac |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 33.0ms | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 18.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 5.0ms | (- (* b b) (* (* 4 a) c)) |
| 5.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 221.0ms | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 95.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 81.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 53.0ms | (- (* b b) (* (* 4 a) c)) |
| 1.2s | (/ (* (* (* (* -4 c) a) (* (* -4 c) a)) (* (* -4 c) a)) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (+ (sqrt (- (* b b) (* (* 4 a) c))) b) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)))) |
| 1.1s | (/ (* (* (* (* -4 c) a) (* (* -4 c) a)) (* (* -4 c) a)) (* (* (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)))) |
| 552.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 420.0ms | (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (- (* b b) (* (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 230.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
8 alts after pruning (4 fresh and 4 done)
Merged error: 4.7b
90.7% (2.6b remaining)
Error of 7.9b against oracle of 5.2b and baseline of 33.4b
| 2.0s | 5854× | body | 80 | valid |
| 1.7s | 1062× | body | 2560 | valid |
| 1.2s | 906× | body | 1280 | valid |
| 699.0ms | 2978× | body | 80 | nan |
| 446.0ms | 451× | body | 640 | valid |
| 271.0ms | 262× | body | 320 | valid |
| 182.0ms | 94× | body | 5120 | valid |
| 56.0ms | 122× | body | 160 | valid |