Details

rewrite9.0m (34.7%)

Algorithm
157×rewrite-expression-head
Calls

482 calls:

1.9m
(* (* (* (* (/ (sqrt (* PI 2)) (exp (+ 0.5 z))) (exp -6)) (+ (+ (+ (+ (/ 771.3234287776531 (+ z 2)) (/ -1259.1392167224028 (+ 1 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -176.6150291621406 (+ 3 z)))) (+ (+ (/ -0.13857109526572012 (+ z 5)) (/ 12.507343278686905 (- z -4))) (/ 1.5056327351493116e-07 (+ 7 z)))) (/ 9.984369578019572e-06 (+ 6 z)))) (pow (+ 0.5 (- z -6)) (- z 1))) (pow (+ 0.5 (- z -6)) 0.5))
1.9m
(* (* (* (/ (sqrt (* PI 2)) (exp (+ 0.5 (- z -6)))) (+ (+ (+ (+ (/ 771.3234287776531 (+ z 2)) (/ -1259.1392167224028 (+ 1 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -176.6150291621406 (+ 3 z)))) (+ (+ (/ -0.13857109526572012 (+ z 5)) (/ 12.507343278686905 (- z -4))) (/ 1.5056327351493116e-07 (+ 7 z)))) (/ 9.984369578019572e-06 (+ 6 z)))) (pow (+ 0.5 (- z -6)) (- z 1))) (pow (+ 0.5 (- z -6)) 0.5))
57.8s
(* (* (/ (sqrt (* PI 2)) (exp (+ 0.5 (- z -6)))) (+ (+ (+ (+ (/ 771.3234287776531 (+ z 2)) (/ -1259.1392167224028 (+ 1 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -176.6150291621406 (+ 3 z)))) (+ (+ (/ -0.13857109526572012 (+ z 5)) (/ 12.507343278686905 (- z -4))) (/ 1.5056327351493116e-07 (+ 7 z)))) (/ 9.984369578019572e-06 (+ 6 z)))) (pow (+ 0.5 (- z -6)) (- z (- 1 0.5))))
50.1s
(* (/ (sqrt (* PI 2)) (exp (+ 0.5 (- z -6)))) (+ (+ (+ (+ (/ 771.3234287776531 (+ z 2)) (/ -1259.1392167224028 (+ 1 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -176.6150291621406 (+ 3 z)))) (+ (+ (/ -0.13857109526572012 (+ z 5)) (/ 12.507343278686905 (- z -4))) (/ 1.5056327351493116e-07 (+ 7 z)))) (/ 9.984369578019572e-06 (+ 6 z))))
44.6s
(* (* d1 (* (* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1)) d1)) d1)
Rules
7450×frac-add
7270×*-un-lft-identity
6006×times-frac
5944×add-sqr-sqrt
5281×frac-times
4718×pow1
4237×add-cube-cbrt
3647×associate-*l/
3181×add-exp-log
2586×flip-+ flip3-+
2384×add-cbrt-cube
2316×pow-prod-up
1928×sqrt-prod
1914×associate-*r/
1580×cbrt-div
1479×distribute-lft-out
1451×prod-exp
1096×pow-plus
1084×associate-*r*
1079×associate-*l*
971×prod-diff
946×cbrt-unprod
921×add-log-exp
643×pow-prod-down
619×un-div-inv
606×div-exp
568×cbrt-prod
516×associate-/l*
486×log1p-expm1-u
485×expm1-log1p-u
482×insert-posit16
453×pow-sub
392×cbrt-undiv
389×associate-/r*
378×pow-sqr
347×unswap-sqr
333×difference-of-squares
321×associate-/r/
311×distribute-lft-out--
300×pow2
294×div-inv
268×sum-log
220×exp-sum
219×fma-def
173×flip3-- flip--
172×swap-sqr
165×unpow-prod-down
162×sqrt-pow1
153×*-commutative
149×fma-neg
145×distribute-rgt-in distribute-lft-in
143×log-pow
111×sqrt-div
104×1-exp
96×pow1/2
95×pow-exp
94×rec-exp
92×pow1/3 frac-2neg clear-num
88×cube-unmult
87×fma-udef
79×associate--r+
77×pow-to-exp
71×sqr-pow
69×associate-/l/
61×log-prod
53×sub-neg
45×distribute-rgt-out
44×pow-flip
42×diff-log
38×inv-pow
37×rem-sqrt-square
36×+-commutative
30×distribute-rgt1-in
28×frac-sub pow-unpow
26×unpow-prod-up
18×cos-sum exp-prod
17×exp-diff div-sub
16×associate-+r- hypot-def
12×distribute-lft-neg-in distribute-rgt-neg-in
11×associate-+l+ pow-div
10×pow3 pow-pow count-2
sin-sum
expm1-log1p expm1-udef cube-prod log1p-expm1 distribute-neg-in log1p-udef
associate--l+ sqrt-unprod
neg-mul-1 associate-+l-
sqrt-undiv associate-+r+ distribute-rgt-neg-out hypot-udef rem-square-sqrt unpow1/3 associate--r-
unpow3 sub0-neg cube-mult rem-cbrt-cube rem-log-exp
distribute-lft1-in distribute-lft-neg-out associate--l- sqrt-pow2 remove-posit16 unsub-neg log-div

simplify8.6m (32.9%)

sample3.9m (14.8%)

Algorithm
100×intervals
Results
1.2m15594×body10240exit
1.1m493302×body80valid
30.2s29992×body1280valid
15.7s22854×body640valid
7.3s71455×body80nan
5.7s11914×body320valid
5.2s3001×body1280nan
3.5s10239×body160valid
3.4s2456×body640nan
2.2s49536×pre80true
1.3s1509×body2560valid
976.0ms1237×body320nan
426.0ms634×body5120valid
350.0ms701×body160nan

prune3.0m (11.4%)

series1.3m (5.0%)

Calls

482 calls:

1.3s
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))
1.2s
(* (pow (* -2 (log u1)) 0.5) (sqrt 1/6))
1.1s
(- 1 (* (/ (/ 1 (fma 0.3275911 (fabs x) 1)) (exp (* (fabs x) (fabs x)))) (+ (/ (/ 1.421413741 (fma 0.3275911 (fabs x) 1)) (fma 0.3275911 (fabs x) 1)) (- (+ 0.254829592 (/ 1.061405429 (* (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1)) (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1))))) (+ (/ 0.284496736 (fma 0.3275911 (fabs x) 1)) (/ (/ 1.453152027 (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1))) (fma 0.3275911 (fabs x) 1)))))))
1.1s
(* (* (pow (* -2 (log u1)) 0.5) (sqrt 1/6)) (sqrt 1/6))
1.1s
(/ (pow (* -2 (log u1)) (/ 0.5 2)) 6)

regimes9.7s (0.6%)

Accuracy

Total 16.5b remaining (12.3%)

Threshold costs 6.7b (5%)

2.6b0%_divideComplex, real part
2.6b61.9%Octave 3.8, jcobi/4
2.4b0%_divideComplex, imaginary part
1.2b51.3%Octave 3.8, jcobi/3
1.0b0%Jmat.Real.lambertw, newton loop step

localize7.4s (0.5%)

bsearch1.3s (0.1%)

end1.0ms (0.0%)