Details

simplify8.7m (37.8%)

rewrite6.7m (29.2%)

Algorithm
159×rewrite-expression-head
Calls

469 calls:

1.5m
(* (/ (* (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)) (sqrt (* PI 2.0))) (exp (+ (+ (- z 1.0) 7.0) 0.5))) (+ (+ (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0))) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (+ (/ -176.6150291621406 (+ (- z 1.0) 4.0)) (+ (+ (/ -0.13857109526572012 (+ (- z 1.0) 6.0)) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (+ (/ 676.5203681218851 z) (+ (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) 0.9999999999998099))))))
1.4m
(/ (+ (/ -1259.1392167224028 (+ (- z) 2.0)) (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (/ -176.6150291621406 (+ 4.0 (- z)))) (/ 771.3234287776531 (+ 3.0 (- z)))) (+ (/ 12.507343278686905 (+ 5.0 (- z))) (+ (/ 1.5056327351493116e-07 (+ 8.0 (- z))) (+ (/ 9.984369578019572e-06 (+ (- z) 7.0)) (/ -0.13857109526572012 (+ 6.0 (- z)))))))) (exp (+ 0.5 (+ (- z) 7.0))))
45.4s
(+ (+ (+ (+ (/ 771.3234287776531 (- 3.0 z)) (+ (/ -176.6150291621406 (- 4.0 z)) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))))) (/ 12.507343278686905 (- 5.0 z))) (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))))) (/ -1259.1392167224028 (+ (- z) 2.0)))
27.5s
(- (log (+ (+ (+ (+ (/ 771.3234287776531 (- 3.0 z)) (+ (/ -176.6150291621406 (- 4.0 z)) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))))) (/ 12.507343278686905 (- 5.0 z))) (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))))) (/ -1259.1392167224028 (+ (- z) 2.0)))) (+ 0.5 (+ (- z) 7.0)))
25.4s
(* (/ (+ (+ (+ (+ (+ 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
10046×*-un-lft-identity
9572×times-frac
7683×add-sqr-sqrt
6041×add-cube-cbrt
2373×add-exp-log
2289×sqrt-prod
2135×distribute-lft-out
1848×add-cbrt-cube
1446×associate-*r*
1402×frac-add
1323×pow1
1243×prod-diff
929×associate-*l*
835×add-log-exp
828×prod-exp
641×associate-/l*
617×unpow-prod-down
595×flip-+ flip3-+
580×cbrt-unprod
562×cbrt-prod
522×difference-of-squares
501×associate-/r* div-exp
470×log1p-expm1-u expm1-log1p-u
392×div-inv
376×cube-prod
358×pow-prod-down
349×distribute-lft-out--
348×cbrt-undiv
321×associate-*r/
276×associate-/r/
249×unswap-sqr
242×distribute-rgt-in distribute-lft-in
212×log-pow sum-log
194×frac-times
172×fma-neg
163×sqr-pow
151×fma-def
138×exp-sum
136×*-commutative
135×associate-/l/
132×associate-*l/
126×sqrt-div
125×exp-prod
102×pow1/3 flip3-- flip--
95×clear-num frac-2neg
93×pow-exp
91×pow-unpow
76×swap-sqr log-prod
75×associate--l+
72×associate--r+
71×fma-udef
66×pow-to-exp
64×sqrt-pow1
62×cbrt-div sub-neg
61×1-exp rec-exp
54×unpow3 cube-mult
49×diff-log associate--l-
44×pow-prod-up log-div
42×pow1/2
37×distribute-rgt-neg-in
33×pow-pow
32×pow-sqr
28×distribute-lft-neg-in +-commutative
22×rem-cube-cbrt
21×rem-sqrt-square pow-plus
18×cube-div
17×div-sub
16×pow2
14×associate-+l-
12×hypot-udef
11×pow-flip
10×associate-+l+ un-div-inv associate-+r+
inv-pow
sin-sum unpow2 unpow-prod-up
neg-mul-1 rem-log-exp
pow-sub sqrt-undiv cos-sum
exp-diff frac-sub
expm1-log1p expm1-udef neg-sub0 hypot-def associate--r-
associate-+r- pow3 exp-neg distribute-neg-in sub-div rem-cbrt-cube distribute-rgt-neg-out
rem-exp-log sqrt-pow2 log1p-udef sqrt-unprod lft-mult-inverse log1p-expm1 rem-square-sqrt neg-log distribute-rgt-out--

sample3.7m (16.1%)

Algorithm
100×intervals
Results
1.0m15801×body10240exit
51.3s441895×body80valid
33.8s30693×body1280valid
15.7s23350×body640valid
6.9s3035×body1280nan
6.9s12132×body320valid
6.5s68545×body80nan
3.7s10283×body160valid
3.7s2596×body640nan
1.7s49536×pre80true
1.6s1262×body320nan
1.3s1451×body2560valid
396.0ms655×body160nan
371.0ms578×body5120valid

prune1.9m (8.2%)

series1.7m (7.4%)

Calls

469 calls:

3.0s
(* (/ i (* (+ (sqrt 1.0) (fma 2.0 i (+ alpha beta))) (/ (fma 2.0 i (+ alpha beta)) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i))))))) (/ (/ (+ (+ alpha beta) i) (/ (fma 2.0 i (+ alpha beta)) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i)))))) (- (fma 2.0 i (+ alpha beta)) (sqrt 1.0))))
1.7s
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ (fma -0.284496736 (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) 0.254829592) (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (- (* 1.0 1.0) (* (* 0.3275911 (fabs x)) (* 0.3275911 (fabs x))))) (* (- 1.0 (* 0.3275911 (fabs x))) 1.061405429)))))) (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x))))))) (exp (- (* (fabs x) (fabs x))))))
1.7s
(- (pow (* d1 (+ d4 (- d2 d3))) 3) (pow d1 6))
1.5s
(expm1 (log1p (* (/ 1.0 6.0) (pow (* -2.0 (log u1)) 0.5))))
1.3s
(log1p (* (/ 1.0 6.0) (pow (* -2.0 (log u1)) 0.5)))

regimes8.5s (0.6%)

Accuracy

Total 22.0b remaining (17.5%)

Threshold costs 7.2b (5.7%)

6.8b39%_divideComplex, real part
3.1b77.1%Octave 3.8, jcobi/1
2.6b34%Octave 3.8, jcobi/4
1.2b55.9%Octave 3.8, jcobi/3
1.2b0%Jmat.Real.lambertw, newton loop step

localize6.2s (0.4%)

bsearch1.9s (0.1%)

end1.0ms (0.0%)