Details

rewrite13.9m (47.9%)

Algorithm
161×rewrite-expression-head
Calls

492 calls:

1.1m
(/ (+ (/ 9.984369578019572e-06 (+ 7.0 (- z))) (+ (+ (+ (/ 771.3234287776531 (+ 3.0 (- z))) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))) (+ (/ -1259.1392167224028 (+ (- z) 2.0)) (/ -176.6150291621406 (+ (- z) 4.0)))) (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ 12.507343278686905 (+ 5.0 (- z)))) (/ 1.5056327351493116e-07 (+ 8.0 (- z)))))) (exp (+ (+ 0.5 7.0) (- z))))
39.7s
(* (* 2.0 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))
25.2s
(* (/ 1 (sqrt (+ (+ (+ (+ (+ (+ 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))))) (* (/ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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))
19.2s
(* (/ (+ (+ (+ (+ (+ 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)
18.6s
(* (/ (+ (+ (+ (+ (+ 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
17745×times-frac
13297×*-un-lft-identity
11234×add-sqr-sqrt
9431×add-cube-cbrt
5385×add-exp-log
5084×pow1
4328×add-cbrt-cube
4034×sqrt-prod
3121×prod-exp
2492×cbrt-unprod
2414×pow-prod-down
1916×distribute-lft-out
1812×associate-*r*
1391×associate-*l*
1026×prod-diff
883×pow-prod-up
834×div-exp
784×add-log-exp
697×div-inv
675×associate-/l*
617×associate-/r*
594×frac-add
524×cbrt-undiv
514×flip-+ flip3-+
492×log1p-expm1-u
491×expm1-log1p-u
453×sqrt-pow1
377×associate-/r/ unswap-sqr
366×difference-of-squares
351×associate-*r/
329×distribute-lft-out--
323×associate-*l/
320×distribute-rgt-in distribute-lft-in
318×cbrt-prod
315×pow-plus
305×frac-times
298×swap-sqr
295×pow-sqr
292×log-pow
225×pow1/2
220×pow1/3
196×log-prod
195×exp-sum
179×*-commutative
177×flip3-- flip--
172×associate-/l/
157×fma-def
149×fma-neg
145×1-exp
143×rec-exp
140×pow2
139×sum-log
115×sqrt-div
109×cbrt-div
98×frac-2neg clear-num
70×un-div-inv
68×exp-prod
67×sub-neg
62×associate--l+
61×fma-udef
59×unpow-prod-down
55×distribute-rgt-out
52×sqr-pow
51×diff-log
50×associate-+l+
43×cube-unmult
36×+-commutative
30×rem-sqrt-square
23×pow-exp div-sub
22×pow-flip distribute-rgt1-in frac-sub distribute-rgt-neg-in
19×distribute-lft-neg-in
16×sqrt-unprod
14×exp-diff pow-to-exp hypot-udef
11×associate--r+ rem-log-exp
associate-+l- hypot-def
distribute-lft1-in
pow3 inv-pow log-div
pow-unpow cos-sum
associate--l- associate-+r+ neg-mul-1
expm1-log1p expm1-udef sin-sum pow-pow rem-exp-log count-2 pow-sub
neg-sub0 rem-square-sqrt exp-neg
distribute-rgt-out-- associate--r-
unpow-prod-up associate-+r- distribute-rgt-neg-out log1p-expm1 exp-to-pow distribute-neg-in rem-cbrt-cube unsub-neg log1p-udef

simplify6.4m (21.9%)

sample3.5m (12.1%)

Algorithm
99×intervals
Results
1.1m12220×body10240exit
52.3s459403×body80valid
25.5s27609×body1280valid
15.0s21307×body640valid
7.6s68964×body80nan
4.7s10566×body320valid
4.5s1645×body1280nan
2.8s9559×body160valid
2.1s49536×pre80true
2.0s1398×body640nan
836.0ms1531×body2560valid
601.0ms602×body320nan
482.0ms337×body160nan
378.0ms597×body5120valid

prune2.9m (9.8%)

series2.0m (6.8%)

Calls

490 calls:

11.5s
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (/ (- (* 0.254829592 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))))))) (* (/ 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))))))))) (log (/ (exp 0.254829592) (exp (* (/ 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))))))
7.6s
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (/ (- (* 0.254829592 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))))))) (* (/ 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))))))))) (- 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))))))
5.4s
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (/ (- (* 0.254829592 0.254829592) (/ (* (* (* (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma 1.061405429 (/ 1.0 (fma (fabs x) 0.3275911 1.0)) -1.453152027) 1.421413741) -0.284496736) 1.0) (/ 1.0 (fma (fabs x) 0.3275911 1.0))) (fma (* (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma 1.061405429 (/ 1.0 (fma (fabs x) 0.3275911 1.0)) -1.453152027) 1.421413741) (/ 1.0 (fma (fabs x) 0.3275911 1.0))) (* (* (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma 1.061405429 (/ 1.0 (fma (fabs x) 0.3275911 1.0)) -1.453152027) 1.421413741) (/ 1.0 (fma (fabs x) 0.3275911 1.0))) (* (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma 1.061405429 (/ 1.0 (fma (fabs x) 0.3275911 1.0)) -1.453152027) 1.421413741) (/ 1.0 (fma (fabs x) 0.3275911 1.0)))) (* (* -0.284496736 -0.284496736) -0.284496736))) (* (fma (* (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma 1.061405429 (/ 1.0 (fma (fabs x) 0.3275911 1.0)) -1.453152027) 1.421413741) (/ 1.0 (fma (fabs x) 0.3275911 1.0))) (- (* (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma 1.061405429 (/ 1.0 (fma (fabs x) 0.3275911 1.0)) -1.453152027) 1.421413741) (/ 1.0 (fma (fabs x) 0.3275911 1.0))) -0.284496736) (* -0.284496736 -0.284496736)) (fma (fabs x) 0.3275911 1.0)))) (log (/ (exp 0.254829592) (exp (* (/ 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.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.5s
(* (/ (/ (* 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))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))))

regimes13.3s (0.8%)

Accuracy

Total 29.9b remaining (21.1%)

Threshold costs 5.1b (3.6%)

7.4b34.8%_divideComplex, real part
5.8b14.5%Octave 3.8, jcobi/2
3.2b75.5%Octave 3.8, jcobi/1
2.8b0%_divideComplex, imaginary part
2.6b29.4%Octave 3.8, jcobi/4

localize10.2s (0.6%)

bsearch883.0ms (0.1%)

end1.0ms (0.0%)