Details

rewrite10.9m (30.7%)

Algorithm
154×rewrite-expression-head
Calls

451 calls:

22.6s
(* (/ 1.0 (sqrt PI)) (+ (+ (+ (/ (sqrt 1.875) (/ (pow (fabs x) 7) (sqrt 1.875))) (/ (/ 0.5 (fabs x)) (* (fabs x) (fabs x)))) (/ 0.75 (pow (fabs x) 5))) (/ 1.0 (fabs x))))
4.0s
(* (/ 1.0 (sqrt PI)) (+ (+ (+ (/ 1.875 (pow (fabs x) 7)) (/ (/ 0.5 (fabs x)) (* (fabs x) (fabs x)))) (/ 0.75 (pow (fabs x) 5))) (/ 1.0 (fabs x))))
2.9s
(+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (+ (+ (/ 771.3234287776531 (+ 3.0 (- 0 z))) (/ -176.6150291621406 (+ (- 0 z) 4.0))) (/ -1259.1392167224028 (+ (- 0 z) 2.0))))
2.8s
(* (* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1)) d1)
2.7s
(+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (+ (+ (/ 771.3234287776531 (+ 3.0 (- 0 z))) (/ -176.6150291621406 (+ (- 0 z) 4.0))) (/ -1259.1392167224028 (+ (- 0 z) 2.0))))
Rules
47239×frac-times
31197×associate-*r/
18383×associate-*l/
15184×exp-neg
13388×flip-+ flip3-+
13060×frac-add
13022×exp-diff
13015×neg-sub0
6935×*-un-lft-identity
6662×times-frac
6342×add-sqr-sqrt
4839×pow1
4721×add-cube-cbrt
3023×add-exp-log
2323×add-cbrt-cube
2177×un-div-inv
1751×pow-prod-up
1668×sqrt-prod
1349×prod-diff
1319×prod-exp
1311×associate-*r*
996×associate-*l*
989×pow-prod-down
946×pow-plus
909×cbrt-unprod
850×distribute-lft-out
703×div-exp
643×add-log-exp
596×unpow-prod-down
529×associate-/l*
503×pow-sqr
463×associate-/r*
455×difference-of-squares
451×log1p-expm1-u
449×expm1-log1p-u
421×cube-prod
392×cbrt-undiv
362×div-inv
342×pow2
327×distribute-lft-out--
305×cbrt-prod
294×associate-/r/ unswap-sqr
263×pow-exp
256×cbrt-div
197×distribute-rgt-in distribute-lft-in
186×sqr-pow
178×swap-sqr
163×fma-neg
161×fma-def
157×exp-sum
155×*-commutative
150×flip3-- flip--
148×1-exp
140×rec-exp
132×log-pow
118×pow-to-exp
108×pow1/3
99×unpow3 cube-mult
98×sqrt-div
97×associate-/l/
91×frac-2neg clear-num
90×cube-div
79×sum-log
63×fma-udef
60×sqrt-pow1 sub-neg
55×cube-unmult
53×log-prod
48×pow1/2
47×associate--l+
41×pow3
34×diff-log
30×+-commutative
24×pow-unpow
23×associate-+l+
22×frac-sub distribute-rgt-out
21×pow-pow div-sub
20×pow-flip rem-cube-cbrt
19×unpow-prod-up
18×distribute-rgt1-in
17×rem-sqrt-square
16×exp-prod
13×pow-sub
12×sub0-neg cos-sum
10×inv-pow associate-+l- distribute-rgt-out--
9×hypot-def
8×associate-+r+
6×distribute-lft-neg-in associate-+r- distribute-rgt-neg-in unsub-neg
5×associate--r+ rem-exp-log
4×sin-sum neg-log
3×rem-cbrt-cube log1p-udef
2×expm1-log1p expm1-udef log1p-expm1 sqrt-unprod sub-div difference-cubes hypot-udef rem-square-sqrt
1×neg-mul-1 rem-log-exp distribute-neg-in associate--r-

simplify9.7m (27.3%)

prune5.8m (16.3%)

series5.3m (14.8%)

Calls

445 calls:

4.6s
(- (pow 1.0 3) (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
4.2s
(+ (* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
3.6s
(* (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (sqrt (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))) 1)) (/ (sqrt (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))))
3.1s
(* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))))
2.8s
(+ (* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))

sample3.6m (10.0%)

Algorithm
98×intervals
Results
1.2m11978×body10240exit
53.7s450253×body80valid
24.7s27048×body1280valid
12.9s21157×body640valid
6.8s70443×body80nan
4.3s10846×body320valid
3.8s1557×body1280nan
2.8s1342×body640nan
2.0s49536×pre80true
1.4s1395×body2560valid
1.3s5664×body160valid
1.2s693×body320nan
657.0ms327×body160nan
370.0ms616×body5120valid

regimes11.3s (0.5%)

Accuracy

Total 23.6b remaining (17.4%)

Threshold costs 5.0b (3.7%)

7.3b34%_divideComplex, imaginary part
3.3b-14.4%_divideComplex, real part
3.1b75.1%Octave 3.8, jcobi/1
2.3b36.3%Octave 3.8, jcobi/4
1.5b40.8%Octave 3.8, jcobi/3

localize8.3s (0.4%)

bsearch1.5s (0.1%)

end1.0ms (0.0%)