


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
Results
if b < -1.2165750385518375e+144Initial program 57.3
Taylor expanded around -inf 11.5
if -1.2165750385518375e+144 < b < 3.0061788106952535e-55Initial program 13.1
Simplified13.1
rmApplied add-cube-cbrt13.4
rmApplied add-cube-cbrt13.5
rmApplied add-sqr-sqrt13.5
Applied rem-sqrt-square13.5
Simplified13.1
if 3.0061788106952535e-55 < b Initial program 53.4
Taylor expanded around inf 18.3
Final simplification14.8
herbie shell --seed 2019022
(FPCore (a b c d)
:name "Cubic critical"
(/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)))
Time bar (total: 15.0s)Debug log
| 37.0ms | 39× | body | 2560 | valid |
| 22.0ms | 181× | body | 80 | valid |
| 14.0ms | 23× | body | 1280 | valid |
| 12.0ms | 105× | body | 80 | nan |
| 9.0ms | 20× | body | 640 | valid |
| 4.0ms | 7× | body | 320 | valid |
| 3.0ms | 4× | body | 5120 | valid |
| 1.0ms | 3× | body | 160 | valid |
| 61.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 31.8b
Found 4 expressions with local error:
| 21.8b | (sqrt (- (* b b) (* (* 3 a) c))) |
| 1.0b | (- (* b b) (* (* 3 a) c)) |
| 1.0b | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 0.2b | (* (* 3 a) c) |
| 8× | add-log-exp |
| 7× | *-un-lft-identity |
| 6× | add-sqr-sqrt |
| 5× | add-cube-cbrt |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | sqrt-prod |
| 2× | flip-- |
| 2× | sqrt-div |
| 2× | flip3-- |
| 1× | distribute-lft-out |
| 1× | flip-+ |
| 1× | associate-+l- |
| 1× | *-commutative |
| 1× | associate-*l* |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | flip3-+ |
| 1× | sub-neg |
| 1× | sum-log |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | neg-sub0 |
| 10.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 4.0ms | (* (* 3 a) c) |
| 4.0ms | (- (* b b) (* (* 3 a) c)) |
| 3.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 59.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 55.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 39.0ms | (- (* b b) (* (* 3 a) c)) |
| 11.0ms | (* (* 3 a) c) |
| 401.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 3 a) c) 3))) |
| 150.0ms | (- b (* 3/2 (/ (* a c) b))) |
| 131.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 3 a) c) (* (* 3 a) c)))) |
| 85.0ms | (* (exp (- b)) (exp (sqrt (- (* b b) (* (* 3 a) c))))) |
| 76.0ms | (/ (exp (* b b)) (exp (* (* 3 a) c))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 12.6b
Found 4 expressions with local error:
| 21.8b | (sqrt (+ (* -3 (* c a)) (* b b))) |
| 1.0b | (+ (* -3 (* c a)) (* b b)) |
| 1.0b | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 0.2b | (* -3 (* c a)) |
| 6× | add-log-exp |
| 6× | add-sqr-sqrt |
| 5× | add-cube-cbrt |
| 5× | *-un-lft-identity |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | sqrt-prod |
| 2× | flip-+ |
| 2× | sqrt-div |
| 2× | flip3-+ |
| 1× | flip-- |
| 1× | *-commutative |
| 1× | associate-*r* |
| 1× | pow1/2 |
| 1× | flip3-- |
| 1× | sub-neg |
| 1× | sum-log |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 10.0ms | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 3.0ms | (+ (* -3 (* c a)) (* b b)) |
| 2.0ms | (sqrt (+ (* -3 (* c a)) (* b b))) |
| 1.0ms | (* -3 (* c a)) |
| 70.0ms | (sqrt (+ (* -3 (* c a)) (* b b))) |
| 61.0ms | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 37.0ms | (+ (* -3 (* c a)) (* b b)) |
| 16.0ms | (* -3 (* c a)) |
| 269.0ms | (sqrt (- (* (* -3 (* c a)) (* -3 (* c a))) (* (* b b) (* b b)))) |
| 162.0ms | (sqrt (+ (pow (* -3 (* c a)) 3) (pow (* b b) 3))) |
| 38.0ms | (- (pow b 2) (* 3 (* a c))) |
| 34.0ms | (- (pow b 2) (* 3 (* a c))) |
| 22.0ms | (- (pow b 2) (* 3 (* a c))) |
14 alts after pruning (13 fresh and 1 done)
Merged error: 12.6b
Found 4 expressions with local error:
| 21.8b | (sqrt (+ (* (* (cbrt (* -3 (* c a))) (cbrt (* -3 (* c a)))) (cbrt (* -3 (* c a)))) (* b b))) |
| 16.4b | (cbrt (* -3 (* c a))) |
| 16.4b | (cbrt (* -3 (* c a))) |
| 16.4b | (cbrt (* -3 (* c a))) |
| 6× | add-sqr-sqrt |
| 5× | add-cube-cbrt |
| 5× | *-un-lft-identity |
| 4× | add-log-exp |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | sqrt-prod |
| 3× | pow1/3 |
| 3× | cbrt-prod |
| 2× | sqrt-div |
| 1× | flip-+ |
| 1× | pow1/2 |
| 1× | flip3-+ |
| 1× | rem-sqrt-square |
| 19.0ms | (sqrt (+ (* (* (cbrt (* -3 (* c a))) (cbrt (* -3 (* c a)))) (cbrt (* -3 (* c a)))) (* b b))) |
| 1.0ms | (cbrt (* -3 (* c a))) |
| 1.0ms | (cbrt (* -3 (* c a))) |
| 1.0ms | (cbrt (* -3 (* c a))) |
| 126.0ms | (sqrt (+ (* (* (cbrt (* -3 (* c a))) (cbrt (* -3 (* c a)))) (cbrt (* -3 (* c a)))) (* b b))) |
| 124.0ms | (cbrt (* -3 (* c a))) |
| 120.0ms | (cbrt (* -3 (* c a))) |
| 110.0ms | (cbrt (* -3 (* c a))) |
| 344.0ms | (sqrt (+ (pow (* (* (cbrt (* -3 (* c a))) (cbrt (* -3 (* c a)))) (cbrt (* -3 (* c a)))) 3) (pow (* b b) 3))) |
| 156.0ms | (* (exp (* -1/3 (+ (log (/ 1 a)) (log (/ 1 c))))) (cbrt -3)) |
| 141.0ms | (sqrt (- (* (* (* (cbrt (* -3 (* c a))) (cbrt (* -3 (* c a)))) (cbrt (* -3 (* c a)))) (* (* (cbrt (* -3 (* c a))) (cbrt (* -3 (* c a)))) (cbrt (* -3 (* c a))))) (* (* b b) (* b b)))) |
| 128.0ms | (* (exp (* -1/3 (+ (log (/ -1 c)) (log (/ -1 a))))) (cbrt -3)) |
| 113.0ms | (* (exp (* -1/3 (+ (log (/ -1 c)) (log (/ -1 a))))) (cbrt -3)) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 12.6b
Found 4 expressions with local error:
| 21.8b | (sqrt (+ (* (* (cbrt (* -3 (* c a))) (* (* (cbrt (cbrt (* -3 (* c a)))) (cbrt (cbrt (* -3 (* c a))))) (cbrt (cbrt (* -3 (* c a)))))) (cbrt (* -3 (* c a)))) (* b b))) |
| 16.4b | (cbrt (* -3 (* c a))) |
| 16.4b | (cbrt (* -3 (* c a))) |
| 16.4b | (cbrt (* -3 (* c a))) |
| 6× | add-sqr-sqrt |
| 5× | add-cube-cbrt |
| 5× | *-un-lft-identity |
| 4× | add-log-exp |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | sqrt-prod |
| 3× | pow1/3 |
| 3× | cbrt-prod |
| 2× | sqrt-div |
| 1× | flip-+ |
| 1× | pow1/2 |
| 1× | flip3-+ |
| 1× | rem-sqrt-square |
| 128.0ms | (sqrt (+ (* (* (cbrt (* -3 (* c a))) (* (* (cbrt (cbrt (* -3 (* c a)))) (cbrt (cbrt (* -3 (* c a))))) (cbrt (cbrt (* -3 (* c a)))))) (cbrt (* -3 (* c a)))) (* b b))) |
| 1.0ms | (cbrt (* -3 (* c a))) |
| 1.0ms | (cbrt (* -3 (* c a))) |
| 1.0ms | (cbrt (* -3 (* c a))) |
| 139.0ms | (cbrt (* -3 (* c a))) |
| 124.0ms | (sqrt (+ (* (* (cbrt (* -3 (* c a))) (* (* (cbrt (cbrt (* -3 (* c a)))) (cbrt (cbrt (* -3 (* c a))))) (cbrt (cbrt (* -3 (* c a)))))) (cbrt (* -3 (* c a)))) (* b b))) |
| 119.0ms | (cbrt (* -3 (* c a))) |
| 116.0ms | (cbrt (* -3 (* c a))) |
| 547.0ms | (sqrt (- (* (* (* (cbrt (* -3 (* c a))) (* (* (cbrt (cbrt (* -3 (* c a)))) (cbrt (cbrt (* -3 (* c a))))) (cbrt (cbrt (* -3 (* c a)))))) (cbrt (* -3 (* c a)))) (* (* (cbrt (* -3 (* c a))) (* (* (cbrt (cbrt (* -3 (* c a)))) (cbrt (cbrt (* -3 (* c a))))) (cbrt (cbrt (* -3 (* c a)))))) (cbrt (* -3 (* c a))))) (* (* b b) (* b b)))) |
| 331.0ms | (sqrt (+ (pow (* (* (cbrt (* -3 (* c a))) (* (* (cbrt (cbrt (* -3 (* c a)))) (cbrt (cbrt (* -3 (* c a))))) (cbrt (cbrt (* -3 (* c a)))))) (cbrt (* -3 (* c a)))) 3) (pow (* b b) 3))) |
| 166.0ms | (* (exp (* -1/3 (+ (log (/ 1 a)) (log (/ 1 c))))) (cbrt -3)) |
| 164.0ms | (* (exp (* -1/3 (+ (log (/ 1 a)) (log (/ 1 c))))) (cbrt -3)) |
| 159.0ms | (sqrt (sqrt (+ (* (* (cbrt (* -3 (* c a))) (* (* (cbrt (cbrt (* -3 (* c a)))) (cbrt (cbrt (* -3 (* c a))))) (cbrt (cbrt (* -3 (* c a)))))) (cbrt (* -3 (* c a)))) (* b b)))) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 12.6b
86.8% (2.8b remaining)
Error of 14.8b against oracle of 12.0b and baseline of 33.2b
| 1.2s | 5812× | body | 80 | valid |
| 874.0ms | 1043× | body | 2560 | valid |
| 719.0ms | 968× | body | 1280 | valid |
| 464.0ms | 2937× | body | 80 | nan |
| 204.0ms | 452× | body | 640 | valid |
| 125.0ms | 104× | body | 5120 | valid |
| 107.0ms | 201× | body | 320 | valid |
| 28.0ms | 125× | body | 160 | valid |