Details

rewrite10.9m (65.6%)

Algorithm
155×rewrite-expression-head
Calls

462 calls:

9.4s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
2.4s
(* (/ 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.8s
(* (/ 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
119964×times-frac
81839×*-un-lft-identity
45284×add-sqr-sqrt
42800×add-cube-cbrt
39556×sqrt-prod
20778×distribute-lft-out
6134×associate-/r/
3446×div-inv
3169×flip-+ flip3-+
2416×add-exp-log
1818×prod-diff
1802×add-cbrt-cube
1365×pow1
1301×fma-def
1078×associate-/r*
1062×associate-/l*
964×distribute-lft-out--
955×difference-of-squares
816×prod-exp
790×associate-*r*
732×associate-*l*
701×add-log-exp
592×div-exp
543×cbrt-unprod
459×distribute-rgt-in log1p-expm1-u distribute-lft-in expm1-log1p-u
357×cbrt-undiv
349×pow-prod-down
295×cbrt-prod
271×log-pow
245×associate-*r/
243×fma-neg
241×log-prod
224×unpow-prod-down
210×associate-*l/
209×sqrt-div
208×associate--l+ unswap-sqr
198×frac-add
168×flip3-- flip--
155×exp-sum
118×1-exp
116×associate-/l/ *-commutative
113×frac-times
111×distribute-rgt-neg-in distribute-lft-neg-in
108×swap-sqr
102×frac-2neg sub-neg clear-num
92×diff-log
91×exp-prod
89×sqrt-pow1
77×pow1/3
73×pow1/2
71×rec-exp pow-exp associate--r+
67×cube-prod
64×cbrt-div pow-unpow
61×sqr-pow
53×associate-+l+
49×associate-+r+
45×pow-to-exp
44×fma-udef sum-log
38×neg-mul-1 pow-prod-up
30×pow-sqr
29×rem-sqrt-square
26×+-commutative
24×pow-flip pow-pow
21×inv-pow
20×rem-log-exp
18×associate--l- hypot-udef
17×rem-exp-log
16×pow-plus div-sub pow2
13×unpow3 cube-mult
12×hypot-def
10×distribute-rgt-out log-div
frac-sub
log-rec
exp-diff distribute-neg-frac
rem-cube-cbrt distribute-rgt-out--
sqrt-unprod sin-sum rem-square-sqrt neg-sub0 distribute-rgt-neg-out unpow2
un-div-inv distribute-lft-neg-out associate-+r- unsub-neg associate-+l-
cube-div 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

prune2.4m (14.2%)

sample1.8m (11.0%)

Algorithm
99×intervals
Results
44.0s11927×body10240exit
23.0s517040×body80valid
10.8s27104×body1280valid
5.0s21470×body640valid
2.5s69203×body80nan
1.7s10603×body320valid
1.6s1583×body1280nan
1.2s9530×body160valid
954.0ms1343×body640nan
597.0ms49536×pre80true
342.0ms1468×body2560valid
326.0ms654×body320nan
141.0ms593×body5120valid
119.0ms358×body160nan

series53.4s (5.4%)

Calls

459 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))))
843.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))
737.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)))))
727.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))))))))
701.0ms
(* (sqrt (/ 1.0 6.0)) (* (sqrt (/ 1.0 6.0)) (pow (* -2.0 (log u1)) 0.5)))

simplify28.7s (2.9%)

Algorithm
253×egg-herbie

regimes5.4s (0.5%)

Accuracy

Total 34.2b remaining (34.4%)

Threshold costs 6.4b (6.4%)

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.6s (0.3%)

bsearch1.1s (0.1%)

end1.0ms (0.0%)