Details

rewrite20.1m (46.7%)

Algorithm
154×rewrite-expression-head
Calls

449 calls:

2.4m
(+ (* (* (- (* (/ 771.3234287776531 (+ 3.0 (- z))) (/ 771.3234287776531 (+ 3.0 (- z)))) (* (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))))) (- (* (/ -1259.1392167224028 (+ (- z) 2.0)) (/ -1259.1392167224028 (+ (- z) 2.0))) (* (/ -176.6150291621406 (+ (- z) 4.0)) (/ -176.6150291621406 (+ (- z) 4.0))))) (* (* (+ (- z) 2.0) (+ (- z) 2.0)) (+ (- z) 4.0))) (* (+ (/ -1259.1392167224028 (+ (- z) 2.0)) (/ -176.6150291621406 (+ (- z) 4.0))) (* (- (/ 771.3234287776531 (+ 3.0 (- z))) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))) (- (* (* -1259.1392167224028 -1259.1392167224028) (+ (- z) 4.0)) (* (* (+ (- z) 2.0) (+ (- z) 2.0)) (* (/ -176.6150291621406 (+ (- z) 4.0)) -176.6150291621406))))))
1.5m
(* (* (/ 1.0 (sqrt PI)) (exp (* (fabs x) (fabs x)))) (+ (+ (+ (/ 1.0 (fabs x)) (* (/ 1.0 2.0) (* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))))) (* (/ 3.0 4.0) (* (* (* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))))) (* (/ 15.0 8.0) (* (* (* (* (* (* (/ 1.0 (fabs x)) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))))))
1.0m
(/ (+ (/ 9.984369578019572e-06 (+ 7.0 (- z))) (+ (+ (+ (/ 771.3234287776531 (+ 3.0 (- z))) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))) (+ (/ -1259.1392167224028 (+ (- z) 2.0)) (/ -176.6150291621406 (+ (- z) 4.0)))) (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ 12.507343278686905 (+ 5.0 (- z)))) (/ 1.5056327351493116e-07 (+ 8.0 (- z)))))) (exp (+ (+ 0.5 7.0) (- z))))
20.1s
(* (/ (* (cbrt (+ (+ (+ (+ (+ 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 (+ (+ (+ (+ (+ 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) (* (/ (cbrt (+ (+ (+ (+ (+ 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))
19.2s
(* (/ (+ (+ (+ (+ (+ 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)
Rules
3605908×frac-times
2873924×associate-*l/
1771191×associate-*r/
1582267×frac-add
38683×flip-+ flip3-+
15953×times-frac
14664×frac-sub
14360×*-un-lft-identity
10600×add-sqr-sqrt
8864×add-cube-cbrt
3820×flip3-- flip--
3648×sqrt-prod
3318×add-exp-log
3099×pow1
2711×distribute-lft-out
2354×prod-diff
2126×add-cbrt-cube
1389×associate-*l*
1242×prod-exp
1142×pow-prod-up
1131×associate-*r*
1069×associate-+l+
993×div-exp
927×distribute-rgt-in distribute-lft-in
735×add-log-exp
711×cbrt-unprod
689×associate-/r*
667×associate-/l*
637×div-inv
603×pow-plus
522×pow-prod-down
485×associate-/r/
470×cbrt-undiv
444×log1p-expm1-u
443×expm1-log1p-u
441×difference-of-squares
388×cbrt-prod
324×distribute-lft-out--
303×distribute-rgt-neg-in
298×unswap-sqr
238×pow-sqr
193×fma-neg
192×swap-sqr
191×pow2
154×1-exp
153×associate-/l/
148×exp-sum rec-exp
145×log-pow
144×*-commutative
130×sum-log
129×cbrt-div
118×fma-def
110×sqrt-pow1
108×pow1/3
99×pow-exp
92×associate-+r+
84×frac-2neg clear-num
82×sqrt-div sub-neg
76×pow1/2
73×fma-udef
72×associate--r+
66×associate--l+
62×log-prod
61×sqr-pow
58×unpow-prod-down diff-log
54×cube-unmult
51×pow-to-exp
31×rem-sqrt-square
29×distribute-rgt1-in distribute-rgt-out
24×exp-prod
23×+-commutative
18×un-div-inv
17×div-sub
16×pow3
10×pow-flip cos-sum
distribute-lft-neg-in pow-unpow hypot-def
inv-pow
cube-prod associate--l- pow-pow rem-exp-log
sqrt-unprod associate-+l-
exp-diff sin-sum neg-log
rem-square-sqrt
neg-sub0 rem-log-exp unpow3 cube-mult neg-mul-1 exp-to-pow rem-cbrt-cube associate--r-
unpow-prod-up expm1-log1p expm1-udef associate-+r- distribute-lft1-in sqrt-undiv distribute-rgt-neg-out log1p-expm1 exp-neg distribute-rgt-out-- distribute-neg-in unsub-neg pow-sub log1p-udef

simplify12.3m (28.6%)

prune6.6m (15.5%)

sample2.1m (4.9%)

Algorithm
97×intervals
Results
45.6s429950×body80valid
26.6s27103×body1280valid
13.1s21421×body640valid
6.9s8255×body10240exit
6.2s68393×body80nan
3.7s9536×body160valid
3.6s10854×body320valid
1.7s49536×pre80true
740.0ms626×body5120valid
719.0ms1387×body2560valid
166.0ms108×body1280nan
109.0ms85×body640nan
42.0ms43×body320nan
11.0ms20×body160nan

series1.5m (3.5%)

Calls

443 calls:

3.6s
(/ (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))
3.4s
(* (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))) (/ 1 (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))))
3.2s
(* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))))
1.3s
(* (* (/ 1 (cbrt 6.0)) (/ 1 (cbrt 6.0))) (/ (pow (* -2.0 (log u1)) 0.5) (/ (cbrt 6.0) 1.0)))
1.2s
(* (* (pow (* -2.0 (log u1)) 0.5) (sqrt (/ 1 (/ 6.0 1.0)))) (sqrt (/ 1 (/ 6.0 1.0))))

regimes10.7s (0.4%)

Accuracy

Total 30.3b remaining (21.6%)

Threshold costs 5.6b (4%)

9.6b18.3%_divideComplex, real part
7.7b30.8%_divideComplex, imaginary part
3.1b75.7%Octave 3.8, jcobi/1
2.4b34.1%Octave 3.8, jcobi/4
2.0b20.6%Octave 3.8, jcobi/3

localize8.9s (0.3%)

bsearch1.1s (0.0%)

end1.0ms (0.0%)