Details

rewrite11.1m (72.0%)

Algorithm
149×rewrite-expression-head
Calls

430 calls:

15.0s
(* (/ (+ (+ (+ (+ (+ 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)))) (+ (+ (+ (+ (+ (+ 1.0 (* 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 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x)
8.8s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
8.4s
(/ (+ (+ (+ (+ (+ 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)))) (+ (+ (+ (+ (+ (+ 1.0 (* 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 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))
4.4s
(* (/ (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))))) 1) (/ (* 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)))))) (fma (* 2.0 0.0001789971) (* (pow x 3) (* (pow x 3) (* (pow x 2) (* x (pow x 3))))) (fma 0.0008327945 (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3))))) (fma 0.0140005442 (* (pow x 2) (* (pow x 2) (* x (pow x 3)))) (fma 0.0694555761 (pow x 6) (fma 0.2909738639 (pow x 4) (fma (* 0.7715471019 x) x 1.0))))))))
1.6s
(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
17307×times-frac
15326×*-un-lft-identity
10827×add-sqr-sqrt
8419×add-cube-cbrt
2779×add-exp-log
2593×distribute-lft-out
2031×prod-diff
1873×add-cbrt-cube
1538×fma-def
1531×pow1
1079×distribute-rgt-in distribute-lft-in
1003×prod-exp
897×distribute-lft-out--
894×difference-of-squares
811×associate-*l*
772×div-exp
732×sqrt-prod
725×associate-*r*
707×associate-/l*
627×add-log-exp
620×associate-/r/
613×associate-/r*
609×cbrt-unprod
520×div-inv
448×cbrt-prod
438×pow-prod-down
433×log1p-expm1-u
430×expm1-log1p-u
406×cbrt-undiv
374×flip-+ flip3-+
290×unpow-prod-down
239×log-pow
234×unswap-sqr
225×associate-*r/
189×1-exp
187×fma-neg
176×sqrt-pow1
168×flip3-- flip--
158×associate-*l/
155×associate--l+
147×sqrt-div
132×cube-prod
121×*-commutative
119×exp-sum
115×cbrt-div
106×pow-exp
103×log-prod
100×sub-neg
99×rec-exp swap-sqr
97×frac-times
95×pow1/3
85×pow1/2
79×pow-prod-up
78×frac-2neg clear-num
75×diff-log sqr-pow
74×pow-unpow
73×associate-/l/
71×associate--r+
66×pow-to-exp
61×exp-prod
55×pow-sqr
52×fma-udef
50×frac-add
36×rem-sqrt-square distribute-lft-neg-in distribute-rgt-neg-in
34×sum-log
33×inv-pow pow-flip
32×pow-plus
31×associate-+l+
29×neg-mul-1
27×pow2
19×pow-pow +-commutative
18×unpow3 cube-mult
17×hypot-def
15×associate--l-
14×cube-div
13×rem-cube-cbrt rem-exp-log
12×hypot-udef div-sub rem-log-exp
un-div-inv
neg-log
rem-cbrt-cube log-div unpow2
sin-sum frac-sub associate-+l-
exp-diff log1p-udef sqrt-unprod log1p-expm1 expm1-log1p rem-square-sqrt neg-sub0 log-rec expm1-udef
sqrt-undiv associate-+r+ distribute-rgt-neg-out
exp-sqrt sqrt-pow2 distribute-lft-neg-out exp-neg distribute-neg-in unpow-prod-up pow-sub distribute-frac-neg cos-sum associate--r- exp-to-pow mul-1-neg

sample1.8m (11.6%)

Algorithm
99×intervals
Results
42.5s12064×body10240exit
21.5s495844×body80valid
11.2s27503×body1280valid
5.0s21253×body640valid
2.4s68070×body80nan
1.7s1658×body1280nan
1.6s10646×body320valid
977.0ms9277×body160valid
971.0ms1391×body640nan
585.0ms49536×pre80true
378.0ms1552×body2560valid
329.0ms646×body320nan
146.0ms587×body5120valid
99.0ms301×body160nan

series1.1m (7.0%)

Calls

430 calls:

3.1s
(pow (* (/ (+ alpha beta) (+ (+ alpha beta) (fma 2.0 i 2.0))) (/ (- beta alpha) (+ (fma 2.0 i beta) alpha))) 3)
2.3s
(/ (* (* (/ i (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i))))) (/ (fma i 2.0 (+ alpha beta)) (/ (sqrt (fma beta alpha (* i (+ (+ alpha beta) i)))) (fma i 2.0 (+ alpha beta)))))
2.2s
(/ (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta))))
2.0s
(/ (* (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (fma i 2.0 (+ alpha beta))) (fma beta alpha (* i (+ (+ alpha beta) i))))
858.0ms
(/ (* i (+ (+ alpha beta) i)) (/ (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta)))))

prune1.1m (7.0%)

simplify13.7s (1.5%)

regimes5.1s (0.6%)

Accuracy

Total 44.0b remaining (35.1%)

Threshold costs 6.3b (5%)

9.0b26.4%_divideComplex, imaginary part
8.2b24.7%_divideComplex, real part
6.3b0%Octave 3.8, jcobi/2
4.7b19.5%math.sqrt on complex, real part
4.4b26.6%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.6s (0.3%)

bsearch751.0ms (0.1%)

end1.0ms (0.0%)