Details

rewrite12.2m (49.0%)

Algorithm
158×rewrite-expression-head
Calls

478 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.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)))))
6.6s
(/ (/ 1 (/ (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)) 1)) (/ (/ (fma i 2.0 (+ alpha beta)) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta)))) (/ i (/ (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)) (+ (+ alpha beta) i)))))
1.2s
(* (/ (* (/ (+ alpha beta) (fma i 2.0 (+ alpha beta))) (/ (- beta alpha) (- (+ alpha beta) (* 2.0 i)))) (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (/ (/ (- (+ alpha beta) (* 2.0 i)) (sqrt (* (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))))) (sqrt (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)))))
1.1s
(- (pow (pow (/ alpha (+ (+ alpha beta) 2.0)) 3) 3) (pow (pow 1.0 3) 3))
Rules
1880117×times-frac
1004090×*-un-lft-identity
762902×add-sqr-sqrt
615622×add-cube-cbrt
180215×distribute-lft-out
73684×difference-of-squares
55778×div-inv
39433×sqrt-prod
36497×distribute-lft-out--
29475×associate-/r/
7900×flip-+ flip3-+
5802×unpow-prod-down
5709×associate-/r*
5620×cube-prod
5057×prod-diff
3111×add-exp-log
2108×add-cbrt-cube
1555×pow1
1538×associate-*r*
1228×div-exp
1053×associate-/l*
835×associate-*l*
834×cbrt-prod
794×prod-exp
656×add-log-exp
641×sqr-pow
629×cbrt-undiv
499×unpow3 cube-mult
497×cbrt-unprod
478×expm1-log1p-u
477×log1p-expm1-u
365×pow-prod-down
311×fma-neg
274×sqrt-div
252×distribute-rgt-in unswap-sqr distribute-lft-in
235×sqrt-pow1
220×cbrt-div
204×log-pow
188×1-exp
166×associate-*r/
147×fma-def
144×pow1/2
138×exp-sum
133×flip3-- *-commutative flip--
131×pow-prod-up
129×swap-sqr
121×log-prod
113×frac-times
112×pow-exp
110×associate--l+
100×associate-*l/
96×frac-2neg clear-num
95×rec-exp
93×pow1/3
89×pow-unpow
79×sub-neg
73×distribute-rgt-neg-in distribute-lft-neg-in
71×associate--r+
70×pow-to-exp
69×associate-/l/
66×diff-log
65×pow-sqr
54×fma-udef
51×exp-prod
42×neg-mul-1
35×pow-pow pow-plus
33×rem-sqrt-square sum-log
29×cube-div
25×pow2
19×div-sub
18×rem-cube-cbrt hypot-udef
17×pow3 +-commutative hypot-def
16×unpow2
15×associate-+l+ pow-flip
14×sqrt-unprod
12×inv-pow
11×associate--l-
10×rem-log-exp
rem-exp-log
frac-sub
neg-log rem-cbrt-cube
rem-square-sqrt neg-sub0
sin-sum distribute-neg-in expm1-log1p expm1-udef
sqrt-pow2 log1p-udef log1p-expm1 distribute-frac-neg associate-+l- log-div
exp-diff difference-cubes distribute-lft-neg-out associate-+r- cos-sum distribute-rgt-neg-out exp-to-pow
exp-neg unpow-prod-up pow-sub unsub-neg associate--r- mul-1-neg

prune7.0m (28.1%)

simplify2.7m (11.0%)

sample1.8m (7.1%)

Algorithm
98×intervals
Results
45.1s12259×body10240exit
19.3s456095×body80valid
11.1s27591×body1280valid
4.9s20923×body640valid
2.6s70684×body80nan
1.7s1665×body1280nan
1.6s10607×body320valid
1.0s9407×body160valid
956.0ms1335×body640nan
493.0ms41536×pre80true
363.0ms683×body320nan
357.0ms1488×body2560valid
145.0ms607×body5120valid
109.0ms324×body160nan

series1.0m (4.2%)

Calls

476 calls:

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.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))))
1.0s
(/ (sqrt (fma 1.0 2.0 (+ alpha beta))) (* (cbrt (+ beta (+ alpha (fma alpha beta 1.0)))) (cbrt (+ beta (+ alpha (fma alpha beta 1.0))))))
859.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)))))
845.0ms
(/ (/ i (/ (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (+ (+ alpha beta) i))) (/ (fma i 2.0 (+ alpha beta)) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta)))))

regimes5.1s (0.3%)

Accuracy

Total 37.2b remaining (39.7%)

Threshold costs 5.6b (6%)

9.0b29.4%_divideComplex, imaginary part
7.4b37.6%_divideComplex, real part
5.5b13.3%Octave 3.8, jcobi/2
4.6b27.4%math.sqrt on complex, imaginary part, im greater than 0 branch
4.1b30.7%math.sqrt on complex, real part

localize2.5s (0.2%)

bsearch800.0ms (0.1%)

end1.0ms (0.0%)