


Bits error versus x
if x < -2312.388795942151 or 12416.843762949418 < x Initial program 58.2
Simplified58.2
Taylor expanded around inf 0.0
Simplified0.0
if -2312.388795942151 < x < 12416.843762949418Initial program 0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019002 +o rules:numerics
(FPCore (x)
:name "Jmat.Real.dawson"
(* (/ (+ (+ (+ (+ (+ 1 (* 0.1049934947 (* x x))) (* 0.0424060604 (* (* x x) (* x x)))) (* 0.0072644182 (* (* (* x x) (* x x)) (* x x)))) (* 0.0005064034 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (+ (+ (+ (+ (+ (+ 1 (* 0.7715471019 (* x x))) (* 0.2909738639 (* (* x x) (* x x)))) (* 0.0694555761 (* (* (* x x) (* x x)) (* x x)))) (* 0.0140005442 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0008327945 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (* 2 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x))
Time bar (total: 30.3s)Debug log
1 calls. Slowest were:
| 82.0ms | (* (/ (+ (+ (+ (+ (+ 1 (* 0.1049934947 (* x x))) (* 0.0424060604 (* (* x x) (* x x)))) (* 0.0072644182 (* (* (* x x) (* x x)) (* x x)))) (* 0.0005064034 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (+ (+ (+ (+ (+ (+ 1 (* 0.7715471019 (* x x))) (* 0.2909738639 (* (* x x) (* x x)))) (* 0.0694555761 (* (* (* x x) (* x x)) (* x x)))) (* 0.0140005442 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0008327945 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (* 2 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 26.8b
Found 4 expressions with local error:
| 26.8b | (/ (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) |
| 13.0b | (* (/ (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) x) |
| 0.2b | (* (* x x) 0.0008327945) |
| 0.1b | (* 0.0072644182 (* x x)) |
4 calls. Slowest were:
| 302.0ms | (/ (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) |
| 212.0ms | (* (/ (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) x) |
| 5.0ms | (* (* x x) 0.0008327945) |
4 calls. Slowest were:
| 146.0ms | (* (/ (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) x) |
| 115.0ms | (/ (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) |
| 13.0ms | (* 0.0072644182 (* x x)) |
| 12.0ms | (* (* x x) 0.0008327945) |
66 calls. Slowest were:
| 584.0ms | (+ (* (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971))))) (- (* (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761))))) (* (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761))))))) |
| 560.0ms | (/ (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) |
| 558.0ms | (/ (sqrt (+ (fma 0.0005064034 (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* 0.0072644182 (* x x)) (* (* x x) (* x x))) (fma (* (* x x) (* x x)) 0.0424060604 (fma 0.1049934947 (* x x) 1)))) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0001789971)))) (+ (+ (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* x x) 0.0008327945)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (* (* (* x x) (* x x)) (* 2 0.0001789971)))) (fma (* (* (* x x) (* x x)) (* (* x x) (* x x))) 0.0140005442 (+ (fma 0.2909738639 (* (* x x) (* x x)) (fma (* x x) 0.7715471019 1)) (* (* (* x x) (* x x)) (* (* x x) 0.0694555761)))))) |
2 alts after pruning (1 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.1b | (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5) |
| 0.1b | (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)) |
| 0.0b | (/ 0.15298196345929327 (pow x 5)) |
| 0.0b | (+ (/ 0.15298196345929327 (pow x 5)) (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5))) |
4 calls. Slowest were:
| 12.0ms | (+ (/ 0.15298196345929327 (pow x 5)) (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5))) |
| 11.0ms | (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)) |
| 1.0ms | (/ 0.15298196345929327 (pow x 5)) |
4 calls. Slowest were:
| 49.0ms | (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)) |
| 47.0ms | (+ (/ 0.15298196345929327 (pow x 5)) (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5))) |
| 37.0ms | (/ 0.15298196345929327 (pow x 5)) |
| 27.0ms | (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5) |
33 calls. Slowest were:
| 305.0ms | (* (exp (/ 0.15298196345929327 (pow x 5))) (exp (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)))) |
| 227.0ms | (+ (* 0.15298196345929327 x) (* (pow x 5) (* 1 (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)))) |
| 131.0ms | (+ (* 0.15298196345929327 (/ 1 (pow x 5))) (+ (* 0.2514179000665375 (/ 1 (pow x 3))) (* 0.5 (/ 1 x)))) |
2 alts after pruning (0 fresh and 2 done)
Merged error: 0b
100% (0.0b remaining)
Error of 0.0b against oracle of 0.0b and baseline of 28.7b