



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.1 |
|---|---|
| Target | 20.7 |
| Herbie | 10.3 |
if b < -2.8671071746705253e+127Initial program 51.2
Simplified51.2
Taylor expanded around -inf 2.8
if -2.8671071746705253e+127 < b < 1.354454736393541e-170Initial program 11.3
Simplified11.4
rmApplied clear-num11.5
if 1.354454736393541e-170 < b < 9.22274670563189e-75Initial program 24.9
Simplified24.9
rmApplied flip--25.0
Applied associate-/l/31.2
Simplified25.9
rmApplied times-frac18.2
Simplified18.2
if 9.22274670563189e-75 < b < 1.849572225033813e-22Initial program 37.0
Simplified37.0
rmApplied add-sqr-sqrt37.5
Applied times-frac37.5
if 1.849572225033813e-22 < b Initial program 54.3
Simplified54.3
Taylor expanded around inf 6.8
Simplified6.8
Final simplification10.3
herbie shell --seed 2019005
(FPCore (a b c)
:name "quadp (p42, positive)"
:herbie-target
(if (< b 0) (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) (/ c (* a (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))))
(/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
Time bar (total: 20.0s)Debug log
1 calls. Slowest were:
| 70.0ms | (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 32.2b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.1b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 0.4b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
4 calls. Slowest were:
| 40.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 21.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 100.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 77.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 48.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 33.0ms | (- (* b b) (* (* 4 a) c)) |
32 calls. Slowest were:
| 377.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 301.0ms | (- (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (* b b)) |
| 165.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 6.5b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.1b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 0.4b | (- (* b b) (* (* 4 a) c)) |
| 0.2b | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* 4 a) c))) b))) |
4 calls. Slowest were:
| 10.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 6.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 2.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 93.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 74.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 54.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 25.0ms | (- (* b b) (* (* 4 a) c)) |
30 calls. Slowest were:
| 371.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 187.0ms | (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (+ (* b b) (* (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 114.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 6.4b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.1b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 1.1b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
4 calls. Slowest were:
| 19.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 13.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 67.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 58.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 47.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 41.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
24 calls. Slowest were:
| 402.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 334.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 137.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 6.4b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 16.3b | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 1.1b | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 0.4b | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 12.0ms | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 8.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 164.0ms | (/ (* (* -4 c) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 73.0ms | (+ (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 72.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 43.0ms | (- (* b b) (* (* 4 a) c)) |
29 calls. Slowest were:
| 941.0ms | (/ (* (* (* (* -4 c) a) (* (* -4 c) a)) (* (* -4 c) a)) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (+ (sqrt (- (* b b) (* (* 4 a) c))) b) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)))) |
| 814.0ms | (/ (* (* (* (* -4 c) a) (* (* -4 c) a)) (* (* -4 c) a)) (* (* (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b))) (* (* 2 a) (+ (sqrt (- (* b b) (* (* 4 a) c))) b)))) |
| 442.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 6.4b
84.2% (4.3b remaining)
Error of 10.3b against oracle of 6.0b and baseline of 33.2b