Details

simplify5.4m (30.5%)

rewrite4.9m (27.7%)

Algorithm
161×rewrite-expression-head
Calls

482 calls:

3.2m
(* (/ (* (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))) (+ (+ (+ (/ 12.507343278686905 (+ (- z 1.0) 5.0)) (/ -0.13857109526572012 (+ (- z 1.0) 6.0))) (+ (+ (+ (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (/ -176.6150291621406 (+ (- z 1.0) 4.0))) (+ (/ 676.5203681218851 z) 0.9999999999998099))) (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)))))
34.0s
(+ (+ (/ 1.5056327351493116e-07 (+ (- z) 8.0)) (/ 12.507343278686905 (+ 5.0 (- z)))) (+ (/ -1259.1392167224028 (+ 2.0 (- z))) (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (/ -176.6150291621406 (+ (- z) 4.0))) (/ 771.3234287776531 (+ 3.0 (- z)))) (+ (/ 9.984369578019572e-06 (+ (- z) 7.0)) (/ -0.13857109526572012 (+ 6.0 (- z)))))))
14.4s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
13.5s
(+ (/ -1259.1392167224028 (+ 2.0 (- z))) (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (/ -176.6150291621406 (+ (- z) 4.0))) (/ 771.3234287776531 (+ 3.0 (- z)))) (+ (/ 9.984369578019572e-06 (+ (- z) 7.0)) (/ -0.13857109526572012 (+ 6.0 (- z))))))
4.4s
(+ (* (- (/ (+ (/ 1.061405429 (pow (fma (fabs x) 0.3275911 1.0) 4)) (- (- (+ (/ 1.421413741 (pow (fma (fabs x) 0.3275911 1.0) 2)) 0.254829592) (/ 1.453152027 (pow (fma (fabs x) 0.3275911 1.0) 3))) (/ 0.284496736 (fma 0.3275911 (fabs x) 1.0)))) (exp (pow (fabs x) 2)))) (/ 1.0 (fma (fabs x) 0.3275911 1.0))) 1.0)
Rules
9368×times-frac
8956×*-un-lft-identity
8436×add-sqr-sqrt
6421×add-cube-cbrt
2877×add-exp-log
2307×sqrt-prod
2169×add-cbrt-cube
1931×prod-diff
1485×pow1
1427×frac-add
1214×associate-*l*
1108×distribute-lft-out
1073×prod-exp
830×associate-*r*
816×add-log-exp
742×associate-/l*
725×cbrt-unprod
704×div-exp
643×flip-+ flip3-+
532×associate-*r/
529×associate-/r*
494×frac-times
484×expm1-log1p-u
483×log1p-expm1-u
477×difference-of-squares
474×cbrt-undiv
466×cbrt-prod
445×div-inv
432×unpow-prod-down
356×unswap-sqr
342×pow-prod-down
314×distribute-rgt-in distribute-lft-in
305×exp-sum
302×associate-*l/
267×distribute-lft-out--
253×log-pow
229×swap-sqr
205×cube-prod
204×fma-neg
194×associate-/r/
183×sum-log
175×pow-exp
167×sqrt-pow1
151×fma-def
139×*-commutative
137×exp-prod
136×pow1/3
132×distribute-rgt-neg-in pow-prod-up
122×sqr-pow flip3-- flip--
117×distribute-lft-neg-in
108×associate-+r+
103×clear-num frac-2neg
94×pow-to-exp
90×pow1/2 pow-plus
84×associate-/l/
78×fma-udef sub-neg
76×sqrt-div
64×pow-sqr
60×cube-div
57×associate--l+
55×pow-unpow
54×unpow3 cube-mult rec-exp log-prod cbrt-div 1-exp
45×diff-log
34×pow2
30×neg-mul-1
29×+-commutative
26×rem-sqrt-square
18×div-sub
16×frac-sub
15×pow-pow hypot-udef
12×unpow2
11×associate-+l+
10×rem-cube-cbrt un-div-inv pow-flip
associate-+l-
sqrt-unprod pow3 inv-pow sin-sum associate--r+
distribute-rgt-neg-out
exp-diff exp-neg distribute-neg-in neg-sub0 rem-log-exp
rem-square-sqrt cos-sum expm1-udef associate--l-
log1p-udef distribute-lft-neg-out log1p-expm1 expm1-log1p
sqrt-undiv associate-+r- pow-sub unsub-neg neg-log associate--r-
rem-exp-log sqrt-pow2 unpow-prod-up distribute-frac-neg distribute-rgt-out hypot-def log-div

sample3.8m (21.1%)

Algorithm
100×intervals
Results
1.1m15701×body10240exit
52.6s445672×body80valid
33.8s30483×body1280valid
15.8s23434×body640valid
6.7s70918×body80nan
5.1s11584×body320valid
4.9s3075×body1280nan
3.7s2577×body640nan
3.0s10052×body160valid
1.5s49536×pre80true
1.5s1237×body320nan
837.0ms1477×body2560valid
804.0ms630×body160nan
802.0ms634×body5120valid

prune1.9m (10.6%)

series1.5m (8.7%)

Calls

482 calls:

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

regimes8.8s (0.8%)

Accuracy

Total 35.3b remaining (24.8%)

Threshold costs 7.9b (5.6%)

10.9b24.3%_divideComplex, real part
10.3b29.8%_divideComplex, imaginary part
3.7b3.1%Octave 3.8, jcobi/4
1.2b54.5%Octave 3.8, jcobi/3
1.1b0%Jmat.Real.gamma, branch z greater than 0.5

localize5.0s (0.5%)

bsearch1.5s (0.1%)

end1.0ms (0.0%)