



Bits error versus a




Bits error versus b




Bits error versus c
| Original | 32.8 |
|---|---|
| Target | 20.2 |
| Herbie | 8.8 |
if b < -4.2485149476691936e+130Initial program 60.8
Taylor expanded around inf 60.8
Simplified60.8
Taylor expanded around -inf 2.0
Simplified2.0
if -4.2485149476691936e+130 < b < -9.532534227544595e-103Initial program 42.3
rmApplied flip--42.3
Applied associate-/l/44.3
Simplified16.8
rmApplied associate-/r*13.1
Simplified13.1
if -9.532534227544595e-103 < b < 4.62148638684652e+78Initial program 11.5
Taylor expanded around inf 11.5
Simplified11.5
rmApplied div-sub11.5
if 4.62148638684652e+78 < b Initial program 40.4
Taylor expanded around inf 40.4
Simplified40.4
Taylor expanded around inf 5.1
Final simplification8.8
herbie shell --seed 2018362 +o rules:numerics
(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: 24.2s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 30.7b
Found 4 expressions with local error:
| 21.0b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 0.9b | (- (- 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:
| 34.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 14.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 5.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
4 calls. Slowest were:
| 69.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 60.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 44.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 28.0ms | (- (* b b) (* 4 (* a c))) |
98 calls. Slowest were:
| 457.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 411.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 382.0ms | (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 4.1b
Found 4 expressions with local error:
| 21.0b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 13.2b | (/ (* (* c 4) a) (* (* 2 a) (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 1.8b | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.4b | (- (* b b) (* 4 (* a c))) |
4 calls. Slowest were:
| 22.0ms | (/ (* (* c 4) a) (* (* 2 a) (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 16.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 6.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
4 calls. Slowest were:
| 192.0ms | (/ (* (* c 4) a) (* (* 2 a) (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 63.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 56.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 29.0ms | (- (* b b) (* 4 (* a c))) |
35 calls. Slowest were:
| 841.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)))))))) |
| 815.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)))))))) |
| 363.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 4.1b
Found 4 expressions with local error:
| 21.0b | (sqrt (fma -4 (* a c) (* b b))) |
| 10.4b | (/ (* (* c 4) a) (* 2 a)) |
| 1.8b | (- (sqrt (fma -4 (* a c) (* b b))) b) |
| 0.4b | (fma -4 (* a c) (* b b)) |
4 calls. Slowest were:
| 10.0ms | (/ (* (* c 4) a) (* 2 a)) |
| 10.0ms | (- (sqrt (fma -4 (* a c) (* b b))) b) |
| 2.0ms | (sqrt (fma -4 (* a c) (* b b))) |
4 calls. Slowest were:
| 56.0ms | (sqrt (fma -4 (* a c) (* b b))) |
| 49.0ms | (fma -4 (* a c) (* b b)) |
| 45.0ms | (- (sqrt (fma -4 (* a c) (* b b))) b) |
| 5.0ms | (/ (* (* c 4) a) (* 2 a)) |
26 calls. Slowest were:
| 755.0ms | (/ (* (* (* (* c 4) a) (* (* c 4) a)) (* (* c 4) a)) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
| 723.0ms | (- (log (* (* c 4) a)) (log (* 2 a))) |
| 41.0ms | (- (pow b 2) (* 4 (* a c))) |
8 alts after pruning (6 fresh and 2 done)
Merged error: 4.1b
Found 4 expressions with local error:
| 21.0b | (sqrt (fma (* a -4) c (* b b))) |
| 0.9b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.1b | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
| 0.0b | (fma (* a -4) c (* b b)) |
4 calls. Slowest were:
| 18.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
| 9.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 2.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 62.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
| 58.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 58.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 45.0ms | (fma (* a -4) c (* b b)) |
93 calls. Slowest were:
| 347.0ms | (- (* (- b) (- b)) (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b))))) |
| 184.0ms | (- (pow (- b) 3) (pow (sqrt (fma (* a -4) c (* b b))) 3)) |
| 119.0ms | (/ (* (* (- (- b) (sqrt (fma (* a -4) c (* b b)))) (- (- b) (sqrt (fma (* a -4) c (* b b))))) (- (- b) (sqrt (fma (* a -4) c (* b b))))) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 4.1b
86.6% (3.5b remaining)
Error of 8.8b against oracle of 5.3b and baseline of 31.3b