


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
if b < -1.127123378506132e+57Initial program 38.3
rmApplied add-cube-cbrt38.3
rmApplied add-sqr-sqrt38.3
Applied rem-sqrt-square38.3
Simplified38.3
Taylor expanded around inf 10.7
Simplified5.5
if -1.127123378506132e+57 < b < 2.0611675014481223Initial program 17.0
rmApplied add-cube-cbrt17.3
rmApplied add-sqr-sqrt17.3
Applied rem-sqrt-square17.3
Simplified17.1
if 2.0611675014481223 < b Initial program 54.8
rmApplied associate-*l*54.8
Taylor expanded around inf 16.6
Final simplification14.6
herbie shell --seed 2019007 +o rules:numerics
(FPCore (a b c d)
:name "Cubic critical"
(/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)))
Time bar (total: 16.9s)Debug log
1 calls. Slowest were:
| 108.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 35.3b
Found 4 expressions with local error:
| 25.7b | (sqrt (- (* b b) (* (* 3 a) c))) |
| 1.5b | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 0.8b | (- (* b b) (* (* 3 a) c)) |
| 0.2b | (* (* 3 a) c) |
4 calls. Slowest were:
| 11.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 6.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)))) |
| 55.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 30.0ms | (- (* b b) (* (* 3 a) c)) |
| 11.0ms | (* (* 3 a) c) |
25 calls. Slowest were:
| 356.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 3 a) c) 3))) |
| 181.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 3 a) c) (* (* 3 a) c)))) |
| 113.0ms | (/ (exp (* b b)) (exp (* (* 3 a) c))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 11.4b
Found 4 expressions with local error:
| 25.7b | (sqrt (- (* b b) (* 3 (* a c)))) |
| 1.5b | (+ (- b) (sqrt (- (* b b) (* 3 (* a c))))) |
| 0.8b | (- (* b b) (* 3 (* a c))) |
| 0.2b | (* 3 (* a c)) |
4 calls. Slowest were:
| 12.0ms | (+ (- b) (sqrt (- (* b b) (* 3 (* a c))))) |
| 6.0ms | (sqrt (- (* b b) (* 3 (* a c)))) |
| 4.0ms | (- (* b b) (* 3 (* a c))) |
4 calls. Slowest were:
| 62.0ms | (+ (- b) (sqrt (- (* b b) (* 3 (* a c))))) |
| 49.0ms | (sqrt (- (* b b) (* 3 (* a c)))) |
| 43.0ms | (- (* b b) (* 3 (* a c))) |
| 18.0ms | (* 3 (* a c)) |
25 calls. Slowest were:
| 258.0ms | (sqrt (- (pow (* b b) 3) (pow (* 3 (* a c)) 3))) |
| 205.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 3 (* a c)) (* 3 (* a c))))) |
| 116.0ms | (- b (* 3/2 (/ (* a c) b))) |
13 alts after pruning (12 fresh and 1 done)
Merged error: 11.4b
Found 4 expressions with local error:
| 25.7b | (sqrt (- (* b b) (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c))))) |
| 13.1b | (cbrt (* (* 3 a) c)) |
| 13.1b | (cbrt (* (* 3 a) c)) |
| 13.1b | (cbrt (* (* 3 a) c)) |
4 calls. Slowest were:
| 6.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c))))) |
| 3.0ms | (cbrt (* (* 3 a) c)) |
| 3.0ms | (cbrt (* (* 3 a) c)) |
4 calls. Slowest were:
| 129.0ms | (cbrt (* (* 3 a) c)) |
| 127.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c))))) |
| 116.0ms | (cbrt (* (* 3 a) c)) |
| 96.0ms | (cbrt (* (* 3 a) c)) |
21 calls. Slowest were:
| 413.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c))) 3))) |
| 311.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c))) (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c)))))) |
| 186.0ms | (* (exp (* -1/3 (+ (log (/ 1 a)) (log (/ 1 c))))) (cbrt 3)) |
12 alts after pruning (10 fresh and 2 done)
Merged error: 11.4b
Found 4 expressions with local error:
| 25.7b | (sqrt (- (* b b) (* a (* 3 c)))) |
| 1.5b | (+ (- b) (fabs (sqrt (- (* b b) (* a (* 3 c)))))) |
| 0.8b | (- (* b b) (* a (* 3 c))) |
| 0.2b | (* a (* 3 c)) |
4 calls. Slowest were:
| 4.0ms | (+ (- b) (fabs (sqrt (- (* b b) (* a (* 3 c)))))) |
| 4.0ms | (sqrt (- (* b b) (* a (* 3 c)))) |
| 2.0ms | (- (* b b) (* a (* 3 c))) |
4 calls. Slowest were:
| 44.0ms | (sqrt (- (* b b) (* a (* 3 c)))) |
| 42.0ms | (- (* b b) (* a (* 3 c))) |
| 14.0ms | (* a (* 3 c)) |
| 10.0ms | (+ (- b) (fabs (sqrt (- (* b b) (* a (* 3 c)))))) |
25 calls. Slowest were:
| 299.0ms | (sqrt (- (pow (* b b) 3) (pow (* a (* 3 c)) 3))) |
| 171.0ms | (sqrt (- (* (* b b) (* b b)) (* (* a (* 3 c)) (* a (* 3 c))))) |
| 127.0ms | (- b (* 3/2 (/ (* a c) b))) |
12 alts after pruning (9 fresh and 3 done)
Merged error: 10.6b
83.8% (3.7b remaining)
Error of 14.6b against oracle of 10.9b and baseline of 33.7b