Details

rewrite10.2m (67.8%)

Algorithm
151×rewrite-expression-head
Calls

440 calls:

494.0ms
(/ (+ (+ (* (* x (* (pow (* x x) 3) (pow x 3))) (+ 0.0008327945 (* (* x x) (* 2.0 0.0001789971)))) (fma (* 0.2909738639 x) (pow x 3) (fma (* 0.7715471019 x) x 1.0))) (* (pow x 6) (+ 0.0694555761 (* (* x x) 0.0140005442)))) (+ (+ (* (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971))) (fma (* 0.1049934947 x) x 1.0)) (* (pow x 4) (+ 0.0424060604 (* (* x x) 0.0072644182)))))
422.0ms
(- (+ (/ (/ x (+ wj 1)) (exp wj)) wj) (/ wj (+ wj 1)))
355.0ms
(+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0)
328.0ms
(* (sqrt (* PI 2.0)) (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)))
288.0ms
(- (/ beta (+ (+ alpha beta) 2.0)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0))
Rules
9311×times-frac
8458×*-un-lft-identity
6660×add-sqr-sqrt
5606×add-cube-cbrt
2840×add-exp-log
1725×add-cbrt-cube
1357×pow1 distribute-lft-out
1255×prod-diff
1033×prod-exp
835×div-exp
657×associate-*r*
653×associate-/l*
626×add-log-exp
616×associate-*l*
568×sqrt-prod
521×div-inv
514×cbrt-unprod
510×associate-/r/
449×associate-/r*
439×log1p-expm1-u
437×expm1-log1p-u
403×cbrt-prod
399×flip-+ flip3-+
388×cbrt-undiv
317×pow-prod-down
265×unpow-prod-down
245×difference-of-squares
242×associate-*r/
236×1-exp
229×distribute-lft-out--
214×unswap-sqr
198×log-pow
194×fma-neg
171×exp-sum
162×distribute-rgt-in distribute-lft-in
160×log-prod
154×fma-def
136×frac-add
134×rec-exp
133×sqrt-pow1
128×associate--l+
117×*-commutative
115×sqrt-div
110×cube-prod
108×flip3-- flip--
95×sqr-pow swap-sqr
89×clear-num pow1/3 frac-2neg
88×pow-exp
86×associate-/l/
77×associate-*l/
73×pow-unpow sub-neg
69×pow1/2 pow-prod-up associate--r+
68×pow-to-exp diff-log
66×cbrt-div
57×fma-udef
55×frac-times
46×exp-prod
45×pow-sqr
41×pow-plus
37×pow-pow sum-log
29×rem-sqrt-square
26×pow2
25×distribute-rgt-neg-in distribute-lft-neg-in
24×+-commutative
23×unpow3 cube-mult
22×log-div
21×div-sub
18×hypot-def
17×neg-mul-1
15×pow-flip
13×rem-cube-cbrt
12×frac-sub associate--l-
11×inv-pow rem-log-exp
10×exp-diff
rem-exp-log
unpow2
associate-+l+
cube-div pow3 sin-sum expm1-log1p unpow1/2 distribute-rgt-neg-out expm1-udef
un-div-inv sqrt-unprod rem-cbrt-cube associate-+l-
sqrt-pow2 exp-neg associate-+r+ rem-square-sqrt mul-1-neg
distribute-lft-neg-out log1p-udef log1p-expm1 hypot-udef unsub-neg distribute-frac-neg neg-sub0 log-rec cos-sum associate--r- exp-to-pow

sample2.6m (17.3%)

Algorithm
99×intervals
Results
54.1s500333×body80valid
46.9s7909×body10240exit
15.5s27438×body1280valid
7.5s20626×body640valid
5.5s69552×body80nan
2.6s10745×body320valid
2.6s1512×body1280nan
2.0s9367×body160valid
1.6s1272×body640nan
1.1s49536×pre80true
696.0ms1516×body2560valid
608.0ms643×body320nan
367.0ms4307×body80overflowed
293.0ms582×body5120valid
212.0ms316×body160nan

prune59.5s (6.6%)

series53.4s (5.9%)

Calls

437 calls:

1.9s
(/ (* (* (- (* (+ (+ 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))))
773.0ms
(/ (* i (+ (+ alpha beta) i)) (/ (* (* (- (* (+ (+ 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)))))
771.0ms
(/ 6.0 (pow (* -2.0 (log u1)) 0.5))
724.0ms
(pow (fma (+ alpha beta) (/ (/ (- beta alpha) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0) 3)
721.0ms
(/ 1 (/ 6.0 (pow (* -2.0 (log u1)) 0.5)))

simplify13.5s (1.5%)

Algorithm
249×egg-herbie

regimes5.2s (0.6%)

Accuracy

Total 42.6b remaining (39%)

Threshold costs 5.1b (4.7%)

10.5b17%_divideComplex, imaginary part
7.2b-7%Octave 3.8, jcobi/2
7.1b38%_divideComplex, real part
4.1b27.2%math.sqrt on complex, imaginary part, im greater than 0 branch
3.1b54.3%Octave 3.8, jcobi/4

localize2.4s (0.3%)

bsearch976.0ms (0.1%)

end1.0ms (0.0%)