


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
if b < -2.5788381684060953e+153Initial program 60.8
Taylor expanded around -inf 2.4
Simplified2.4
if -2.5788381684060953e+153 < b < 1.2378604166113681e-144Initial program 11.0
rmApplied associate-/r*11.0
Taylor expanded around -inf 11.0
if 1.2378604166113681e-144 < b Initial program 48.9
Taylor expanded around inf 13.3
Final simplification11.1
herbie shell --seed 2018365 +o rules:numerics
(FPCore (a b c d)
:name "Cubic critical"
(/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)))
Time bar (total: 13.7s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 32.2b
Found 4 expressions with local error:
| 24.6b | (sqrt (- (* b b) (* (* 3 a) c))) |
| 1.0b | (- (* b b) (* (* 3 a) c)) |
| 0.8b | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 0.3b | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
4 calls. Slowest were:
| 17.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
| 6.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 6.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
4 calls. Slowest were:
| 79.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
| 54.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 54.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 36.0ms | (- (* b b) (* (* 3 a) c)) |
40 calls. Slowest were:
| 413.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 3 a) c) 3))) |
| 383.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 3 a) c))) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 196.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 3 a) c) (* (* 3 a) c)))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 24.6b | (sqrt (- (* b b) (* (* 3 a) c))) |
| 1.0b | (- (* b b) (* (* 3 a) c)) |
| 0.8b | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 0.2b | (* (* 3 a) c) |
4 calls. Slowest were:
| 6.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 3.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 2.0ms | (- (* b b) (* (* 3 a) c)) |
4 calls. Slowest were:
| 75.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 61.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 39.0ms | (- (* b b) (* (* 3 a) c)) |
| 18.0ms | (* (* 3 a) c) |
25 calls. Slowest were:
| 361.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 3 a) c) 3))) |
| 147.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 3 a) c) (* (* 3 a) c)))) |
| 141.0ms | (- b (* 3/2 (/ (* a c) b))) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 24.6b | (sqrt (- (* b b) (* 3 (* a c)))) |
| 1.0b | (- (* b b) (* 3 (* a c))) |
| 0.8b | (+ (- b) (sqrt (- (* b b) (* 3 (* a c))))) |
| 0.2b | (* 3 (* a c)) |
4 calls. Slowest were:
| 9.0ms | (+ (- b) (sqrt (- (* b b) (* 3 (* a c))))) |
| 3.0ms | (sqrt (- (* b b) (* 3 (* a c)))) |
| 2.0ms | (- (* b b) (* 3 (* a c))) |
4 calls. Slowest were:
| 61.0ms | (+ (- b) (sqrt (- (* b b) (* 3 (* a c))))) |
| 53.0ms | (sqrt (- (* b b) (* 3 (* a c)))) |
| 42.0ms | (- (* b b) (* 3 (* a c))) |
| 19.0ms | (* 3 (* a c)) |
25 calls. Slowest were:
| 387.0ms | (sqrt (- (pow (* b b) 3) (pow (* 3 (* a c)) 3))) |
| 228.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 3 (* a c)) (* 3 (* a c))))) |
| 171.0ms | (- b (* 3/2 (/ (* a c) b))) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 24.6b | (sqrt (fma (* -3 a) c (* b b))) |
| 0.8b | (- (sqrt (fma (* -3 a) c (* b b))) b) |
| 0.3b | (* 1/3 (/ (- (sqrt (fma (* -3 a) c (* b b))) b) a)) |
| 0.1b | (/ (- (sqrt (fma (* -3 a) c (* b b))) b) a) |
4 calls. Slowest were:
| 20.0ms | (/ (- (sqrt (fma (* -3 a) c (* b b))) b) a) |
| 14.0ms | (* 1/3 (/ (- (sqrt (fma (* -3 a) c (* b b))) b) a)) |
| 10.0ms | (- (sqrt (fma (* -3 a) c (* b b))) b) |
4 calls. Slowest were:
| 86.0ms | (/ (- (sqrt (fma (* -3 a) c (* b b))) b) a) |
| 80.0ms | (* 1/3 (/ (- (sqrt (fma (* -3 a) c (* b b))) b) a)) |
| 70.0ms | (- (sqrt (fma (* -3 a) c (* b b))) b) |
| 65.0ms | (sqrt (fma (* -3 a) c (* b b))) |
30 calls. Slowest were:
| 814.0ms | (- (* (sqrt (fma (* -3 a) c (* b b))) (sqrt (fma (* -3 a) c (* b b)))) (* b b)) |
| 128.0ms | (- (pow (sqrt (fma (* -3 a) c (* b b))) 3) (pow b 3)) |
| 31.0ms | (* (cbrt (- (sqrt (fma (* -3 a) c (* b b))) b)) (cbrt (- (sqrt (fma (* -3 a) c (* b b))) b))) |
13 alts after pruning (10 fresh and 3 done)
Merged error: 6.6b
80.9% (5.2b remaining)
Error of 11.1b against oracle of 5.9b and baseline of 33.2b