


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
Results
if b < -1.127123378506132e+57Initial program 38.3
Taylor expanded around -inf 10.8
if -1.127123378506132e+57 < b < 2.0611675014481223Initial program 17.0
Simplified17.1
rmApplied associate-*r*17.1
if 2.0611675014481223 < b Initial program 54.8
Taylor expanded around inf 16.6
Final simplification15.7
herbie shell --seed 2019007
(FPCore (a b c d)
:name "Cubic critical"
(/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)))
Time bar (total: 14.2s)Debug log
1 calls. Slowest were:
| 102.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:
| 5.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 2.0ms | (* (* 3 a) c) |
| 2.0ms | (- (* b b) (* (* 3 a) c)) |
4 calls. Slowest were:
| 71.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 67.0ms | (sqrt (- (* b b) (* (* 3 a) c))) |
| 33.0ms | (- (* b b) (* (* 3 a) c)) |
| 12.0ms | (* (* 3 a) c) |
21 calls. Slowest were:
| 422.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 3 a) c) 3))) |
| 167.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 3 a) c) (* (* 3 a) c)))) |
| 133.0ms | (- b (* 3/2 (/ (* a c) b))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 11.4b
Found 4 expressions with local error:
| 25.7b | (sqrt (+ (* -3 (* c a)) (* b b))) |
| 1.5b | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 0.8b | (+ (* -3 (* c a)) (* b b)) |
| 0.2b | (* -3 (* c a)) |
4 calls. Slowest were:
| 21.0ms | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 5.0ms | (+ (* -3 (* c a)) (* b b)) |
| 4.0ms | (sqrt (+ (* -3 (* c a)) (* b b))) |
4 calls. Slowest were:
| 52.0ms | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 51.0ms | (sqrt (+ (* -3 (* c a)) (* b b))) |
| 49.0ms | (+ (* -3 (* c a)) (* b b)) |
| 13.0ms | (* -3 (* c a)) |
19 calls. Slowest were:
| 199.0ms | (sqrt (- (* (* -3 (* c a)) (* -3 (* c a))) (* (* b b) (* b b)))) |
| 178.0ms | (sqrt (+ (pow (* -3 (* c a)) 3) (pow (* b b) 3))) |
| 46.0ms | (- (pow b 2) (* 3 (* a c))) |
12 alts after pruning (11 fresh and 1 done)
Merged error: 11.4b
Found 4 expressions with local error:
| 25.7b | (sqrt (+ (* (* -3 c) a) (* b b))) |
| 1.5b | (- (sqrt (+ (* (* -3 c) a) (* b b))) b) |
| 0.8b | (+ (* (* -3 c) a) (* b b)) |
| 0.2b | (* (* -3 c) a) |
4 calls. Slowest were:
| 29.0ms | (- (sqrt (+ (* (* -3 c) a) (* b b))) b) |
| 7.0ms | (+ (* (* -3 c) a) (* b b)) |
| 4.0ms | (* (* -3 c) a) |
4 calls. Slowest were:
| 70.0ms | (sqrt (+ (* (* -3 c) a) (* b b))) |
| 59.0ms | (- (sqrt (+ (* (* -3 c) a) (* b b))) b) |
| 42.0ms | (+ (* (* -3 c) a) (* b b)) |
| 15.0ms | (* (* -3 c) a) |
19 calls. Slowest were:
| 207.0ms | (sqrt (+ (pow (* (* -3 c) a) 3) (pow (* b b) 3))) |
| 198.0ms | (sqrt (- (* (* (* -3 c) a) (* (* -3 c) a)) (* (* b b) (* b b)))) |
| 50.0ms | (- (pow b 2) (* 3 (* a c))) |
11 alts after pruning (9 fresh and 2 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:
| 4.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:
| 136.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c))))) |
| 131.0ms | (cbrt (* (* 3 a) c)) |
| 123.0ms | (cbrt (* (* 3 a) c)) |
| 110.0ms | (cbrt (* (* 3 a) c)) |
21 calls. Slowest were:
| 327.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* (cbrt (* (* 3 a) c)) (cbrt (* (* 3 a) c))) (cbrt (* (* 3 a) c))) 3))) |
| 228.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)))))) |
| 178.0ms | (* (exp (* -1/3 (+ (log (/ 1 a)) (log (/ 1 c))))) (cbrt 3)) |
11 alts after pruning (8 fresh and 3 done)
Merged error: 11.4b
83.3% (3.6b remaining)
Error of 15.7b against oracle of 12.0b and baseline of 33.7b