


Bits error versus a



Bits error versus b



Bits error versus c
Results
if b < -4.4165345932310753e+135Initial program 54.4
Taylor expanded around -inf 10.4
if -4.4165345932310753e+135 < b < 7.564155873244899e+99Initial program 9.0
rmApplied add-cube-cbrt9.2
Applied sqrt-prod9.2
Simplified9.2
rmApplied add-cube-cbrt9.2
Applied cbrt-prod9.3
if 7.564155873244899e+99 < b Initial program 30.9
Taylor expanded around inf 7.1
Final simplification8.9
herbie shell --seed 2019007
(FPCore (a b c)
:name "jeff quadratic root 2"
(if (>= b 0) (/ (* 2 c) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a))))
Time bar (total: 24.0s)Debug log
1 calls. Slowest were:
| 107.0ms | (if (>= b 0) (/ (* 2 c) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 22.1b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.7b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 1.5b | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
4 calls. Slowest were:
| 7.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 6.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 5.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
4 calls. Slowest were:
| 77.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 74.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 53.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 46.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
32 calls. Slowest were:
| 419.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 371.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 153.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 9.4b
Found 4 expressions with local error:
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 9.0b | (+ (- b) (* (fabs (cbrt (- (* b b) (* (* 4 a) c)))) (sqrt (cbrt (- (* b b) (* (* 4 a) c)))))) |
4 calls. Slowest were:
| 8.0ms | (+ (- b) (* (fabs (cbrt (- (* b b) (* (* 4 a) c)))) (sqrt (cbrt (- (* b b) (* (* 4 a) c)))))) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 3.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 115.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 108.0ms | (+ (- b) (* (fabs (cbrt (- (* b b) (* (* 4 a) c)))) (sqrt (cbrt (- (* b b) (* (* 4 a) c)))))) |
| 100.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 66.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
29 calls. Slowest were:
| 489.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 378.0ms | (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 291.0ms | (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 9.4b
Found 4 expressions with local error:
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 1.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 1.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 1.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 130.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 109.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 96.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 96.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
32 calls. Slowest were:
| 476.0ms | (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 371.0ms | (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 362.0ms | (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 9.4b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 68.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 67.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 64.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 55.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
36 calls. Slowest were:
| 698.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 441.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 406.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 9.4b
91.8% (1.0b remaining)
Error of 8.9b against oracle of 7.9b and baseline of 20.5b