Details

rewrite10.8m (70.6%)

Algorithm
155×rewrite-expression-head
Calls

463 calls:

9.0s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
2.3s
(* (/ 1 (sqrt (+ (+ (* (* 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)))))) (/ x (/ (- (sqrt (+ (+ (* (* 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 4)) (fma x (* x 0.0072644182) 0.0424060604)) (- (fma (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971)) (fma (* 0.1049934947 x) x 1.0)))))))
1.7s
(* (/ 1 (sqrt (+ (+ (* (* 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)))))) (/ x (/ (sqrt (+ (+ (* (* 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)))))))
1.3s
(+ (/ (/ 1 (/ (/ (fma i 2.0 (+ alpha beta)) (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0)
1.1s
(/ (- (sqrt (+ (+ (* (* 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 4)) (fma x (* x 0.0072644182) 0.0424060604)) (- (fma (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971)) (fma (* 0.1049934947 x) x 1.0)))))
Rules
22850×times-frac
18031×*-un-lft-identity
13506×add-sqr-sqrt
10947×add-cube-cbrt
2750×distribute-lft-out
2400×add-exp-log
1842×prod-diff
1783×add-cbrt-cube
1301×fma-def
1273×pow1
1225×associate-*l*
1010×associate-/l*
964×distribute-lft-out--
955×difference-of-squares
903×associate-/r/
849×div-inv
817×associate-*r*
780×prod-exp
719×associate-/r*
698×add-log-exp
671×sqrt-prod
610×div-exp
549×flip-+ flip3-+
512×cbrt-unprod
476×distribute-rgt-in distribute-lft-in
460×log1p-expm1-u expm1-log1p-u
373×cbrt-undiv
305×pow-prod-down
304×cbrt-prod
270×unpow-prod-down
249×log-pow
244×log-prod
243×fma-neg
234×unswap-sqr
231×associate-*r/
214×associate--l+
189×flip3-- flip--
164×exp-sum
160×frac-add
154×associate-*l/
132×sqrt-div
118×1-exp
117×distribute-rgt-neg-in distribute-lft-neg-in
115×*-commutative
111×associate-/l/
108×sub-neg
107×exp-prod
99×frac-2neg clear-num
98×cube-prod
94×swap-sqr frac-times
92×diff-log
82×pow-exp sqrt-pow1
80×sqr-pow
78×rec-exp
77×pow1/3
71×associate--r+
67×pow-unpow
66×pow1/2
64×cbrt-div
56×pow-to-exp
53×associate-+l+
49×associate-+r+
47×fma-udef
42×sum-log
40×neg-mul-1
32×pow-prod-up
31×rem-sqrt-square
26×pow-sqr
25×pow-pow
24×+-commutative
23×pow-flip
22×inv-pow
20×rem-exp-log rem-log-exp
19×unpow3 cube-mult
18×associate--l- hypot-udef
16×div-sub
14×pow-plus pow2
12×hypot-def
10×distribute-rgt-out log-div
frac-sub
rem-cube-cbrt log-rec
exp-diff distribute-neg-frac
distribute-rgt-out-- unpow2
un-div-inv sin-sum neg-sub0 distribute-rgt-neg-out
cube-div sqrt-unprod associate-+r- unsub-neg rem-square-sqrt associate-+l-
distribute-neg-out distribute-lft-neg-out sub-div distribute-frac-neg neg-log rem-cbrt-cube exp-to-pow
sqrt-pow2 sum-cubes expm1-log1p cos-sum associate--r- expm1-udef

sample1.8m (11.8%)

Algorithm
99×intervals
Results
43.6s11927×body10240exit
22.6s517040×body80valid
10.7s27104×body1280valid
4.9s21470×body640valid
2.4s69203×body80nan
1.6s1583×body1280nan
1.6s10603×body320valid
990.0ms9530×body160valid
960.0ms1343×body640nan
608.0ms49536×pre80true
333.0ms654×body320nan
332.0ms1468×body2560valid
137.0ms593×body5120valid
128.0ms358×body160nan

prune1.3m (8.6%)

series57.9s (6.3%)

Calls

460 calls:

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))))
868.0ms
(/ (fma (* (cos (* (* 2.0 PI) u2)) (cos (* (* 2.0 PI) u2))) (* (/ (* 1.0 1.0) (* 6.0 6.0)) (pow (* -2.0 (log u1)) (* 2 0.5))) (- (* 0.5 0.5))) (- (* (* (/ 1.0 6.0) (pow (* -2.0 (log u1)) 0.5)) (cos (* (* 2.0 PI) u2))) 0.5))
852.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)))))
769.0ms
(/ (/ 1 (sqrt (+ (+ (* (* 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)))))) (* (cbrt (/ (sqrt (+ (+ (* (* 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)))))) (cbrt (/ (sqrt (+ (+ (* (* 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))))))))
702.0ms
(* (sqrt (/ 1.0 6.0)) (* (sqrt (/ 1.0 6.0)) (pow (* -2.0 (log u1)) 0.5)))

simplify14.6s (1.6%)

Algorithm
253×egg-herbie

regimes5.2s (0.6%)

Accuracy

Total 34.2b remaining (34.3%)

Threshold costs 6.4b (6.5%)

7.3b32.6%_divideComplex, real part
5.7b18.7%Octave 3.8, jcobi/2
4.3b22.2%math.sqrt on complex, real part
4.2b31.9%math.sqrt on complex, imaginary part, im greater than 0 branch
3.2b74.8%Octave 3.8, jcobi/1

localize2.5s (0.3%)

bsearch1.1s (0.1%)

end1.0ms (0.0%)