


Bits error versus a



Bits error versus b



Bits error versus c
Initial program 52.3
rmApplied flip-+52.3
Applied associate-/l/52.3
Simplified0.4
rmApplied associate-/r*0.2
Simplified0.2
Taylor expanded around 0 0.1
Final simplification0.1
herbie shell --seed 2019004 +o rules:numerics
(FPCore (a b c)
:name "Quadratic roots, wide range"
:pre (and (< 4.930380657631324e-32 a 2.028240960365167e+31) (< 4.930380657631324e-32 b 2.028240960365167e+31) (< 4.930380657631324e-32 c 2.028240960365167e+31))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 1.2m)Debug log
1 calls. Slowest were:
| 100.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 51.6b
Found 4 expressions with local error:
| 27.2b | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.1b | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 0.1b | (- (* b b) (* (* 4 a) c)) |
| 0.0b | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 27.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 16.0ms | (- (* b b) (* (* 4 a) c)) |
| 7.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
4 calls. Slowest were:
| 78.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 74.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 68.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 29.0ms | (- (* b b) (* (* 4 a) c)) |
40 calls. Slowest were:
| 429.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 376.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 211.0ms | (+ (pow (- b) 3) (pow (sqrt (- (* b b) (* (* 4 a) c))) 3)) |
3 alts after pruning (3 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 0.4b | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 0.1b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.1b | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
4 calls. Slowest were:
| 27.0ms | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 10.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 7.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
4 calls. Slowest were:
| 151.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 96.0ms | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 56.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 28.0ms | (- (* b b) (* (* 4 a) c)) |
143 calls. Slowest were:
| 1.1s | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))))) |
| 821.0ms | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 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))))))) |
| 595.0ms | (* (fma (- (cbrt (sqrt (- (* b b) (* (* 4 a) c))))) (* (cbrt (sqrt (- (* b b) (* (* 4 a) c)))) (cbrt (sqrt (- (* b b) (* (* 4 a) c))))) (* (cbrt (sqrt (- (* b b) (* (* 4 a) c)))) (* (cbrt (sqrt (- (* b b) (* (* 4 a) c)))) (cbrt (sqrt (- (* b b) (* (* 4 a) c))))))) (* 2 a)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.1b | (/ (* 4 (* c a)) (* 2 a)) |
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.1b | (/ (/ (* 4 (* c a)) (* 2 a)) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
4 calls. Slowest were:
| 17.0ms | (/ (/ (* 4 (* c a)) (* 2 a)) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 8.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 3.0ms | (/ (* 4 (* c a)) (* 2 a)) |
4 calls. Slowest were:
| 174.0ms | (/ (/ (* 4 (* c a)) (* 2 a)) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 88.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 52.0ms | (fma (* a -4) c (* b b)) |
| 7.0ms | (/ (* 4 (* c a)) (* 2 a)) |
149 calls. Slowest were:
| 1.2s | (/ (* (cbrt (/ (* 4 (* c a)) (* 2 a))) (cbrt (/ (* 4 (* c a)) (* 2 a)))) (+ (sqrt (- b)) (sqrt (sqrt (fma (* a -4) c (* b b)))))) |
| 1.1s | (/ (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* 2 a) (* 2 a)) (* 2 a))) (* (* (- (- b) (sqrt (fma (* a -4) c (* b b)))) (- (- b) (sqrt (fma (* a -4) c (* b b))))) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 1.1s | (/ (sqrt (/ (* 4 (* c a)) (* 2 a))) (sqrt (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
5 alts after pruning (5 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.1b | (/ (* 2 c) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 0.0b | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 7.0ms | (/ (* 2 c) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 5.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 1.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 142.0ms | (/ (* 2 c) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 53.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 39.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 29.0ms | (fma (* a -4) c (* b b)) |
94 calls. Slowest were:
| 481.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 73.0ms | (fma 1 (- b) (- (* (sqrt (cbrt (fma (* a -4) c (* b b)))) (sqrt (* (cbrt (fma (* a -4) c (* b b))) (cbrt (fma (* a -4) c (* b b)))))))) |
| 66.0ms | (/ (* (* (* 2 c) (* 2 c)) (* 2 c)) (* (* (- (- b) (sqrt (fma (* a -4) c (* b b)))) (- (- b) (sqrt (fma (* a -4) c (* b b))))) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
5 alts after pruning (4 fresh and 1 done)
Merged error: 0b
0% (0.1b remaining)
Error of 0.1b against oracle of 0.0b and baseline of 0.1b