


Bits error versus a



Bits error versus b



Bits error versus c
Results
if b < -8.422430181155895e+146Initial program 57.7
Initial simplification57.7
Taylor expanded around -inf 2.3
if -8.422430181155895e+146 < b < -1.347277676039659e-308Initial program 8.4
Initial simplification8.4
Taylor expanded around -inf 8.4
if -1.347277676039659e-308 < b < 4.331216038455128e+110Initial program 33.8
Initial simplification33.8
rmApplied flip--33.9
Applied associate-/l/38.2
Simplified21.5
rmApplied associate-/r*16.1
Taylor expanded around 0 9.0
if 4.331216038455128e+110 < b Initial program 59.2
Initial simplification59.2
Taylor expanded around inf 1.8
Simplified1.8
Final simplification6.3
herbie shell --seed 2018360
(FPCore (a b c)
:name "Quadratic roots, full range"
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 15.6s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.6b
Found 4 expressions with local error:
| 22.5b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.9b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.4b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
4 calls. Slowest were:
| 41.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 20.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 77.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 74.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 64.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 32.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 240.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 173.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 172.0ms | (- (pow (sqrt (- (* b b) (* (* c a) 4))) 3) (pow b 3)) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.5b | (sqrt (- (* b b) (* (* c a) 4))) |
| 13.3b | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 0.8b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.4b | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 24.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 15.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 173.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 57.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 57.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 37.0ms | (- (* b b) (* (* c a) 4)) |
29 calls. Slowest were:
| 820.0ms | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)))) |
| 780.0ms | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (+ (sqrt (- (* b b) (* (* c a) 4))) b) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) (+ (sqrt (- (* b b) (* (* c a) 4))) b)))) |
| 505.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 4.7b
Found 4 expressions with local error:
| 22.5b | (sqrt (- (* b b) (* (* c a) 4))) |
| 12.4b | (/ (* (* c -4) a) (* 2 a)) |
| 0.8b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.4b | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 8.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 5.0ms | (/ (* (* c -4) a) (* 2 a)) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 60.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 46.0ms | (- (* b b) (* (* c a) 4)) |
| 43.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 8.0ms | (/ (* (* c -4) a) (* 2 a)) |
24 calls. Slowest were:
| 553.0ms | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
| 352.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 338.0ms | (- (log (* (* c -4) a)) (log (* 2 a))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 4.7b
Found 4 expressions with local error:
| 22.5b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.8b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.4b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
4 calls. Slowest were:
| 22.0ms | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 15.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 152.0ms | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 58.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 57.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 36.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 313.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 199.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (- (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 142.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 4.7b
96.6% (0.8b remaining)
Error of 6.3b against oracle of 5.5b and baseline of 29.8b