


Bits error versus F



Bits error versus B



Bits error versus x
if F < 6.32901813840737e+172Initial program 9.7
Simplified9.6
rmApplied add-sqr-sqrt9.6
Applied unpow-prod-down9.6
if 6.32901813840737e+172 < F Initial program 39.9
Simplified39.8
rmApplied add-sqr-sqrt39.8
Applied unpow-prod-down39.8
Taylor expanded around inf 34.4
Simplified34.4
Final simplification13.0
herbie shell --seed 2019007 +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: 31.7s)Debug log
1 calls. Slowest were:
| 167.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: 13.6b
Found 4 expressions with local error:
| 19.7b | (pow (fma 2 x (fma F F 2)) -1/2) |
| 7.1b | (fma (pow (fma 2 x (fma F F 2)) -1/2) (/ F (sin B)) (/ (- x) (tan B))) |
| 0.2b | (/ (- x) (tan B)) |
| 0.1b | (/ F (sin B)) |
4 calls. Slowest were:
| 3.0ms | (/ (- x) (tan B)) |
| 3.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 0.0ms | (/ F (sin B)) |
4 calls. Slowest were:
| 279.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))) |
| 146.0ms | (/ (- x) (tan B)) |
| 61.0ms | (/ F (sin B)) |
54 calls. Slowest were:
| 249.0ms | (- (/ (* F (sqrt 1/2)) B) (+ (/ x B) (/ (* x (* F (sqrt 1/8))) B))) |
| 237.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)))))) |
| 213.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)))))) |
15 alts after pruning (15 fresh and 0 done)
Merged error: 13.0b
Found 4 expressions with local error:
| 19.7b | (pow (fma 2 x (fma F F 2)) -1/2) |
| 7.1b | (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:
| 6.0ms | (/ (- x) (tan B)) |
| 3.0ms | (* F (/ 1 (sin B))) |
| 2.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
4 calls. Slowest were:
| 248.0ms | (fma (pow (fma 2 x (fma F F 2)) -1/2) (* F (/ 1 (sin B))) (/ (- x) (tan B))) |
| 193.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 114.0ms | (/ (- x) (tan B)) |
| 63.0ms | (* F (/ 1 (sin B))) |
57 calls. Slowest were:
| 273.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))) |
| 222.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)))))) |
16 alts after pruning (16 fresh and 0 done)
Merged error: 13.0b
Found 4 expressions with local error:
| 19.7b | (pow (fma 2 x (fma F F 2)) -1/2) |
| 7.1b | (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:
| 12.0ms | (/ (* (cos B) (- x)) (sin B)) |
| 2.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 2.0ms | (* F (/ 1 (sin B))) |
4 calls. Slowest were:
| 233.0ms | (fma (pow (fma 2 x (fma F F 2)) -1/2) (* F (/ 1 (sin B))) (/ (* (cos B) (- x)) (sin B))) |
| 189.0ms | (pow (fma 2 x (fma F F 2)) -1/2) |
| 72.0ms | (/ (* (cos B) (- x)) (sin B)) |
| 63.0ms | (* F (/ 1 (sin B))) |
38 calls. Slowest were:
| 265.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))) |
| 257.0ms | (- (/ (* F (sqrt 1/2)) B) (+ (/ x B) (/ (* x (* F (sqrt 1/8))) B))) |
| 241.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)))))) |
17 alts after pruning (16 fresh and 1 done)
Merged error: 13.0b
Found 4 expressions with local error:
| 19.7b | (sqrt (fma 2 x (fma F F 2))) |
| 19.7b | (sqrt (fma 2 x (fma F F 2))) |
| 7.1b | (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:
| 7.0ms | (* (pow (sqrt (fma 2 x (fma F F 2))) -1/2) (pow (sqrt (fma 2 x (fma F F 2))) -1/2)) |
| 3.0ms | (sqrt (fma 2 x (fma F F 2))) |
| 1.0ms | (sqrt (fma 2 x (fma F F 2))) |
4 calls. Slowest were:
| 181.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))) |
| 64.0ms | (sqrt (fma 2 x (fma F F 2))) |
| 59.0ms | (sqrt (fma 2 x (fma F F 2))) |
| 35.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:
| 207.0ms | (- (/ (* F (sqrt 1/2)) B) (+ (/ x B) (/ (* x (* F (sqrt 1/8))) B))) |
| 189.0ms | (- (+ (sqrt 2) (* x (sqrt 1/2))) (* 1/2 (* (pow x 2) (sqrt 1/8)))) |
| 156.0ms | (- (+ (sqrt 2) (* x (sqrt 1/2))) (* 1/2 (* (pow x 2) (sqrt 1/8)))) |
17 alts after pruning (16 fresh and 1 done)
Merged error: 11.4b
33% (1.5b remaining)
Error of 13.0b against oracle of 11.5b and baseline of 13.7b