Details

simplify9.7m (34.7%)

sample5.1m (18.3%)

Algorithm
100×intervals
Results
1.6m480596×body80valid
1.5m15179×body10240exit
43.2s30002×body1280valid
21.1s22692×body640valid
10.6s70927×body80nan
9.4s2903×body1280nan
6.8s11733×body320valid
5.4s2420×body640nan
4.1s10381×body160valid
2.3s1225×body320nan
2.1s49536×pre80true
1.4s1470×body2560valid
642.0ms640×body160nan
559.0ms586×body5120valid

prune4.7m (16.7%)

rewrite4.3m (15.2%)

Algorithm
160×rewrite-expression-head
Calls

503 calls:

27.7s
(* (* (/ (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))))
27.5s
(* (* (* (/ (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 0.5))))
23.6s
(* (/ (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))))
19.1s
(* (* (/ (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))))
17.3s
(/ (* (* (* (sqrt (* PI 2)) (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)) (* (exp (+ 0.5 z)) (pow (+ 0.5 (- z -6)) (- 1 0.5))))
Rules
25436×frac-times
16759×associate-*r/
11790×frac-add
10377×associate-*l/
8639×flip-+ flip3-+
7882×exp-neg
7720×*-un-lft-identity
6759×exp-diff
6756×neg-sub0
6645×times-frac
6256×add-sqr-sqrt
6185×pow1
4583×add-cube-cbrt
3945×sqrt-pow1
2910×add-exp-log
2770×pow-prod-up
2381×sqrt-prod
2158×add-cbrt-cube
2061×pow1/2
1604×distribute-lft-out
1465×log-pow
1364×prod-exp
1304×pow-prod-down
1236×prod-diff
1152×associate-*l*
1136×un-div-inv
1097×pow-sqr
1085×associate-*r*
964×add-log-exp
788×cbrt-unprod
624×associate-/l*
608×div-exp
511×log1p-expm1-u
509×pow-plus
504×expm1-log1p-u
503×insert-posit16
480×swap-sqr
462×associate-/r*
421×unpow-prod-down
382×associate-/r/
371×unswap-sqr
335×cbrt-undiv
323×cube-prod
305×pow-sub
291×sum-log
290×div-inv
276×pow2
263×cbrt-prod
247×difference-of-squares
232×distribute-lft-out--
226×distribute-rgt-in distribute-lft-in
190×sqr-pow
179×fma-neg
178×cbrt-div
166×fma-def
160×*-commutative
147×associate-/l/
140×exp-sum
135×flip3-- flip--
126×pow-exp
123×cube-div
113×sqrt-div
111×log-prod
103×1-exp
101×rec-exp
92×frac-2neg clear-num
89×cube-mult
88×unpow3
80×pow-to-exp
74×pow1/3 fma-udef
64×pow-div
59×sub-neg
52×pow-flip
48×inv-pow
47×unpow-prod-up
43×pow-unpow
42×diff-log
36×cos-sum
32×rem-sqrt-square
30×cube-unmult distribute-rgt-out
29×distribute-rgt1-in +-commutative
26×rem-cube-cbrt pow-pow
25×sqrt-unprod
22×hypot-udef
21×sin-sum
20×frac-sub
19×associate-+l+ hypot-def
16×associate-+r+
15×pow3 rem-square-sqrt
14×sqrt-undiv
13×div-sub exp-prod
expm1-log1p expm1-udef distribute-lft-neg-in distribute-rgt-neg-in
rem-log-exp
log1p-expm1 associate--r+ distribute-neg-in log1p-udef
associate-+l-
associate-+r- associate--l+
distribute-rgt-neg-out neg-mul-1 associate--r-
sub0-neg difference-cubes rem-cbrt-cube
remove-posit16 associate--l- distribute-rgt-out-- exp-to-pow unsub-neg log-div

series3.9m (14.0%)

Calls

503 calls:

2.3m
(- (pow (pow 1 3) 3) (pow (* (sqrt (pow (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3)) (sqrt (pow (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))) 3))
4.0s
(- (pow 1 3) (pow (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
3.8s
(* (+ (sqrt (+ (- (* d1 d2) (* d1 d3)) (* d4 d1))) d1) (- (sqrt (+ (- (* d1 d2) (* d1 d3)) (* d4 d1))) d1))
3.6s
(+ (* (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1 (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
3.0s
(- (pow 1 3) (* (sqrt (pow (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3)) (sqrt (pow (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))))

regimes9.1s (0.5%)

Accuracy

Total 20.7b remaining (17.5%)

Threshold costs 7.3b (6.2%)

7.1b36%_divideComplex, real part
2.7b62.6%Octave 3.8, jcobi/4
2.4b0%_divideComplex, imaginary part
1.1b54%Octave 3.8, jcobi/3
0.8b0%Jmat.Real.gamma, branch z greater than 0.5

localize7.9s (0.5%)

bsearch1.2s (0.1%)

end1.0ms (0.0%)