



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.3 |
|---|---|
| Target | 20.7 |
| Herbie | 6.7 |
if b < -1.0041518941750358e+133Initial program 60.5
Taylor expanded around -inf 1.9
Simplified1.9
if -1.0041518941750358e+133 < b < 3.9279776760338723e-290Initial program 32.8
rmApplied div-inv32.8
rmApplied flip--32.9
Applied associate-*l/33.0
Simplified14.9
Taylor expanded around 0 8.9
if 3.9279776760338723e-290 < b < 2.653900512555386e+58Initial program 8.9
Taylor expanded around -inf 8.9
if 2.653900512555386e+58 < b Initial program 37.0
Taylor expanded around -inf 37.0
Taylor expanded around inf 4.8
Final simplification6.7
herbie shell --seed 2019010
(FPCore (a b c)
:name "quadm (p42, negative)"
:herbie-target
(if (< b 0) (/ c (* a (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))) (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
(/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
Time bar (total: 23.1s)Debug log
| 85.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:
| 23.4b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 0.3b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.2b | (- (* b b) (* 4 (* a c))) |
| 0.1b | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 20× | add-sqr-sqrt |
| 13× | *-un-lft-identity |
| 8× | add-log-exp |
| 7× | add-cube-cbrt |
| 6× | difference-of-squares |
| 6× | associate-/l* |
| 6× | sqrt-prod |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | times-frac |
| 4× | flip-- |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | distribute-lft-out-- |
| 2× | associate-/l/ |
| 2× | diff-log |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 30.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 11.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 4.0ms | (- (* b b) (* 4 (* a c))) |
| 4.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 81.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 67.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 63.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 33.0ms | (- (* b b) (* 4 (* a c))) |
| 319.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 292.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 201.0ms | (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) |
| 174.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 133.0ms | (/ (* (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (- (- b) (sqrt (- (* b b) (* 4 (* a c)))))) (- (- b) (sqrt (- (* b b) (* 4 (* a c)))))) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 6.2b
Found 4 expressions with local error:
| 23.4b | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 0.3b | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 0.2b | (- (pow b 2) (* 4 (* a c))) |
| 0.1b | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 20× | add-sqr-sqrt |
| 13× | *-un-lft-identity |
| 8× | add-log-exp |
| 7× | add-cube-cbrt |
| 6× | difference-of-squares |
| 6× | associate-/l* |
| 6× | sqrt-prod |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | times-frac |
| 4× | flip-- |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | distribute-lft-out-- |
| 2× | associate-/l/ |
| 2× | diff-log |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 36.0ms | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 13.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 6.0ms | (- (pow b 2) (* 4 (* a c))) |
| 5.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 88.0ms | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 71.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 70.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 28.0ms | (- (pow b 2) (* 4 (* a c))) |
| 843.0ms | (- (* (- b) (- b)) (* (sqrt (- (pow b 2) (* 4 (* a c)))) (sqrt (- (pow b 2) (* 4 (* a c)))))) |
| 401.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 180.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 170.0ms | (- (pow (- b) 3) (pow (sqrt (- (pow b 2) (* 4 (* a c)))) 3)) |
| 125.0ms | (- b (* 2 (/ (* a c) b))) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 6.2b
Found 4 expressions with local error:
| 23.4b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 0.3b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.2b | (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (/ 1 (* 2 a))) |
| 0.2b | (- (* b b) (* 4 (* a c))) |
| 16× | add-sqr-sqrt |
| 11× | *-un-lft-identity |
| 8× | add-log-exp |
| 8× | add-exp-log |
| 7× | add-cube-cbrt |
| 6× | associate-*l* |
| 6× | add-cbrt-cube |
| 6× | pow1 |
| 5× | flip-- |
| 5× | sqrt-prod |
| 5× | flip3-- |
| 4× | difference-of-squares |
| 4× | associate-*r* |
| 2× | distribute-lft-out-- |
| 2× | prod-exp |
| 2× | associate-*l/ |
| 2× | diff-log |
| 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× | pow-prod-down |
| 1× | pow1/2 |
| 1× | rec-exp |
| 1× | rem-sqrt-square |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 19.0ms | (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (/ 1 (* 2 a))) |
| 5.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 5.0ms | (- (* b b) (* 4 (* a c))) |
| 2.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 90.0ms | (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (/ 1 (* 2 a))) |
| 77.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 48.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 27.0ms | (- (* b b) (* 4 (* a c))) |
| 958.0ms | (* (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) 1) |
| 436.0ms | (* (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) (/ 1 (* 2 a))) |
| 414.0ms | (* (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) 1) |
| 414.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 337.0ms | (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (/ 1 (* 2 a))) |
8 alts after pruning (6 fresh and 2 done)
Merged error: 6.2b
Found 4 expressions with local error:
| 23.4b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 15.3b | (* (* (* 4 a) c) (/ 1/2 a)) |
| 1.7b | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.2b | (- (* b b) (* 4 (* a c))) |
| 8× | add-log-exp |
| 8× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 6× | add-cube-cbrt |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | pow1 |
| 4× | associate-*r* |
| 3× | sqrt-prod |
| 2× | flip-- |
| 2× | sqrt-div |
| 2× | flip3-- |
| 1× | distribute-lft-out |
| 1× | div-inv |
| 1× | flip-+ |
| 1× | cbrt-unprod |
| 1× | associate-+l- |
| 1× | *-commutative |
| 1× | associate-*r/ |
| 1× | prod-exp |
| 1× | pow-prod-down |
| 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 |
| 15.0ms | (* (* (* 4 a) c) (/ 1/2 a)) |
| 5.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 4.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 2.0ms | (- (* b b) (* 4 (* a c))) |
| 53.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 48.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 38.0ms | (- (* b b) (* 4 (* a c))) |
| 7.0ms | (* (* (* 4 a) c) (/ 1/2 a)) |
| 368.0ms | (+ (log (* (* 4 a) c)) (log (/ 1/2 a))) |
| 366.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 199.0ms | (* (* (* 4 a) c) (/ 1/2 a)) |
| 189.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 125.0ms | (- b (* 2 (/ (* a c) b))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 6.2b
95.1% (1.2b remaining)
Error of 6.7b against oracle of 5.5b and baseline of 29.6b