



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.6 |
|---|---|
| Target | 20.6 |
| Herbie | 10.0 |
if b < -1.584464648766991e+89Initial program 58.8
Taylor expanded around -inf 3.0
Simplified3.0
if -1.584464648766991e+89 < b < -1.6242069554645545e-96Initial program 41.9
rmApplied flip--42.0
Applied associate-/l/44.7
Simplified18.7
if -1.6242069554645545e-96 < b < 8.095749828230367e+95Initial program 12.6
rmApplied *-un-lft-identity12.6
Applied associate-/l*12.7
if 8.095749828230367e+95 < b Initial program 42.9
rmApplied *-un-lft-identity42.9
Applied associate-/l*43.0
Taylor expanded around 0 3.9
Simplified3.9
Final simplification10.0
herbie shell --seed 2019005
(FPCore (a b c)
:name "quadm (p42, negative)"
:herbie-target
(if (< b 0) (/ c (* 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))))) (* 2 a)))
Time bar (total: 31.3s)Debug log
1 calls. Slowest were:
| 73.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.4b
Found 4 expressions with local error:
| 24.3b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 1.1b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.4b | (- (* b b) (* 4 (* a c))) |
| 0.1b | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
4 calls. Slowest were:
| 32.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 11.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 4.0ms | (- (* b b) (* 4 (* a c))) |
4 calls. Slowest were:
| 89.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 55.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 50.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 33.0ms | (- (* b b) (* 4 (* a c))) |
46 calls. Slowest were:
| 348.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 251.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 239.0ms | (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 6.4b
Found 4 expressions with local error:
| 24.3b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 16.9b | (/ (* (* c 4) a) (* (* 2 a) (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 1.1b | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.4b | (- (* b b) (* 4 (* a c))) |
4 calls. Slowest were:
| 11.0ms | (/ (* (* c 4) a) (* (* 2 a) (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 6.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 2.0ms | (- (* b b) (* 4 (* a c))) |
4 calls. Slowest were:
| 157.0ms | (/ (* (* c 4) a) (* (* 2 a) (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 76.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 64.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 36.0ms | (- (* b b) (* 4 (* a c))) |
31 calls. Slowest were:
| 770.0ms | (/ (* (* (* (* c 4) a) (* (* c 4) a)) (* (* c 4) 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)))))))) |
| 756.0ms | (/ (* (* (* (* c 4) a) (* (* c 4) a)) (* (* c 4) 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)))))))) |
| 313.0ms | (+ (* (- b) (- b)) (- (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c))))) (* (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 6.4b
Found 4 expressions with local error:
| 25.1b | (cbrt (- (* b b) (* 4 (* a c)))) |
| 25.1b | (cbrt (- (* b b) (* 4 (* a c)))) |
| 25.1b | (cbrt (- (* b b) (* 4 (* a c)))) |
| 24.3b | (sqrt (* (* (cbrt (- (* b b) (* 4 (* a c)))) (cbrt (- (* b b) (* 4 (* a c))))) (cbrt (- (* b b) (* 4 (* a c)))))) |
4 calls. Slowest were:
| 11.0ms | (sqrt (* (* (cbrt (- (* b b) (* 4 (* a c)))) (cbrt (- (* b b) (* 4 (* a c))))) (cbrt (- (* b b) (* 4 (* a c)))))) |
| 3.0ms | (cbrt (- (* b b) (* 4 (* a c)))) |
| 3.0ms | (cbrt (- (* b b) (* 4 (* a c)))) |
4 calls. Slowest were:
| 88.0ms | (cbrt (- (* b b) (* 4 (* a c)))) |
| 81.0ms | (cbrt (- (* b b) (* 4 (* a c)))) |
| 77.0ms | (cbrt (- (* b b) (* 4 (* a c)))) |
| 59.0ms | (sqrt (* (* (cbrt (- (* b b) (* 4 (* a c)))) (cbrt (- (* b b) (* 4 (* a c))))) (cbrt (- (* b b) (* 4 (* a c)))))) |
54 calls. Slowest were:
| 896.0ms | (sqrt (* (* (cbrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) (cbrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3)))) (cbrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))))) |
| 574.0ms | (sqrt (* (* (cbrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) (cbrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3)))) (cbrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))))) |
| 568.0ms | (sqrt (* (* (cbrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) (cbrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3)))) (cbrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))))) |
8 alts after pruning (6 fresh and 2 done)
Merged error: 6.4b
Found 4 expressions with local error:
| 24.3b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 1.1b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.4b | (- (* b b) (* 4 (* a c))) |
| 0.2b | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
4 calls. Slowest were:
| 6.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 6.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 2.0ms | (- (* b b) (* 4 (* a c))) |
4 calls. Slowest were:
| 94.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 75.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 54.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 42.0ms | (- (* b b) (* 4 (* a c))) |
39 calls. Slowest were:
| 426.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c))))) (* (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 318.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 147.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) |
8 alts after pruning (5 fresh and 3 done)
Merged error: 6.4b
86% (4.1b remaining)
Error of 10.0b against oracle of 6.0b and baseline of 35.0b