


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
Results
Initial program 28.4
Simplified28.4
rmApplied flip--28.4
Applied associate-/l/28.4
Simplified0.6
rmApplied times-frac0.6
Simplified0.5
rmApplied pow10.5
Applied pow10.5
Applied pow-prod-down0.5
Simplified0.3
rmApplied flip3-+0.3
Final simplification0.3
herbie shell --seed 2019004
(FPCore (a b c d)
:name "Cubic critical, narrow range"
:pre (and (< 1.0536712127723509e-08 a 94906265.62425156) (< 1.0536712127723509e-08 b 94906265.62425156) (< 1.0536712127723509e-08 c 94906265.62425156))
(/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)))
Time bar (total: 15.8s)Debug log
1 calls. Slowest were:
| 46.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 27.4b
Found 4 expressions with local error:
| 27.3b | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 0.4b | (/ (- (sqrt (+ (* -3 (* c a)) (* b b))) b) (* 3 a)) |
| 0.3b | (* -3 (* c a)) |
| 0.3b | (+ (* -3 (* c a)) (* b b)) |
4 calls. Slowest were:
| 23.0ms | (/ (- (sqrt (+ (* -3 (* c a)) (* b b))) b) (* 3 a)) |
| 22.0ms | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 3.0ms | (+ (* -3 (* c a)) (* b b)) |
4 calls. Slowest were:
| 64.0ms | (/ (- (sqrt (+ (* -3 (* c a)) (* b b))) b) (* 3 a)) |
| 55.0ms | (+ (* -3 (* c a)) (* b b)) |
| 52.0ms | (- (sqrt (+ (* -3 (* c a)) (* b b))) b) |
| 13.0ms | (* -3 (* c a)) |
26 calls. Slowest were:
| 587.0ms | (- (* (sqrt (+ (* -3 (* c a)) (* b b))) (sqrt (+ (* -3 (* c a)) (* b b)))) (* b b)) |
| 101.0ms | (/ (* (* (- (sqrt (+ (* -3 (* c a)) (* b b))) b) (- (sqrt (+ (* -3 (* c a)) (* b b))) b)) (- (sqrt (+ (* -3 (* c a)) (* b b))) b)) (* (* (* 3 a) (* 3 a)) (* 3 a))) |
| 58.0ms | (- (pow (sqrt (+ (* -3 (* c a)) (* b b))) 3) (pow b 3)) |
5 alts after pruning (5 fresh and 0 done)
Merged error: 0.6b
Found 4 expressions with local error:
| 0.4b | (/ (* (* a c) -3) (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) |
| 0.3b | (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)) |
| 0.3b | (* (* a c) -3) |
| 0.3b | (* -3 (* c a)) |
4 calls. Slowest were:
| 19.0ms | (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)) |
| 12.0ms | (/ (* (* a c) -3) (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) |
| 4.0ms | (* (* a c) -3) |
4 calls. Slowest were:
| 149.0ms | (/ (* (* a c) -3) (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) |
| 119.0ms | (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)) |
| 13.0ms | (* -3 (* c a)) |
| 11.0ms | (* (* a c) -3) |
32 calls. Slowest were:
| 698.0ms | (/ (* (* (* (* a c) -3) (* (* a c) -3)) (* (* a c) -3)) (* (* (* (* 3 a) (* 3 a)) (* 3 a)) (* (* (+ (sqrt (+ (* -3 (* c a)) (* b b))) b) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)))) |
| 647.0ms | (/ (* (* (* (* a c) -3) (* (* a c) -3)) (* (* a c) -3)) (* (* (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)) (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) (* (* 3 a) (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)))) |
| 568.0ms | (- (log (* (* a c) -3)) (+ (log (* 3 a)) (log (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)))) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.4b | (* (* 1/3 c) (/ -3 (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) |
| 0.3b | (* 1/3 c) |
| 0.3b | (* -3 (* c a)) |
| 0.3b | (+ (* -3 (* c a)) (* b b)) |
4 calls. Slowest were:
| 9.0ms | (* (* 1/3 c) (/ -3 (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) |
| 5.0ms | (+ (* -3 (* c a)) (* b b)) |
| 2.0ms | (* -3 (* c a)) |
4 calls. Slowest were:
| 179.0ms | (* (* 1/3 c) (/ -3 (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) |
| 35.0ms | (+ (* -3 (* c a)) (* b b)) |
| 10.0ms | (* -3 (* c a)) |
| 8.0ms | (* 1/3 c) |
23 calls. Slowest were:
| 527.0ms | (* (* (* (* 1/3 c) (* 1/3 c)) (* 1/3 c)) (* (* (/ -3 (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)) (/ -3 (+ (sqrt (+ (* -3 (* c a)) (* b b))) b))) (/ -3 (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)))) |
| 128.0ms | (+ (* (sqrt (+ (* -3 (* c a)) (* b b))) (sqrt (+ (* -3 (* c a)) (* b b)))) (- (* b b) (* (sqrt (+ (* -3 (* c a)) (* b b))) b))) |
| 110.0ms | (+ (log (* 1/3 c)) (log (/ -3 (+ (sqrt (+ (* -3 (* c a)) (* b b))) b)))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.3b | (* -3 (* c a)) |
| 0.3b | (+ (* b b) (* -3 (* c a))) |
| 0.2b | (/ (- c) (+ b (sqrt (+ (* b b) (* -3 (* c a)))))) |
| 0.2b | (+ b (sqrt (+ (* b b) (* -3 (* c a))))) |
4 calls. Slowest were:
| 6.0ms | (/ (- c) (+ b (sqrt (+ (* b b) (* -3 (* c a)))))) |
| 2.0ms | (+ (* b b) (* -3 (* c a))) |
| 1.0ms | (+ b (sqrt (+ (* b b) (* -3 (* c a))))) |
4 calls. Slowest were:
| 132.0ms | (/ (- c) (+ b (sqrt (+ (* b b) (* -3 (* c a)))))) |
| 46.0ms | (+ b (sqrt (+ (* b b) (* -3 (* c a))))) |
| 27.0ms | (+ (* b b) (* -3 (* c a))) |
| 11.0ms | (* -3 (* c a)) |
48 calls. Slowest were:
| 318.0ms | (/ (* (* (- c) (- c)) (- c)) (* (* (+ b (sqrt (+ (* b b) (* -3 (* c a))))) (+ b (sqrt (+ (* b b) (* -3 (* c a)))))) (+ b (sqrt (+ (* b b) (* -3 (* c a))))))) |
| 119.0ms | (+ (* b b) (- (* (sqrt (+ (* b b) (* -3 (* c a)))) (sqrt (+ (* b b) (* -3 (* c a))))) (* b (sqrt (+ (* b b) (* -3 (* c a))))))) |
| 110.0ms | (- (* 2 b) (* 3/2 (/ (* a c) b))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0b
0% (0.3b remaining)
Error of 0.3b against oracle of 0.0b and baseline of 0.3b