Details

rewrite13.2m (48.6%)

Algorithm
158×rewrite-expression-head
Calls

460 calls:

9.6m
(/ (* (* (sqrt (* PI 2.0)) (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (* 1 z))) (/ -1259.1392167224028 (+ (- z 1.0) 2.0))) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (/ -176.6150291621406 (+ (- z 1.0) 4.0))) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (/ -0.13857109526572012 (+ (- z 1.0) 6.0))) (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0))) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)))) (exp (+ (+ (- z 1.0) 7.0) 0.5)))
50.2s
(* (/ (sqrt 1) (* (cbrt (/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))))) (cbrt (/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))))))) (* (/ (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))) (cbrt (/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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)))))))) x))
19.3s
(* (/ (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))) (/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))))) x)
18.1s
(/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))))
18.0s
(/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))))
Rules
23564×times-frac
15415×*-un-lft-identity
14230×add-sqr-sqrt
10346×add-cube-cbrt
5903×sqrt-prod
2955×add-exp-log
2436×add-cbrt-cube
1780×distribute-lft-out
1678×associate-*r*
1541×difference-of-squares
1424×pow1
1153×prod-exp
1128×associate-*l*
998×associate-/r/
908×distribute-lft-out--
892×cbrt-unprod
726×div-inv
710×add-log-exp
704×flip-+ flip3-+
677×associate-/r*
667×div-exp
639×unpow-prod-down
626×associate-/l*
513×cbrt-undiv
511×associate-*r/
498×flip3-- flip--
406×frac-times
355×pow-prod-down
332×cbrt-prod
320×associate-*l/
298×unswap-sqr
232×sqrt-div
206×frac-sub
201×associate-/l/
194×log-pow
173×swap-sqr
149×cube-prod
148×cbrt-div
135×*-commutative
134×log-prod
133×pow-exp
126×sqrt-pow1
114×sqr-pow
101×distribute-rgt-in distribute-lft-in
100×frac-add pow-to-exp
98×pow1/3
90×frac-2neg clear-num
88×pow-unpow
87×pow1/2
84×sub-neg
74×diff-log
68×sum-log
61×rem-sqrt-square
60×distribute-neg-frac
40×cube-div +-commutative associate--l+
36×pow-prod-up 1-exp
34×pow-sqr pow-pow
33×rec-exp
32×distribute-rgt-neg-in
31×rem-cube-cbrt
30×distribute-lft-neg-in
25×pow-flip
24×exp-prod
23×pow2 pow-plus
22×unpow3 cube-mult
19×log-div
14×exp-sum associate--r+
12×associate-+l+
11×div-sub
10×inv-pow associate-+r+ neg-mul-1
9×rem-log-exp
8×associate--l-
7×exp-diff associate-+l- unpow2
6×difference-cubes un-div-inv rem-exp-log distribute-neg-in sum-cubes
5×cos-sum
4×sin-sum
3×distribute-rgt-neg-out
2×associate--r-
1×unpow1/3 sqrt-unprod associate-+r- unpow-prod-up pow-sub rem-square-sqrt

series10.6m (39.0%)

Calls

457 calls:

6.1m
(+ (pow (pow (/ (/ (+ alpha beta) (/ (+ (+ alpha beta) (* 2.0 i)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 3) 3) (pow (pow 1.0 3) 3))
2.1m
(pow (/ (/ (+ alpha beta) (/ (+ (+ alpha beta) (* 2.0 i)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 6)
59.5s
(- (pow 1.0 3) (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (/ (+ (pow 0.254829592 3) (pow (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))) 3)) (+ (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))) (- (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))) 0.254829592)) (* 0.254829592 0.254829592)))) (exp (- (* (fabs x) (fabs x))))) 3))
13.9s
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (/ (+ (pow 0.254829592 3) (pow (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))) 3)) (+ (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))) (- (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))) 0.254829592)) (* 0.254829592 0.254829592)))) (exp (- (* (fabs x) (fabs x))))))
4.1s
(+ (pow (/ (/ (+ alpha beta) (/ (+ (+ alpha beta) (* 2.0 i)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 3) (pow 1.0 3))

sample1.7m (6.3%)

Algorithm
99×intervals
Results
42.8s11830×body10240exit
21.1s482608×body80valid
11.3s27380×body1280valid
5.1s20759×body640valid
2.6s69699×body80nan
1.7s10743×body320valid
1.6s1552×body1280nan
1.0s9447×body160valid
962.0ms1326×body640nan
642.0ms49536×pre80true
389.0ms1536×body2560valid
319.0ms616×body320nan
146.0ms598×body5120valid
111.0ms322×body160nan

prune1.2m (4.3%)

simplify16.7s (1.0%)

Algorithm
256×egg-herbie

regimes6.7s (0.4%)

Accuracy

Total 88.7b remaining (31.4%)

Threshold costs 5.5b (2%)

10.8b56.5%math.log/2 on complex, real part
10.5b57.8%math.log10 on complex, real part
10.4b58.4%math.log/1 on complex, real part
9.8b59.6%math.abs on complex
9.6b54.6%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.6s (0.2%)

bsearch1.6s (0.1%)

end1.0ms (0.0%)