Details

rewrite11.1m (71.5%)

Algorithm
149×rewrite-expression-head
Calls

430 calls:

15.9s
(* (/ (+ (+ (+ (+ (+ 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)
9.2s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
8.9s
(/ (+ (+ (+ (+ (+ 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.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))))) 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.7s
(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
17444×times-frac
15001×*-un-lft-identity
10507×add-sqr-sqrt
8097×add-cube-cbrt
3023×add-exp-log
2595×distribute-lft-out
1960×add-cbrt-cube
1537×fma-def
1526×pow1
1346×prod-diff
1065×prod-exp
1027×sqrt-prod
957×div-exp
899×difference-of-squares
898×distribute-lft-out--
800×associate-*l*
779×associate-/l*
708×associate-*r*
649×associate-/r*
636×associate-/r/
620×add-log-exp
598×cbrt-unprod
485×div-inv
451×cbrt-undiv
434×log1p-expm1-u
430×expm1-log1p-u
425×pow-prod-down
397×cbrt-prod
381×flip-+ flip3-+
345×unpow-prod-down
273×1-exp
250×associate-*r/
242×log-pow
229×unswap-sqr
225×sqrt-pow1
218×distribute-rgt-in distribute-lft-in
212×sqrt-div
193×fma-neg
175×flip3-- flip--
156×associate-*l/
155×associate--l+
136×cube-prod
129×rec-exp
125×exp-sum
124×sqr-pow
123×pow1/2
116×*-commutative
115×cbrt-div
112×log-prod
109×frac-times
96×pow-prod-up
88×swap-sqr
84×pow-exp
83×frac-2neg clear-num
81×pow1/3 sub-neg
75×associate-/l/
74×diff-log pow-unpow
71×exp-prod associate--r+
64×pow-sqr
55×pow-to-exp
52×fma-udef
50×frac-add
38×rem-sqrt-square
36×distribute-lft-neg-in distribute-rgt-neg-in
35×pow-plus
33×pow-flip
32×inv-pow sum-log
31×associate-+l+
29×neg-mul-1
28×pow2
19×pow-pow unpow3 cube-mult
18×+-commutative
17×hypot-def
15×associate--l-
14×cube-div
13×rem-cube-cbrt
12×rem-exp-log hypot-udef div-sub rem-log-exp
10×un-div-inv sqrt-unprod
neg-log
rem-cbrt-cube log-div unpow2
sin-sum expm1-log1p rem-square-sqrt frac-sub expm1-udef associate-+l-
log1p-udef log1p-expm1 neg-sub0 log-rec
exp-diff associate-+r+ distribute-rgt-neg-out
sqrt-undiv sqrt-pow2 exp-neg distribute-neg-in unpow-prod-up pow-sub distribute-frac-neg cos-sum associate--r-

sample1.8m (11.5%)

Algorithm
99×intervals
Results
42.9s12064×body10240exit
21.8s495844×body80valid
11.2s27503×body1280valid
5.1s21253×body640valid
2.5s68070×body80nan
1.8s1658×body1280nan
1.6s10646×body320valid
1.0s1391×body640nan
983.0ms9277×body160valid
616.0ms49536×pre80true
373.0ms1552×body2560valid
337.0ms646×body320nan
136.0ms587×body5120valid
105.0ms301×body160nan

series1.2m (7.6%)

Calls

430 calls:

4.9s
(fma (* (cbrt (fma (sqrt 1.0) (sqrt 1.0) (- (* (cbrt z) (* (cbrt z) (cbrt z)))))) (cbrt (fma (sqrt 1.0) (sqrt 1.0) (- (* (cbrt z) (* (cbrt z) (cbrt z))))))) (cbrt (fma (sqrt 1.0) (sqrt 1.0) (- (* (cbrt z) (* (cbrt z) (cbrt z)))))) (fma -1 z (- z 1.0)))
3.3s
(pow (* (/ (+ alpha beta) (+ (+ alpha beta) (fma 2.0 i 2.0))) (/ (- beta alpha) (+ (fma 2.0 i beta) alpha))) 3)
2.4s
(/ (* (* (/ 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.1s
(/ (* (- (* (+ (+ 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.1s
(/ (* (* (- (* (+ (+ 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))))

prune1.1m (7.0%)

simplify13.8s (1.5%)

regimes5.2s (0.6%)

Accuracy

Total 43.3b remaining (34.6%)

Threshold costs 5.6b (4.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.5%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.6s (0.3%)

bsearch725.0ms (0.1%)

end1.0ms (0.0%)