



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.4 |
|---|---|
| Target | 20.6 |
| Herbie | 7.4 |
if b < -2.1317624771383882e+58Initial program 37.4
Simplified37.4
Taylor expanded around -inf 5.4
if -2.1317624771383882e+58 < b < 1.080192248888389e-171Initial program 11.7
Simplified11.7
rmApplied div-sub11.7
if 1.080192248888389e-171 < b < 1.3660877782595927e+23Initial program 32.3
Simplified32.3
rmApplied flip--32.4
Applied associate-/l/37.4
Simplified20.4
rmApplied associate-/r*14.6
Taylor expanded around 0 5.8
rmApplied div-inv6.0
if 1.3660877782595927e+23 < b Initial program 55.5
Simplified55.5
Taylor expanded around inf 4.4
Simplified4.4
Final simplification7.4
herbie shell --seed 2019022
(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: 19.0s)Debug log
| 45.0ms | 40× | body | 2560 | valid |
| 20.0ms | 185× | body | 80 | valid |
| 14.0ms | 27× | body | 1280 | valid |
| 10.0ms | 97× | body | 80 | nan |
| 7.0ms | 15× | body | 640 | valid |
| 3.0ms | 7× | body | 320 | valid |
| 2.0ms | 3× | body | 5120 | valid |
| 1.0ms | 6× | body | 160 | valid |
| 54.0ms | (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 32.5b
Found 4 expressions with local error:
| 21.7b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 2.0b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 0.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 |
| 44.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 10.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 2.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 2.0ms | (- (* b b) (* (* 4 a) c)) |
| 66.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 59.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 54.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 33.0ms | (- (* b b) (* (* 4 a) c)) |
| 392.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 306.0ms | (- (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (* b b)) |
| 191.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
| 178.0ms | (- (pow (sqrt (- (* b b) (* (* 4 a) c))) 3) (pow b 3)) |
| 111.0ms | (- b (* 2 (/ (* a c) b))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 5.5b
Found 4 expressions with local error:
| 21.7b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 13.5b | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 1.3b | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 0.4b | (- (* b b) (* (* 4 a) c)) |
| 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 |
| 11.0ms | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 8.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 2.0ms | (- (* b b) (* (* 4 a) c)) |
| 2.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 145.0ms | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 57.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 48.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 28.0ms | (- (* b b) (* (* 4 a) c)) |
| 826.0ms | (/ (* (* (* (* -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)))) |
| 798.0ms | (/ (* (* (* (* -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)))) |
| 370.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 294.0ms | (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (- (* b b) (* (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 159.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 5.5b
Found 4 expressions with local error:
| 21.7b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 12.1b | (/ (* (* -4 c) a) (* 2 a)) |
| 1.3b | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 0.4b | (- (* b b) (* (* 4 a) c)) |
| 6× | add-log-exp |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | add-sqr-sqrt |
| 5× | add-cube-cbrt |
| 5× | *-un-lft-identity |
| 4× | pow1 |
| 3× | sqrt-prod |
| 2× | flip-- |
| 2× | sqrt-div |
| 2× | flip3-- |
| 1× | associate-/l* |
| 1× | div-inv |
| 1× | flip-+ |
| 1× | div-exp |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | times-frac |
| 1× | flip3-+ |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 8.0ms | (/ (* (* -4 c) a) (* 2 a)) |
| 7.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 2.0ms | (- (* b b) (* (* 4 a) c)) |
| 62.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 44.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 25.0ms | (- (* b b) (* (* 4 a) c)) |
| 6.0ms | (/ (* (* -4 c) a) (* 2 a)) |
| 672.0ms | (- (log (* (* -4 c) a)) (log (* 2 a))) |
| 414.0ms | (/ (* (* (* (* -4 c) a) (* (* -4 c) a)) (* (* -4 c) a)) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
| 357.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 158.0ms | (- (* 2 b) (* 2 (/ (* a c) b))) |
| 128.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 5.5b
Found 4 expressions with local error:
| 21.7b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.3b | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 0.4b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) |
| 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× | pow1 |
| 3× | sqrt-prod |
| 3× | times-frac |
| 3× | associate-/r* |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | associate-/r/ |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | flip3-+ |
| 1× | associate-/l* |
| 1× | div-inv |
| 1× | div-exp |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 21.0ms | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) |
| 15.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 129.0ms | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) |
| 62.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 47.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 38.0ms | (- (* b b) (* (* 4 a) c)) |
| 421.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 276.0ms | (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (- (* b b) (* (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 133.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
| 100.0ms | (- (* 2 b) (* 2 (/ (* a c) b))) |
| 93.0ms | (/ (* (* (* -2 c) (* -2 c)) (* -2 c)) (* (* (+ (sqrt (- (* b b) (* (* 4 a) c))) b) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 5.5b
92.2% (1.9b remaining)
Error of 7.4b against oracle of 5.5b and baseline of 29.4b
| 1.3s | 1063× | body | 2560 | valid |
| 986.0ms | 5847× | body | 80 | valid |
| 652.0ms | 924× | body | 1280 | valid |
| 609.0ms | 435× | body | 640 | valid |
| 440.0ms | 2979× | body | 80 | nan |
| 95.0ms | 102× | body | 5120 | valid |
| 80.0ms | 222× | body | 320 | valid |
| 32.0ms | 145× | body | 160 | valid |