Details

simplify5.9m (33.4%)

sample5.3m (29.9%)

Algorithm
100×intervals
Results
1.7m15995×body10240exit
1.3m450567×body80valid
46.3s30807×body1280valid
24.6s23906×body640valid
10.9s3147×body1280nan
10.6s69909×body80nan
9.2s2713×body640nan
8.9s11925×body320valid
5.4s10136×body160valid
2.8s49536×pre80true
2.3s1261×body320nan
1.3s1359×body2560valid
789.0ms664×body160nan
582.0ms596×body5120valid

prune3.0m (17.0%)

rewrite1.9m (11.0%)

Algorithm
151×rewrite-expression-head
Calls

449 calls:

16.3s
(* (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x)
14.8s
(* (exp (- (- 0.5) (- z -6))) (* (+ (+ (+ (+ (/ -176.6150291621406 (- z -3)) (/ -1259.1392167224028 (+ 1 z))) (+ (/ 771.3234287776531 (+ z 2)) (+ 0.9999999999998099 (/ 676.5203681218851 z)))) (/ 12.507343278686905 (- z -4))) (+ (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (- z -7))) (/ -0.13857109526572012 (+ z 5)))) (* (pow (+ 0.5 (- z -6)) (- z (- 1 0.5))) (sqrt (* PI 2)))))
14.7s
(* (exp -6) (* (+ (+ (+ (+ (/ -176.6150291621406 (- z -3)) (/ -1259.1392167224028 (+ 1 z))) (+ (/ 771.3234287776531 (+ z 2)) (+ 0.9999999999998099 (/ 676.5203681218851 z)))) (/ 12.507343278686905 (- z -4))) (+ (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (- z -7))) (/ -0.13857109526572012 (+ z 5)))) (* (pow (+ 0.5 (- z -6)) (- z (- 1 0.5))) (sqrt (* PI 2)))))
9.3s
(/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))
8.9s
(sqrt (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x))))
Rules
56236×*-un-lft-identity
53419×times-frac
23984×add-sqr-sqrt
21429×distribute-lft-out
18997×add-cube-cbrt
17566×sqrt-prod
3613×pow1
3369×add-exp-log
2910×associate-*l/
2741×associate-/r*
2739×frac-times
2299×add-cbrt-cube
2070×cbrt-div
2062×frac-add
2023×flip-+ flip3-+
1736×difference-of-squares
1640×associate-/r/
1423×pow-prod-up
1391×prod-exp
1172×associate-*l*
1109×flip3-- flip--
1095×associate-*r/
1059×prod-diff
1020×div-inv
1009×distribute-lft-out--
983×associate-*r*
950×div-exp
818×cbrt-unprod
679×associate-/l*
652×pow-prod-down
634×add-log-exp
607×pow-plus
596×sqrt-div
497×cbrt-undiv
450×expm1-log1p-u
449×insert-posit16 log1p-expm1-u
317×unswap-sqr
307×log-pow
279×cbrt-prod
265×pow-sqr
253×pow-sub
234×fma-def
229×swap-sqr
211×exp-sum
196×1-exp
195×pow2
188×rec-exp
174×fma-neg
169×exp-diff
165×pow1/3
153×distribute-rgt-in distribute-lft-in
147×sqrt-pow1
142×*-commutative
140×log-prod
92×frac-2neg clear-num unpow-prod-down
88×unpow-prod-up
84×pow-flip
82×associate--l+
78×inv-pow
72×sub-neg
69×cube-unmult associate-/l/
68×pow1/2
66×fma-udef
60×pow-exp
52×sum-log
51×diff-log
48×associate-+r+
42×sqr-pow
33×pow-to-exp
31×+-commutative
22×frac-sub rem-sqrt-square distribute-neg-in
21×un-div-inv
16×pow-unpow exp-prod
15×distribute-lft-neg-in distribute-rgt-neg-in distribute-rgt-out
14×distribute-rgt1-in pow3 hypot-udef
12×associate--r+
11×hypot-def
pow-pow
div-sub
difference-of-sqr-1 associate--l- expm1-def rem-log-exp
associate-+l+ cos-sum
neg-mul-1 associate--r- log-div
sin-sum
sqrt-unprod rem-square-sqrt
neg-sub0 expm1-log1p expm1-udef associate-+r- rem-exp-log exp-neg
remove-posit16 remove-double-div log1p-expm1 log1p-udef associate-+l- e-exp-1

series1.2m (7.0%)

Calls

449 calls:

1.3s
(pow (* (pow -2 1.0) (pow (log u1) 1.0)) 0.5)
1.1s
(/ (pow (* -2 (log u1)) 0.5) 6)
1.1s
(* (/ 1 6) (pow (* -2 (log u1)) 0.5))
1.1s
(* (/ (/ (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma 2 i (+ alpha beta))) (+ (fma 2 i (+ alpha beta)) (sqrt 1.0))) (/ (* (+ (+ alpha beta) i) i) (fma 2 i (+ alpha beta))))
986.0ms
(* 1/6 (pow (* (pow -2 1.0) (pow (log u1) 1.0)) 0.5))

regimes8.9s (0.8%)

Accuracy

Total 18.9b remaining (12.5%)

Threshold costs 6.6b (4.4%)

3.3b73.7%Octave 3.8, jcobi/1
3.0b-9.7%_divideComplex, imaginary part
3.0b0%_divideComplex, real part
1.5b46.2%Octave 3.8, jcobi/4
1.5b0%Jmat.Real.lambertw, newton loop step

localize6.9s (0.6%)

bsearch1.1s (0.1%)

end1.0ms (0.0%)