Details

simplify8.4m (35.2%)

rewrite7.0m (29.3%)

Algorithm
158×rewrite-expression-head
Calls

473 calls:

1.4m
(+ (+ (+ (/ 1.5056327351493116e-07 (- 8 z)) (/ 9.984369578019572e-06 (- 7 z))) (+ (+ (/ -0.13857109526572012 (- 6 z)) (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (/ -176.6150291621406 (+ 1 (+ (- 1 z) 2))))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1 z)) (/ -1259.1392167224028 (- 2 z)))))) (/ 12.507343278686905 (- 5 z)))
1.2m
(+ (+ (+ (/ 1.5056327351493116e-07 (- 8 z)) (/ 9.984369578019572e-06 (- 7 z))) (+ (+ (/ -0.13857109526572012 (- 6 z)) (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (/ -176.6150291621406 (+ 1 (+ (- 1 z) 2))))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1 z)) (/ -1259.1392167224028 (- 2 z)))))) (/ 12.507343278686905 (- 5 z)))
56.3s
(* (* (* (/ (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))
27.5s
(* (* (/ (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))))
25.5s
(* (+ (+ (+ (+ (/ 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))))
Rules
13359×times-frac
12598×*-un-lft-identity
9065×add-sqr-sqrt
7410×frac-add
7151×add-cube-cbrt
3060×add-exp-log
2880×distribute-lft-out
2569×flip-+ flip3-+
2534×pow1
2080×add-cbrt-cube
1564×sqrt-prod
1460×frac-times
1338×associate-*l/
1330×associate-*l*
1273×associate-*r*
1246×prod-diff
1213×prod-exp
1179×associate-*r/
857×add-log-exp
844×associate-/l*
769×div-exp
680×pow-sub
673×associate-/r*
667×cbrt-unprod
665×associate-/r/
533×pow-prod-up
479×difference-of-squares
478×div-inv
475×pow-prod-down
474×expm1-log1p-u log1p-expm1-u
473×insert-posit16
397×log-pow
396×cbrt-undiv
320×distribute-lft-out--
313×exp-sum
298×unswap-sqr
264×pow-plus
248×sum-log
242×fma-def
238×log-prod
231×associate-/l/
198×unpow-prod-down
182×1-exp rec-exp
181×distribute-rgt-in distribute-lft-in
168×cbrt-prod
163×sqrt-pow1
159×fma-neg
157×*-commutative
156×unpow-prod-up
152×swap-sqr
135×pow-sqr
105×sqrt-div
103×flip3-- flip--
99×pow-exp
98×pow2
95×fma-udef distribute-rgt-out
94×pow1/2
93×frac-2neg clear-num
92×cbrt-div
91×sqr-pow
84×pow-flip
78×inv-pow
77×pow1/3
76×associate--r+
67×sub-neg
58×distribute-lft1-in
57×cube-prod pow-to-exp
56×distribute-rgt1-in
38×frac-sub
36×rem-sqrt-square
32×+-commutative
31×diff-log
30×pow-unpow cube-mult
29×unpow3
25×cube-unmult exp-prod
24×rem-cube-cbrt pow-div
20×pow-pow un-div-inv
19×count-2
16×associate-+l+
15×sin-sum pow3 div-sub cos-sum
14×hypot-def
12×distribute-lft-neg-in distribute-rgt-neg-in rem-log-exp
10×exp-diff
associate-+r- associate-+r+
cube-div distribute-neg-in
associate--l+ associate--r-
sub-div neg-mul-1
distribute-rgt-neg-out sqrt-unprod associate-+l- rem-cbrt-cube log-div
expm1-log1p expm1-udef log1p-expm1 sub0-neg rem-square-sqrt log1p-udef
distribute-lft-neg-out difference-cubes exp-to-pow distribute-rgt-out-- rem-exp-log unsub-neg

sample3.7m (15.4%)

Algorithm
100×intervals
Results
1.1m502590×body80valid
1.0m15932×body10240exit
31.9s30793×body1280valid
17.4s23526×body640valid
7.7s68774×body80nan
5.6s12053×body320valid
5.5s3121×body1280nan
3.4s10184×body160valid
3.2s2578×body640nan
2.0s49536×pre80true
1.3s1260×body320nan
1.2s1441×body2560valid
370.0ms585×body5120valid
314.0ms679×body160nan

prune3.0m (12.6%)

series1.5m (6.2%)

Calls

473 calls:

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

regimes10.0s (0.7%)

Accuracy

Total 26.5b remaining (19.6%)

Threshold costs 7.2b (5.3%)

9.2b43.3%Octave 3.8, jcobi/4
3.0b75.9%Octave 3.8, jcobi/1
3.0b0%_divideComplex, imaginary part
2.4b0%_divideComplex, real part
1.7b0%Jmat.Real.lambertw, newton loop step

localize7.6s (0.5%)

bsearch857.0ms (0.1%)

end1.0ms (0.0%)