


Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus d
Initial program 28.5
rmApplied flip-+28.5
Applied associate-/l/28.5
Simplified0.6
rmApplied times-frac0.6
Simplified0.5
Simplified0.5
rmApplied pow10.5
Applied pow10.5
Applied pow-prod-down0.5
Simplified0.3
Final simplification0.3
herbie shell --seed 2019007 +o rules:numerics
(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: 53.4s)Debug log
1 calls. Slowest were:
| 59.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 28.4b
Found 4 expressions with local error:
| 28.3b | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 0.4b | (/ (+ (- 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:
| 26.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
| 12.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 2.0ms | (- (* b b) (* (* 3 a) c)) |
4 calls. Slowest were:
| 91.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)) |
| 67.0ms | (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) |
| 27.0ms | (- (* b b) (* (* 3 a) c)) |
| 11.0ms | (* (* 3 a) c) |
34 calls. Slowest were:
| 383.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)) |
| 118.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))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 0.5b
Found 4 expressions with local error:
| 0.4b | (/ (* 3 (* c a)) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
| 0.4b | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 0.3b | (* 3 (* c a)) |
| 0.3b | (* (* 3 a) c) |
4 calls. Slowest were:
| 30.0ms | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 20.0ms | (/ (* 3 (* c a)) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
| 4.0ms | (* (* 3 a) c) |
4 calls. Slowest were:
| 147.0ms | (/ (* 3 (* c a)) (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
| 101.0ms | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 16.0ms | (* 3 (* c a)) |
| 15.0ms | (* (* 3 a) c) |
84 calls. Slowest were:
| 2.5s | (/ (* (* (* 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.6s | (/ (* (* (* 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))))))) |
| 536.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (- (* b b) (* (* 3 a) c))) (sqrt (- (* b b) (* (* 3 a) c)))) (* (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.4b | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 0.3b | (/ (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (* c a)) |
| 0.3b | (* (* 3 a) c) |
| 0.3b | (/ 3 (/ (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (* c a))) |
4 calls. Slowest were:
| 18.0ms | (/ (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (* c a)) |
| 16.0ms | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 2.0ms | (* (* 3 a) c) |
4 calls. Slowest were:
| 122.0ms | (/ 3 (/ (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (* c a))) |
| 102.0ms | (/ (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) (* c a)) |
| 83.0ms | (* (* 3 a) (- (- b) (sqrt (- (* b b) (* (* 3 a) c))))) |
| 19.0ms | (* (* 3 a) c) |
90 calls. Slowest were:
| 625.0ms | (* (fma (- (sqrt (- (* b b) (* (* 3 a) c)))) (sqrt 1) (* (sqrt (- (* b b) (* (* 3 a) c))) (sqrt 1))) (* 3 a)) |
| 509.0ms | (/ (* (* (* (* 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)))))) (* (* (* c a) (* c a)) (* c a))) |
| 456.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (- (* b b) (* (* 3 a) c))) (sqrt (- (* b b) (* (* 3 a) c)))) (* (- b) (sqrt (- (* b b) (* (* 3 a) c)))))) |
9 alts after pruning (9 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.2b | (- (- b) (sqrt (fma -3 (* c a) (* b b)))) |
4 calls. Slowest were:
| 17.0ms | (* (/ 1 a) (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
| 13.0ms | (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b))))) |
| 8.0ms | (- (- b) (sqrt (fma -3 (* c a) (* b b)))) |
4 calls. Slowest were:
| 147.0ms | (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b))))) |
| 128.0ms | (* (/ 1 a) (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
| 58.0ms | (- (- b) (sqrt (fma -3 (* c a) (* b b)))) |
| 32.0ms | (fma -3 (* c a) (* b b)) |
113 calls. Slowest were:
| 517.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)))))) |
| 428.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)))))) |
| 315.0ms | (* (/ 1 a) (/ (* c a) (- (- b) (sqrt (fma -3 (* c a) (* b b)))))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 0b
0% (0.2b remaining)
Error of 0.3b against oracle of 0.1b and baseline of 0.3b