


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
Initial program 44.0
Simplified44.0
rmApplied flip--44.0
Applied associate-/l/44.0
Simplified0.5
rmApplied associate-/r*0.5
Taylor expanded around inf 0.2
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied sqrt-prod0.3
Applied fma-def0.2
Final simplification0.2
herbie shell --seed 2019004 +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: 26.7s)Debug log
1 calls. Slowest were:
| 43.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 40.8b
Found 4 expressions with local error:
| 34.9b | (- (sqrt (fma -3 (* c a) (* b b))) b) |
| 0.3b | (/ (- (sqrt (fma -3 (* c a) (* b b))) b) (* 3 a)) |
| 0.2b | (fma -3 (* c a) (* b b)) |
| 0.1b | (sqrt (fma -3 (* c a) (* b b))) |
4 calls. Slowest were:
| 10.0ms | (/ (- (sqrt (fma -3 (* c a) (* b b))) b) (* 3 a)) |
| 5.0ms | (- (sqrt (fma -3 (* c a) (* b b))) b) |
| 1.0ms | (sqrt (fma -3 (* c a) (* b b))) |
4 calls. Slowest were:
| 53.0ms | (/ (- (sqrt (fma -3 (* c a) (* b b))) b) (* 3 a)) |
| 53.0ms | (sqrt (fma -3 (* c a) (* b b))) |
| 44.0ms | (- (sqrt (fma -3 (* c a) (* b b))) b) |
| 32.0ms | (fma -3 (* c a) (* b b)) |
34 calls. Slowest were:
| 701.0ms | (- (* (sqrt (fma -3 (* c a) (* b b))) (sqrt (fma -3 (* c a) (* b b)))) (* b b)) |
| 111.0ms | (- (pow (sqrt (fma -3 (* c a) (* b b))) 3) (pow b 3)) |
| 81.0ms | (/ (* (* (- (sqrt (fma -3 (* c a) (* b b))) b) (- (sqrt (fma -3 (* c a) (* b b))) b)) (- (sqrt (fma -3 (* c a) (* b b))) b)) (* (* (* 3 a) (* 3 a)) (* 3 a))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 0.6b
Found 4 expressions with local error:
| 0.4b | (/ (* (* a c) -3) (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b))) |
| 0.4b | (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 0.3b | (* (* a c) -3) |
| 0.2b | (fma -3 (* c a) (* b b)) |
4 calls. Slowest were:
| 9.0ms | (/ (* (* a c) -3) (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b))) |
| 8.0ms | (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 2.0ms | (* (* a c) -3) |
4 calls. Slowest were:
| 143.0ms | (/ (* (* a c) -3) (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b))) |
| 100.0ms | (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 30.0ms | (fma -3 (* c a) (* b b)) |
| 20.0ms | (* (* a c) -3) |
32 calls. Slowest were:
| 609.0ms | (/ (* (* (* (* a c) -3) (* (* a c) -3)) (* (* a c) -3)) (* (* (* (* 3 a) (* 3 a)) (* 3 a)) (* (* (+ (sqrt (fma -3 (* c a) (* b b))) b) (+ (sqrt (fma -3 (* c a) (* b b))) b)) (+ (sqrt (fma -3 (* c a) (* b b))) b)))) |
| 591.0ms | (- (log (* (* a c) -3)) (+ (log (* 3 a)) (log (+ (sqrt (fma -3 (* c a) (* b b))) b)))) |
| 570.0ms | (/ (* (* (* (* a c) -3) (* (* a c) -3)) (* (* a c) -3)) (* (* (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b)) (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b))) (* (* 3 a) (+ (sqrt (fma -3 (* c a) (* b b))) b)))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.3b | (/ (* (* a c) -3) (* 3 a)) |
| 0.3b | (* (* a c) -3) |
| 0.2b | (fma -3 (* c a) (* b b)) |
| 0.1b | (/ (/ (* (* a c) -3) (* 3 a)) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
4 calls. Slowest were:
| 17.0ms | (/ (/ (* (* a c) -3) (* 3 a)) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 5.0ms | (/ (* (* a c) -3) (* 3 a)) |
| 2.0ms | (* (* a c) -3) |
4 calls. Slowest were:
| 111.0ms | (/ (/ (* (* a c) -3) (* 3 a)) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 49.0ms | (fma -3 (* c a) (* b b)) |
| 18.0ms | (* (* a c) -3) |
| 10.0ms | (/ (* (* a c) -3) (* 3 a)) |
61 calls. Slowest were:
| 957.0ms | (/ (/ (* (* (* (* a c) -3) (* (* a c) -3)) (* (* a c) -3)) (* (* (* 3 a) (* 3 a)) (* 3 a))) (* (* (+ (sqrt (fma -3 (* c a) (* b b))) b) (+ (sqrt (fma -3 (* c a) (* b b))) b)) (+ (sqrt (fma -3 (* c a) (* b b))) b))) |
| 468.0ms | (/ (sqrt (/ (* (* a c) -3) (* 3 a))) (* (cbrt (+ (sqrt (fma -3 (* c a) (* b b))) b)) (cbrt (+ (sqrt (fma -3 (* c a) (* b b))) b)))) |
| 432.0ms | (/ (* (* (/ (* (* a c) -3) (* 3 a)) (/ (* (* a c) -3) (* 3 a))) (/ (* (* a c) -3) (* 3 a))) (* (* (+ (sqrt (fma -3 (* c a) (* b b))) b) (+ (sqrt (fma -3 (* c a) (* b b))) b)) (+ (sqrt (fma -3 (* c a) (* b b))) b))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.2b | (fma -3 (* c a) (* b b)) |
| 0.1b | (/ (- c) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 0.1b | (+ (sqrt (fma -3 (* c a) (* b b))) b) |
| 0.1b | (sqrt (fma -3 (* c a) (* b b))) |
4 calls. Slowest were:
| 17.0ms | (/ (- c) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 5.0ms | (+ (sqrt (fma -3 (* c a) (* b b))) b) |
| 1.0ms | (sqrt (fma -3 (* c a) (* b b))) |
4 calls. Slowest were:
| 138.0ms | (/ (- c) (+ (sqrt (fma -3 (* c a) (* b b))) b)) |
| 47.0ms | (+ (sqrt (fma -3 (* c a) (* b b))) b) |
| 40.0ms | (sqrt (fma -3 (* c a) (* b b))) |
| 28.0ms | (fma -3 (* c a) (* b b)) |
56 calls. Slowest were:
| 305.0ms | (/ (* (* (- c) (- c)) (- c)) (* (* (+ (sqrt (fma -3 (* c a) (* b b))) b) (+ (sqrt (fma -3 (* c a) (* b b))) b)) (+ (sqrt (fma -3 (* c a) (* b b))) b))) |
| 292.0ms | (+ (* (sqrt (fma -3 (* c a) (* b b))) (sqrt (fma -3 (* c a) (* b b)))) (- (* b b) (* (sqrt (fma -3 (* c a) (* b b))) b))) |
| 49.0ms | (- (pow b 2) (* 3 (* a c))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 0b
0% (0.2b remaining)
Error of 0.2b against oracle of 0.0b and baseline of 0.2b