


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
Initial program 43.8
rmApplied flip-+43.8
Applied associate-/l/43.8
Simplified0.5
rmApplied times-frac0.5
Simplified0.5
Simplified0.5
rmApplied pow10.5
Applied pow10.5
Applied pow-prod-down0.5
Simplified0.2
rmApplied clear-num0.3
Final simplification0.3
herbie shell --seed 2018365 +o rules:numerics
(FPCore (a b c d)
:name "Cubic critical, medium range"
:pre (and (< 1.1102230246251565e-16 a 9007199254740992.0) (< 1.1102230246251565e-16 b 9007199254740992.0) (< 1.1102230246251565e-16 c 9007199254740992.0))
(/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)))
Time bar (total: 39.3s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 44.1b
Found 4 expressions with local error:
| 38.0b | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 0.3b | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
| 0.3b | (* (* 3 a) c) |
| 0.2b | (- (* b b) (* (* 3 a) c)) |
4 calls. Slowest were:
| 36.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
| 12.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 5.0ms | (- (* b b) (* (* 3 a) c)) |
4 calls. Slowest were:
| 85.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 74.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
| 42.0ms | (- (* b b) (* (* 3 a) c)) |
| 18.0ms | (* (* 3 a) c) |
34 calls. Slowest were:
| 469.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 3 a) c))) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 154.0ms | (+ (pow (- b) 3) (pow (sqrt (- (* b b) (* (* 3 a) c))) 3)) |
| 144.0ms | (/ (* (* (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (+ (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (+ (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (* (* (* 3 a) (* 3 a)) (* 3 a))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.2b
Found 4 expressions with local error:
| 0.4b | (/ (* 3 (* c a)) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
| 0.3b | (* 3 (* c a)) |
| 0.3b | (* (* 3 a) c) |
| 0.3b | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
4 calls. Slowest were:
| 28.0ms | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 25.0ms | (/ (* 3 (* c a)) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
| 4.0ms | (* (* 3 a) c) |
4 calls. Slowest were:
| 133.0ms | (/ (* 3 (* c a)) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
| 93.0ms | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 19.0ms | (* (* 3 a) c) |
| 17.0ms | (* 3 (* c a)) |
84 calls. Slowest were:
| 2.4s | (/ (* (* (* 3 (* c a)) (* 3 (* c a))) (* 3 (* c a))) (* (* (* (* 3 a) (* 3 a)) (* 3 a)) (* (* (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))))) |
| 1.9s | (/ (* (* (* 3 (* c a)) (* 3 (* c a))) (* 3 (* c a))) (* (* (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))))) |
| 628.0ms | (* (* 3 a) (fma (- (sqrt (cbrt (- (* b b) (* (* 3 a) c))))) (sqrt (* (cbrt (- (* b b) (* (* 3 a) c))) (cbrt (- (* b b) (* (* 3 a) c))))) (* (sqrt (cbrt (- (* b b) (* (* 3 a) c)))) (sqrt (* (cbrt (- (* b b) (* (* 3 a) c))) (cbrt (- (* b b) (* (* 3 a) c)))))))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.4b | (* (/ 1 a) (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
| 0.3b | (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b))))) |
| 0.2b | (fma -3 (* c a) (* b b)) |
| 0.1b | (- (- b) (sqrt (fma -3 (* c a) (* b b)))) |
4 calls. Slowest were:
| 18.0ms | (* (/ 1 a) (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
| 14.0ms | (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b))))) |
| 9.0ms | (- (- b) (sqrt (fma -3 (* c a) (* b b)))) |
4 calls. Slowest were:
| 142.0ms | (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b))))) |
| 132.0ms | (* (/ 1 a) (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
| 64.0ms | (- (- b) (sqrt (fma -3 (* c a) (* b b)))) |
| 55.0ms | (fma -3 (* c a) (* b b)) |
113 calls. Slowest were:
| 558.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma -3 (* c a) (* b b))) (sqrt (fma -3 (* c a) (* b b)))) (* (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
| 493.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma -3 (* c a) (* b b))) (sqrt (fma -3 (* c a) (* b b)))) (* (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
| 341.0ms | (+ (log (/ 1 a)) (- (log (* c a)) (log (- (- b) (sqrt (fma -3 (* c a) (* b b))))))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.2b | (fma (* a -3) c (* b b)) |
| 0.1b | (- (- b) (sqrt (fma (* a -3) c (* b b)))) |
| 0.1b | (/ c (- (- b) (sqrt (fma (* a -3) c (* b b))))) |
| 0.1b | (sqrt (fma (* a -3) c (* b b))) |
4 calls. Slowest were:
| 9.0ms | (- (- b) (sqrt (fma (* a -3) c (* b b)))) |
| 4.0ms | (/ c (- (- b) (sqrt (fma (* a -3) c (* b b))))) |
| 1.0ms | (sqrt (fma (* a -3) c (* b b))) |
4 calls. Slowest were:
| 124.0ms | (/ c (- (- b) (sqrt (fma (* a -3) c (* b b))))) |
| 63.0ms | (- (- b) (sqrt (fma (* a -3) c (* b b)))) |
| 60.0ms | (sqrt (fma (* a -3) c (* b b))) |
| 33.0ms | (fma (* a -3) c (* b b)) |
80 calls. Slowest were:
| 488.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma (* a -3) c (* b b))) (sqrt (fma (* a -3) c (* b b)))) (* (- b) (sqrt (fma (* a -3) c (* b b)))))) |
| 68.0ms | (fma 1 (- b) (- (* (sqrt (fma (* a -3) c (* b b))) (sqrt 1)))) |
| 65.0ms | (fma (sqrt (- b)) (sqrt (- b)) (- (* (sqrt (fma (* a -3) c (* b b))) (sqrt 1)))) |
6 alts after pruning (6 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