Details

rewrite10.7m (69.9%)

Algorithm
155×rewrite-expression-head
Calls

461 calls:

9.2s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
2.2s
(* (/ 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.6s
(* (/ 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.5s
(+ (/ (/ 1 (/ (/ (fma i 2.0 (+ alpha beta)) (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0)
1.0s
(/ (- (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
22786×times-frac
17970×*-un-lft-identity
13438×add-sqr-sqrt
10890×add-cube-cbrt
2750×distribute-lft-out
2400×add-exp-log
1833×prod-diff
1791×add-cbrt-cube
1304×fma-def
1285×pow1
1205×associate-*l*
1010×associate-/l*
962×distribute-lft-out--
953×difference-of-squares
902×associate-/r/
847×div-inv
797×associate-*r*
762×prod-exp
723×associate-/r*
694×add-log-exp
674×sqrt-prod
633×div-exp
548×flip-+ flip3-+
509×cbrt-unprod
477×distribute-rgt-in distribute-lft-in
458×log1p-expm1-u expm1-log1p-u
382×cbrt-undiv
307×pow-prod-down
289×cbrt-prod
264×unpow-prod-down
248×log-pow
243×fma-neg
238×log-prod
236×associate-*r/
228×unswap-sqr
214×associate--l+
188×flip3-- flip--
160×frac-add
147×exp-sum
146×associate-*l/
132×sqrt-div
127×1-exp
114×*-commutative
111×distribute-rgt-neg-in distribute-lft-neg-in associate-/l/
107×sub-neg
102×swap-sqr
101×frac-2neg clear-num
99×exp-prod
98×cube-prod
94×frac-times
92×diff-log
89×sqrt-pow1
84×rec-exp
78×pow-exp
75×pow1/3
73×pow1/2
71×associate--r+
70×pow-unpow sqr-pow
64×cbrt-div
54×pow-to-exp
53×associate-+l+
49×associate-+r+
48×fma-udef
41×sum-log
38×neg-mul-1
37×pow-prod-up
30×rem-sqrt-square
29×pow-sqr
27×pow-pow
24×pow-flip
23×inv-pow +-commutative
20×rem-log-exp
19×hypot-udef unpow3 cube-mult
18×rem-exp-log associate--l-
16×div-sub
15×pow-plus pow2
12×hypot-def
10×distribute-rgt-out log-div
frac-sub
rem-cube-cbrt log-rec
distribute-neg-frac
exp-diff unpow2
un-div-inv sqrt-unprod sin-sum rem-square-sqrt neg-sub0 distribute-rgt-neg-out
cube-div distribute-rgt-out-- distribute-lft-neg-out associate-+r- unsub-neg associate-+l-
distribute-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.9m (12.2%)

Algorithm
99×intervals
Results
44.2s11927×body10240exit
24.1s517040×body80valid
11.0s27104×body1280valid
5.3s21470×body640valid
2.5s69203×body80nan
1.7s1583×body1280nan
1.6s10603×body320valid
1.0s1343×body640nan
987.0ms9530×body160valid
649.0ms49536×pre80true
359.0ms654×body320nan
343.0ms1468×body2560valid
141.0ms593×body5120valid
130.0ms358×body160nan

prune1.4m (8.8%)

series58.7s (6.4%)

Calls

458 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))))
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)))))
764.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))))))))
763.0ms
(/ 6.0 (* 1.0 (pow (* -2.0 (log u1)) 0.5)))
727.0ms
(* (sqrt (/ 1.0 6.0)) (* (sqrt (/ 1.0 6.0)) (pow (* -2.0 (log u1)) 0.5)))

simplify16.2s (1.8%)

Algorithm
253×egg-herbie

regimes5.4s (0.6%)

Accuracy

Total 34.3b remaining (34.4%)

Threshold costs 6.5b (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.7%Octave 3.8, jcobi/1

localize2.5s (0.3%)

bsearch1.1s (0.1%)

end1.0ms (0.0%)