Details

rewrite11.7m (74.0%)

Algorithm
151×rewrite-expression-head
Calls

441 calls:

15.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)
15.1s
(* (/ (+ (+ (+ (+ (+ 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)
8.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)))))
8.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)))))
8.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)))))
Rules
7416×times-frac
7408×*-un-lft-identity
6045×add-sqr-sqrt
4958×add-cube-cbrt
2734×add-exp-log
2254×add-cbrt-cube
1885×pow1
1302×prod-diff
1193×prod-exp
1176×sqrt-prod
1003×cbrt-unprod
985×distribute-lft-out
848×pow-prod-down
747×associate-*l*
627×associate-/r*
618×add-log-exp
547×div-exp
527×associate-/l*
480×associate-*r*
467×fma-def
439×log1p-expm1-u expm1-log1p-u
406×associate-/r/
370×cbrt-prod
358×unpow-prod-down
351×cbrt-undiv
333×div-inv
316×difference-of-squares
300×distribute-lft-out--
283×exp-sum
273×flip-+ flip3-+
246×log-pow
241×log-prod
225×distribute-rgt-in distribute-lft-in
197×fma-neg
182×flip3-- flip--
176×associate-*r/
166×unswap-sqr
147×cbrt-div
137×associate--l+
125×associate-*l/
123×cube-prod
115×sqrt-div
114×pow-exp
110×frac-times
109×sqr-pow
104×*-commutative
95×associate-/l/ 1-exp
94×sqrt-pow1
92×pow1/3
84×frac-2neg clear-num
82×sub-neg
73×pow-unpow
72×pow-to-exp
70×diff-log
68×swap-sqr
62×pow1/2
59×exp-prod
55×rec-exp
53×fma-udef
50×frac-add
36×rem-sqrt-square
28×pow-pow pow-prod-up
27×sum-log
22×pow-sqr
21×distribute-lft-neg-in neg-mul-1 distribute-rgt-neg-in
20×+-commutative
19×unpow3 cube-mult
18×pow-plus
16×hypot-def
14×rem-cube-cbrt pow2 associate--r+ associate--l-
12×div-sub
11×rem-log-exp
10×rem-exp-log inv-pow pow-flip
exp-diff
log-div
hypot-udef
cube-div frac-sub
sin-sum associate-+r+ unpow2
sqrt-undiv un-div-inv log1p-udef distribute-lft-neg-out log1p-expm1 expm1-log1p expm1-udef
neg-sub0 neg-log distribute-rgt-neg-out
sqrt-pow2 associate-+r- sqrt-unprod pow3 associate-+l+ associate-+l- unsub-neg rem-square-sqrt cos-sum

sample1.8m (11.1%)

Algorithm
99×intervals
Results
43.3s11937×body10240exit
20.9s483540×body80valid
11.0s26791×body1280valid
4.8s20572×body640valid
2.4s69231×body80nan
1.7s1631×body1280nan
1.6s10563×body320valid
982.0ms9421×body160valid
944.0ms1335×body640nan
609.0ms49536×pre80true
357.0ms1546×body2560valid
337.0ms674×body320nan
134.0ms594×body5120valid
118.0ms335×body160nan

series1.1m (6.8%)

Calls

438 calls:

2.1s
(/ (* (* (/ i (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i))))) (/ (fma i 2.0 (+ alpha beta)) (/ (sqrt (fma beta alpha (* i (+ (+ alpha beta) i)))) (fma i 2.0 (+ alpha beta)))))
1.9s
(/ (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta))))
1.8s
(/ (* (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (fma i 2.0 (+ alpha beta))) (fma beta alpha (* i (+ (+ alpha beta) i))))
1.8s
(* (* (/ i (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i)))))
1.1s
(log (sqrt (exp (- 1.0 (* (* (/ 1 (exp (pow (fabs x) 2))) (/ 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)))))))))))))

prune54.7s (5.8%)

simplify13.9s (1.5%)

Algorithm
249×egg-herbie

regimes5.1s (0.5%)

Accuracy

Total 46.5b remaining (35%)

Threshold costs 6.1b (4.6%)

9.6b21.7%_divideComplex, imaginary part
7.5b33.5%_divideComplex, real part
6.7b0%Octave 3.8, jcobi/2
4.8b16.9%math.sqrt on complex, real part
4.0b33.1%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.4s (0.3%)

bsearch797.0ms (0.1%)

end1.0ms (0.0%)