Details

simplify19.6m (41.5%)

rewrite10.9m (23.1%)

Algorithm
161×rewrite-expression-head
Calls

479 calls:

15.4s
(* d1 (* (* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1)) d1))
3.4s
(* (/ 1.0 (sqrt PI)) (+ (+ (+ (* 2.0 (fabs x)) (* (/ 2.0 3.0) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1.0 5.0) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))))
2.9s
(/ (- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0)) 2.0)
2.8s
(* (sqrt (* PI 2.0)) (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)))
2.0s
(* (/ 1.0 5.0) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
Rules
10727×times-frac
9190×*-un-lft-identity
8335×add-sqr-sqrt
6336×add-cube-cbrt
3503×add-exp-log
3421×pow1
2285×add-cbrt-cube
1930×sqrt-prod
1648×prod-exp
1551×associate-*r*
1360×pow-prod-up
1240×distribute-lft-out
1179×associate-*l*
867×prod-diff
789×cbrt-unprod
745×cbrt-prod
674×add-log-exp
662×associate-/l*
649×div-exp
626×pow-plus
520×div-inv
500×pow-prod-down
487×cbrt-undiv
478×log1p-expm1-u expm1-log1p-u
477×difference-of-squares
451×associate-*l/
430×associate-/r*
369×unswap-sqr
352×distribute-lft-out--
344×frac-times
280×cbrt-div
267×flip-+ flip3-+ log-prod
260×associate-*r/
255×distribute-rgt-in flip3-- distribute-lft-in flip--
247×pow-sqr
241×swap-sqr
224×log-pow
214×fma-def
207×pow2
205×pow-exp
193×associate-/r/
176×exp-sum
165×*-commutative
142×sqrt-pow1
127×unpow-prod-down
124×fma-neg
115×pow1/3
106×pow-to-exp
93×clear-num frac-2neg
86×pow1/2
84×sqrt-div
75×associate-/l/
68×sum-log
66×1-exp rec-exp sub-neg
62×fma-udef
52×cube-unmult frac-add
43×diff-log
40×unpow-prod-up
34×+-commutative
33×exp-prod
30×sqr-pow
25×pow-flip
23×distribute-rgt-out associate--r+
22×pow3 frac-sub
20×pow-unpow associate-+l+ distribute-rgt1-in
18×div-sub
17×un-div-inv
16×rem-sqrt-square rem-log-exp
15×distribute-rgt-neg-in distribute-lft-neg-in
14×hypot-udef
12×inv-pow
11×associate-+l-
exp-diff
sqrt-unprod pow-sub associate-+r+
rem-exp-log cos-sum log-div
hypot-def
pow-pow distribute-neg-in associate--l+
pow-div sub0-neg sin-sum neg-log neg-mul-1
count-2 rem-square-sqrt associate--r- distribute-rgt-neg-out rem-cbrt-cube
log1p-udef associate-+r- log1p-expm1 unsub-neg expm1-log1p expm1-udef
e-exp-1 exp-neg neg-sub0 log-rec associate--l-

sample8.6m (18.3%)

Algorithm
99×intervals
Results
2.8m12067×body10240exit
2.7m484080×body80valid
47.8s26965×body1280valid
29.1s21055×body640valid
20.8s69961×body80nan
18.1s1379×body640nan
16.2s1545×body1280nan
8.8s9404×body160valid
7.9s10601×body320valid
4.3s49536×pre80true
3.4s1516×body2560valid
2.2s671×body320nan
1.3s318×body160nan
991.0ms622×body5120valid

prune4.8m (10.1%)

series2.5m (5.4%)

Calls

477 calls:

4.3s
(* (* 0.5 (sin re)) (- (fma (pow im 5) 1/60 (* im (fma 1/3 (* im im) 2)))))
3.4s
(- (* 1.0 1.0) (* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (fma (* (cbrt 0.254829592) (cbrt 0.254829592)) (cbrt 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)))) (fma (* (cbrt 0.254829592) (cbrt 0.254829592)) (cbrt 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.2s
(* (* (cbrt (* 0.5 (cos re))) (cbrt (- (* (* (* im im) im) -1/3) (fma (pow im 5) 1/60 (+ im im))))) (cbrt (* (* 0.5 (cos re)) (- (* (* (* im im) im) -1/3) (fma (pow im 5) 1/60 (+ im im))))))
3.0s
(- (* 1.0 1.0) (* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (fma (* (cbrt 0.254829592) (cbrt 0.254829592)) (cbrt 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)))) (fma (* (cbrt 0.254829592) (cbrt 0.254829592)) (cbrt 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)))))))
2.8s
(- 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))))))

regimes24.3s (0.9%)

Accuracy

Total 41.7b remaining (36.1%)

Threshold costs 6.3b (5.5%)

10.7b27%_divideComplex, real part
9.5b41%Octave 3.8, jcobi/4
7.1b35.8%_divideComplex, imaginary part
5.9b55.5%Octave 3.8, jcobi/1
1.2b52.6%Octave 3.8, jcobi/3

localize18.1s (0.6%)

bsearch3.6s (0.1%)

end11.0ms (0.0%)