


Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 44.0
Simplified44.0
rmApplied flip--44.0
Applied associate-/l/44.0
Simplified0.4
rmApplied associate-/r*0.2
Taylor expanded around 0 0.2
rmApplied flip3--0.2
Applied sqrt-div0.2
Final simplification0.2
herbie shell --seed 2019007
(FPCore (a b c)
:name "Quadratic roots, medium range"
:pre (and (< 1.1102230246251565e-16 a 9007199254740992.0) (< 1.1102230246251565e-16 b 9007199254740992.0) (< 1.1102230246251565e-16 c 9007199254740992.0))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 31.5s)Debug log
1 calls. Slowest were:
| 59.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 45.7b
Found 4 expressions with local error:
| 40.0b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.2b | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 0.1b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 20.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 12.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 3.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 90.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 79.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 40.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 35.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 369.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 242.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 219.0ms | (- (pow (sqrt (- (* b b) (* (* c a) 4))) 3) (pow b 3)) |
5 alts after pruning (5 fresh and 0 done)
Merged error: 0.2b
Found 4 expressions with local error:
| 0.3b | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 0.1b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.1b | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
4 calls. Slowest were:
| 24.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 14.0ms | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 10.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
4 calls. Slowest were:
| 148.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 92.0ms | (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 54.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 27.0ms | (- (* b b) (* (* c a) 4)) |
33 calls. Slowest were:
| 900.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)))) |
| 758.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)))) |
| 320.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (- (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.1b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (/ (* (* c -4) a) (* 2 a)) |
| 0.1b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.1b | (/ (/ (* (* c -4) a) (* 2 a)) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
4 calls. Slowest were:
| 30.0ms | (/ (/ (* (* c -4) a) (* 2 a)) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 16.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 10.0ms | (/ (* (* c -4) a) (* 2 a)) |
4 calls. Slowest were:
| 122.0ms | (/ (/ (* (* c -4) a) (* 2 a)) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 68.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 34.0ms | (- (* b b) (* (* c a) 4)) |
| 5.0ms | (/ (* (* c -4) a) (* 2 a)) |
62 calls. Slowest were:
| 1.1s | (/ (/ (* (* c -4) a) (* 2 a)) (cbrt (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 1.1s | (/ (/ (* (* c -4) a) (* 2 a)) (sqrt (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 1.1s | (/ (/ (* (* (* (* 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))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.1b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.1b | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 0.1b | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 12.0ms | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 8.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 133.0ms | (/ (* -2 c) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 66.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 46.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 37.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 356.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 261.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (- (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 147.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 0b
0% (0.2b remaining)
Error of 0.2b against oracle of 0.0b and baseline of 0.2b