


Bits error versus F



Bits error versus B



Bits error versus x
Initial program 13.0
Simplified13.0
rmApplied add-sqr-sqrt13.0
Applied unpow-prod-down13.0
rmApplied *-commutative13.0
Final simplification13.0
herbie shell --seed 2019002 +o rules:numerics
(FPCore (F B x)
:name "VandenBroeck and Keller, Equation (23)"
(+ (- (* x (/ 1 (tan B)))) (* (/ F (sin B)) (pow (+ (+ (* F F) 2) (* 2 x)) (- (/ 1 2))))))
Time bar (total: 32.0s)Debug log
1 calls. Slowest were:
| 159.0ms | (+ (- (* x (/ 1 (tan B)))) (* (/ F (sin B)) (pow (+ (+ (* F F) 2) (* 2 x)) (- (/ 1 2))))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 12.5b
Found 4 expressions with local error:
| 16.5b | (pow (fma 2 x (fma F F 2)) -1/2) |
| 4.2b | (fma (pow (fma 2 x (fma F F 2)) -1/2) (/ F (sin B)) (/ (- x) (tan B))) |
| 0.2b | (/ (- x) (tan B)) |
| 0.2b | (/ F (sin B)) |
4 calls. Slowest were:
| 6.0ms | (/ (- x) (tan B)) |
| 5.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 1.0ms | (/ F (sin B)) |
4 calls. Slowest were:
| 233.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 220.0ms | (fma (pow (fma 2 x (fma F F 2)) -1/2) (/ F (sin B)) (/ (- x) (tan B))) |
| 140.0ms | (/ (- x) (tan B)) |
| 87.0ms | (/ F (sin B)) |
54 calls. Slowest were:
| 282.0ms | (- (/ (* F (sqrt 1/2)) B) (+ (/ x B) (/ (* x (* F (sqrt 1/8))) B))) |
| 280.0ms | (- (+ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) (* 3/8 (/ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) (pow x 2)))) (* 1/2 (/ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) x))) |
| 258.0ms | (- (/ (* F (exp (* -1/2 (- (log 2) (log (/ 1 x)))))) (sin B)) (+ (/ (* x (cos B)) (sin B)) (* 1/2 (/ (* F (exp (* -1/2 (- (log 2) (log (/ 1 x)))))) (* x (sin B)))))) |
15 alts after pruning (15 fresh and 0 done)
Merged error: 12.1b
Found 4 expressions with local error:
| 16.5b | (pow (fma 2 x (fma F F 2)) -1/2) |
| 4.2b | (fma (pow (fma 2 x (fma F F 2)) -1/2) (* F (/ 1 (sin B))) (/ (- x) (tan B))) |
| 0.2b | (* F (/ 1 (sin B))) |
| 0.2b | (/ (- x) (tan B)) |
4 calls. Slowest were:
| 3.0ms | (/ (- x) (tan B)) |
| 1.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 1.0ms | (* F (/ 1 (sin B))) |
4 calls. Slowest were:
| 234.0ms | (fma (pow (fma 2 x (fma F F 2)) -1/2) (* F (/ 1 (sin B))) (/ (- x) (tan B))) |
| 173.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 122.0ms | (/ (- x) (tan B)) |
| 56.0ms | (* F (/ 1 (sin B))) |
57 calls. Slowest were:
| 264.0ms | (- (/ (* F (sqrt 1/2)) B) (+ (/ x B) (/ (* x (* F (sqrt 1/8))) B))) |
| 251.0ms | (- (+ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) (* 3/8 (/ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) (pow x 2)))) (* 1/2 (/ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) x))) |
| 249.0ms | (- (/ (* F (exp (* -1/2 (- (log 2) (log (/ 1 x)))))) (sin B)) (+ (/ (* x (cos B)) (sin B)) (* 1/2 (/ (* F (exp (* -1/2 (- (log 2) (log (/ 1 x)))))) (* x (sin B)))))) |
14 alts after pruning (14 fresh and 0 done)
Merged error: 12.1b
Found 4 expressions with local error:
| 16.5b | (pow (fma 2 x (fma F F 2)) -1/2) |
| 4.2b | (fma (pow (fma 2 x (fma F F 2)) -1/2) (* F (/ 1 (sin B))) (/ (* (cos B) (- x)) (sin B))) |
| 0.2b | (* F (/ 1 (sin B))) |
| 0.2b | (/ (* (cos B) (- x)) (sin B)) |
4 calls. Slowest were:
| 11.0ms | (/ (* (cos B) (- x)) (sin B)) |
| 1.0ms | (* F (/ 1 (sin B))) |
| 1.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
4 calls. Slowest were:
| 206.0ms | (fma (pow (fma 2 x (fma F F 2)) -1/2) (* F (/ 1 (sin B))) (/ (* (cos B) (- x)) (sin B))) |
| 185.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 84.0ms | (/ (* (cos B) (- x)) (sin B)) |
| 61.0ms | (* F (/ 1 (sin B))) |
38 calls. Slowest were:
| 291.0ms | (- (+ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) (* 3/8 (/ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) (pow x 2)))) (* 1/2 (/ (exp (* -1/2 (- (log -2) (log (/ -1 x))))) x))) |
| 277.0ms | (- (/ (* F (sqrt 1/2)) B) (+ (/ x B) (/ (* x (* F (sqrt 1/8))) B))) |
| 277.0ms | (- (/ (* (exp (* -1/2 (- (log -2) (log (/ -1 x))))) F) (sin B)) (+ (/ (* x (cos B)) (sin B)) (* 1/2 (/ (* (exp (* -1/2 (- (log -2) (log (/ -1 x))))) F) (* x (sin B)))))) |
14 alts after pruning (14 fresh and 0 done)
Merged error: 12.1b
Found 4 expressions with local error:
| 16.5b | (sqrt (fma 2 x (fma F F 2))) |
| 16.5b | (sqrt (fma 2 x (fma F F 2))) |
| 4.2b | (fma (* (pow (sqrt (fma 2 x (fma F F 2))) -1/2) (pow (sqrt (fma 2 x (fma F F 2))) -1/2)) (/ F (sin B)) (/ (- x) (tan B))) |
| 0.7b | (* (pow (sqrt (fma 2 x (fma F F 2))) -1/2) (pow (sqrt (fma 2 x (fma F F 2))) -1/2)) |
4 calls. Slowest were:
| 15.0ms | (* (pow (sqrt (fma 2 x (fma F F 2))) -1/2) (pow (sqrt (fma 2 x (fma F F 2))) -1/2)) |
| 2.0ms | (sqrt (fma 2 x (fma F F 2))) |
| 2.0ms | (sqrt (fma 2 x (fma F F 2))) |
4 calls. Slowest were:
| 140.0ms | (fma (* (pow (sqrt (fma 2 x (fma F F 2))) -1/2) (pow (sqrt (fma 2 x (fma F F 2))) -1/2)) (/ F (sin B)) (/ (- x) (tan B))) |
| 55.0ms | (sqrt (fma 2 x (fma F F 2))) |
| 54.0ms | (sqrt (fma 2 x (fma F F 2))) |
| 32.0ms | (* (pow (sqrt (fma 2 x (fma F F 2))) -1/2) (pow (sqrt (fma 2 x (fma F F 2))) -1/2)) |
50 calls. Slowest were:
| 342.0ms | (- (/ (* F (sqrt 1/2)) B) (+ (/ x B) (/ (* x (* F (sqrt 1/8))) B))) |
| 273.0ms | (- (+ (sqrt 2) (* x (sqrt 1/2))) (* 1/2 (* (pow x 2) (sqrt 1/8)))) |
| 249.0ms | (- (+ (sqrt 2) (* x (sqrt 1/2))) (* 1/2 (* (pow x 2) (sqrt 1/8)))) |
17 alts after pruning (17 fresh and 0 done)
Merged error: 10.9b
0% (2.0b remaining)
Error of 13.0b against oracle of 11.0b and baseline of 13.0b