


Bits error versus x
if x < -1.133500920034836 or 1.1444028190617166 < x Initial program 58.2
Simplified58.1
Taylor expanded around inf 0.3
Simplified0.3
if -1.133500920034836 < x < 1.1444028190617166Initial program 0.0
Simplified0.0
Taylor expanded around 0 0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019008 +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: 31.3s)Debug log
1 calls. Slowest were:
| 93.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) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 28.0b
Found 4 expressions with local error:
| 28.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)))))) |
| 13.7b | (* (/ (+ (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.1b | (* 0.0072644182 (* x x)) |
| 0.1b | (* (* x x) (* x x)) |
4 calls. Slowest were:
| 280.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)))))) |
| 222.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) |
| 3.0ms | (* (* x x) (* x x)) |
4 calls. Slowest were:
| 151.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) |
| 118.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 | (* (* x x) (* x x)) |
| 11.0ms | (* 0.0072644182 (* x x)) |
74 calls. Slowest were:
| 776.0ms | (/ (cbrt (+ (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)))) (sqrt (+ (+ (* (* (* (* 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))))))) |
| 586.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))))))) |
| 562.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 (2 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 11.8b | (- (* (pow x 4) 0.265709700396151) (fma 0.6665536072 (* x x) -1)) |
| 0.1b | (fma 0.6665536072 (* x x) -1) |
| 0.1b | (* (pow x 4) 0.265709700396151) |
| 0.0b | (* (- (* (pow x 4) 0.265709700396151) (fma 0.6665536072 (* x x) -1)) x) |
4 calls. Slowest were:
| 13.0ms | (* (- (* (pow x 4) 0.265709700396151) (fma 0.6665536072 (* x x) -1)) x) |
| 9.0ms | (- (* (pow x 4) 0.265709700396151) (fma 0.6665536072 (* x x) -1)) |
| 4.0ms | (* (pow x 4) 0.265709700396151) |
4 calls. Slowest were:
| 28.0ms | (* (- (* (pow x 4) 0.265709700396151) (fma 0.6665536072 (* x x) -1)) x) |
| 27.0ms | (fma 0.6665536072 (* x x) -1) |
| 24.0ms | (- (* (pow x 4) 0.265709700396151) (fma 0.6665536072 (* x x) -1)) |
| 15.0ms | (* (pow x 4) 0.265709700396151) |
26 calls. Slowest were:
| 277.0ms | (fma (- (fma 0.6665536072 (* x x) -1)) 1 (* (fma 0.6665536072 (* x x) -1) 1)) |
| 209.0ms | (* (- (pow (* (pow x 4) 0.265709700396151) 3) (pow (fma 0.6665536072 (* x x) -1) 3)) x) |
| 193.0ms | (* (- (* (* (pow x 4) 0.265709700396151) (* (pow x 4) 0.265709700396151)) (* (fma 0.6665536072 (* x x) -1) (fma 0.6665536072 (* x x) -1))) x) |
2 alts after pruning (1 fresh and 1 done)
Merged error: 0b
Found 3 expressions with local error:
| 0.1b | (/ 0.15298196345929327 (pow x 5)) |
| 0.1b | (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5) |
| 0.1b | (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)) |
3 calls. Slowest were:
| 4.0ms | (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)) |
| 1.0ms | (/ 0.15298196345929327 (pow x 5)) |
| 0.0ms | (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5) |
3 calls. Slowest were:
| 58.0ms | (* (/ 1 x) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)) |
| 32.0ms | (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5) |
| 23.0ms | (/ 0.15298196345929327 (pow x 5)) |
23 calls. Slowest were:
| 107.0ms | (* (* (* (/ 1 x) (/ 1 x)) (/ 1 x)) (* (* (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5)) (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5))) |
| 81.0ms | (+ (* 0.2514179000665375 (/ 1 (pow x 3))) (* 0.5 (/ 1 x))) |
| 75.0ms | (+ (log (/ 1 x)) (log (fma (/ 1 x) (/ 0.2514179000665375 x) 0.5))) |
2 alts after pruning (0 fresh and 2 done)
Merged error: 0b
100% (0.0b remaining)
Error of 0.2b against oracle of 0.2b and baseline of 31.5b