Details

rewrite10.4m (69.2%)

Algorithm
151×rewrite-expression-head
Calls

453 calls:

9.0m
(/ (* (* (sqrt (* PI 2.0)) (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -1259.1392167224028 (+ (- z 1.0) 2.0))) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (/ -176.6150291621406 (+ (- z 1.0) 4.0))) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (/ -0.13857109526572012 (+ (- z 1.0) 6.0))) (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0))) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)))) (exp (+ (+ (- z 1.0) 7.0) 0.5)))
6.8s
(* (/ (+ (+ (+ (+ (+ 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)
3.9s
(/ (+ (+ (+ (+ (+ 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)))))
3.9s
(/ (+ (+ (+ (+ (+ 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)))))
591.0ms
(log (exp (- (pow (/ alpha (+ (+ alpha beta) 2.0)) 3) (pow 1.0 3))))
Rules
11150×*-un-lft-identity
10846×times-frac
8614×add-sqr-sqrt
7142×add-cube-cbrt
3415×sqrt-prod
2874×prod-diff
2693×add-exp-log
2028×distribute-lft-out
1835×add-cbrt-cube
1639×unpow-prod-down
1509×cube-prod
1278×associate-*r*
1237×pow1
1010×associate-*l*
940×prod-exp
708×log-prod
707×exp-sum
698×div-exp
657×add-log-exp
616×associate-/l*
540×cbrt-unprod
508×div-inv
472×associate-/r/
451×log1p-expm1-u
450×expm1-log1p-u
448×cbrt-prod
406×cbrt-undiv
379×associate-/r*
329×unswap-sqr
324×flip-+ flip3-+
302×sqr-pow
290×pow-prod-down difference-of-squares
252×fma-neg
237×log-pow
227×distribute-lft-out--
217×unpow3 cube-mult
165×distribute-rgt-in distribute-lft-in
152×fma-def
148×1-exp
137×associate--l+
133×*-commutative
129×associate-*r/
125×associate-*l/
121×pow-exp
113×sqrt-div
110×flip3-- flip--
99×pow-unpow
92×rec-exp
90×clear-num frac-2neg
80×pow-to-exp
79×exp-prod
77×frac-times
76×sqrt-pow1
75×sub-neg
69×associate-/l/ associate--r+
67×diff-log
66×swap-sqr
59×pow1/2 pow1/3
47×fma-udef
46×rem-sqrt-square
44×sum-log
42×frac-add
41×distribute-rgt-neg-in distribute-lft-neg-in
40×cbrt-div
34×pow-pow
31×+-commutative
25×neg-mul-1
22×rem-cube-cbrt
21×div-sub
17×pow-flip
16×pow-prod-up pow-sqr
14×hypot-udef
13×associate-+l+
12×cube-div associate--l- inv-pow hypot-def
11×rem-exp-log pow2 pow-plus
10×rem-log-exp
un-div-inv sqrt-unprod frac-sub associate-+l-
exp-diff sin-sum neg-sub0 distribute-rgt-neg-out log-div
difference-cubes distribute-lft-neg-out unpow2
unsub-neg distribute-frac-neg neg-log expm1-udef exp-to-pow
sqrt-undiv sqrt-pow2 log1p-udef associate-+r+ exp-neg log1p-expm1 unpow-prod-up pow-sub expm1-log1p cos-sum associate--r- rem-cbrt-cube

sample2.2m (14.5%)

Algorithm
99×intervals
Results
46.5s8102×body10240exit
34.8s451460×body80valid
14.6s27112×body1280valid
7.1s21052×body640valid
4.2s69832×body80nan
2.9s1601×body1280nan
2.5s10957×body320valid
1.8s1365×body640nan
1.8s9526×body160valid
942.0ms49536×pre80true
637.0ms662×body320nan
591.0ms1472×body2560valid
262.0ms609×body5120valid
208.0ms341×body160nan
192.0ms3985×body80overflowed

series1.1m (7.2%)

Calls

450 calls:

3.5s
(* (* (/ (cbrt (- beta alpha)) (/ (fabs (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (cbrt (- beta alpha)))) (/ (+ alpha beta) (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)))) (/ (/ (cbrt (- beta alpha)) (sqrt (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)))) (fma i 2.0 (+ alpha beta))))
2.1s
(fma (cbrt (+ (* (/ x (+ wj 1)) (/ 1 (exp wj))) wj)) (* (cbrt (+ (* (/ x (+ wj 1)) (/ 1 (exp wj))) wj)) (cbrt (+ (* (/ x (+ wj 1)) (/ 1 (exp wj))) wj))) (/ (- (pow (cbrt wj) 3)) (+ wj 1)))
2.0s
(/ (* (* (- (* (+ (+ 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.1s
(+ (* (/ (+ alpha beta) (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (/ (* (/ (cbrt (- beta alpha)) (/ (fabs (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (cbrt (- beta alpha)))) (/ (cbrt (- beta alpha)) (sqrt (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))))) (fma i 2.0 (+ alpha beta)))) 1.0)
1.1s
(+ (* (* (/ (cbrt (- beta alpha)) (/ (fabs (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (cbrt (- beta alpha)))) (/ (+ alpha beta) (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)))) (/ (/ (cbrt (- beta alpha)) (sqrt (cbrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)))) (fma i 2.0 (+ alpha beta)))) 1.0)

prune59.8s (6.6%)

simplify13.9s (1.5%)

Algorithm
249×egg-herbie

regimes5.0s (0.6%)

Accuracy

Total 44.5b remaining (30.3%)

Threshold costs 5.1b (3.4%)

8.8b40%_divideComplex, imaginary part
6.8b38.5%_divideComplex, real part
5.2b19.2%Octave 3.8, jcobi/2
5.0b18.3%math.sqrt on complex, imaginary part, im greater than 0 branch
4.3b30.4%math.sqrt on complex, real part

localize2.4s (0.3%)

bsearch1.1s (0.1%)

end1.0ms (0.0%)