


Bits error versus a



Bits error versus b



Bits error versus c
Results
if b < -2.835340391058544e+148Initial program 58.7
Simplified58.7
Taylor expanded around -inf 2.1
if -2.835340391058544e+148 < b < 3.587545867948427e-80Initial program 12.5
Simplified12.5
if 3.587545867948427e-80 < b Initial program 52.6
Simplified52.6
Taylor expanded around inf 8.5
Simplified8.5
Final simplification9.7
herbie shell --seed 2019004
(FPCore (a b c)
:name "Quadratic roots, full range"
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 19.1s)Debug log
1 calls. Slowest were:
| 96.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 32.6b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
| 0.3b | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 0.0b | (* (* c a) 4) |
4 calls. Slowest were:
| 21.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 4.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 98.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 51.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 28.0ms | (- (* b b) (* (* c a) 4)) |
| 18.0ms | (* (* c a) 4) |
32 calls. Slowest were:
| 420.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 320.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 193.0ms | (- (pow (sqrt (- (* b b) (* (* c a) 4))) 3) (pow b 3)) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
| 0.4b | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 0.2b | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
4 calls. Slowest were:
| 23.0ms | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 12.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 4.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 168.0ms | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 110.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 67.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 45.0ms | (- (* b b) (* (* c a) 4)) |
43 calls. Slowest were:
| 361.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 171.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
| 156.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))))) |
| 16.3b | (cbrt (* (* c a) 4)) |
| 16.3b | (cbrt (* (* c a) 4)) |
| 16.3b | (cbrt (* (* c a) 4)) |
4 calls. Slowest were:
| 4.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))))) |
| 3.0ms | (cbrt (* (* c a) 4)) |
| 3.0ms | (cbrt (* (* c a) 4)) |
4 calls. Slowest were:
| 141.0ms | (cbrt (* (* c a) 4)) |
| 136.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))))) |
| 133.0ms | (cbrt (* (* c a) 4)) |
| 112.0ms | (cbrt (* (* c a) 4)) |
21 calls. Slowest were:
| 353.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))) 3))) |
| 206.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4)))))) |
| 182.0ms | (* (exp (* -1/3 (+ (log (/ 1 a)) (log (/ 1 c))))) (cbrt 4)) |
9 alts after pruning (6 fresh and 3 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 5.0ms | (- (* b b) (* (* c a) 4)) |
| 2.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 61.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 48.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 40.0ms | (- (* b b) (* (* c a) 4)) |
| 31.0ms | (- (* b b) (* (* c a) 4)) |
26 calls. Slowest were:
| 412.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 393.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 128.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
9 alts after pruning (5 fresh and 4 done)
Merged error: 4.8b
85.1% (4.2b remaining)
Error of 9.7b against oracle of 5.6b and baseline of 33.5b