Details

rewrite10.4m (67.4%)

Algorithm
153×rewrite-expression-head
Calls

462 calls:

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)))))
1.0s
(+ (/ (/ 1 (/ (/ (fma i 2.0 (+ alpha beta)) (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0)
532.0ms
(* (+ (/ 1.0 (/ 2.0 (pow (/ 1.0 (fabs x)) 3))) (* (/ 1.0 (fabs x)) (+ (/ (* (/ 1.0 (pow (fabs x) 6)) 15.0) 8.0) (fma (/ 3.0 4.0) (pow (/ 1.0 (fabs x)) 4) 1)))) (* (/ 1.0 (sqrt PI)) (pow (exp (fabs x)) (/ (fabs x) 2))))
493.0ms
(* (/ 1.0 5.0) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
Rules
20219×times-frac
16117×*-un-lft-identity
12439×add-sqr-sqrt
9590×add-cube-cbrt
2358×distribute-lft-out
2304×add-exp-log
1687×add-cbrt-cube
1390×pow1
1322×prod-diff
1317×fma-def
1293×cbrt-prod
1164×sqrt-prod
1154×associate-*r*
1061×difference-of-squares
1021×associate-*l*
1012×distribute-lft-out--
857×associate-/r/
853×associate-/r*
781×associate-/l*
718×prod-exp
671×add-log-exp
621×div-inv
614×div-exp
461×expm1-log1p-u
460×log1p-expm1-u
446×flip-+ flip3-+
434×cbrt-unprod
364×cbrt-undiv
349×unpow-prod-down
325×pow-prod-down
274×log-pow
263×unswap-sqr
210×fma-neg
203×cube-prod
202×flip3-- flip--
198×exp-prod
184×distribute-rgt-in distribute-lft-in
175×associate-*r/
158×exp-sum
156×log-prod
141×associate--l+
137×1-exp
130×pow1/3
127×sqrt-pow1
126×sqrt-div
122×*-commutative
111×associate-*l/
108×cbrt-div
105×pow-prod-up
102×swap-sqr
95×frac-times
93×pow-exp
89×pow-unpow
85×sqr-pow
83×clear-num frac-2neg
81×sub-neg
77×associate-/l/
75×rec-exp
73×pow1/2
72×diff-log
71×associate--r+
70×pow-pow
68×pow-to-exp
63×pow-sqr
55×distribute-rgt-neg-in
51×distribute-lft-neg-in
47×fma-udef
46×frac-add pow-plus
42×pow-flip
39×sum-log
33×unpow3 cube-mult neg-mul-1
28×rem-sqrt-square +-commutative
26×pow2
23×inv-pow
22×rem-cube-cbrt
18×rem-exp-log
16×exp-diff rem-log-exp
15×hypot-udef log-div
14×hypot-def associate--l-
12×associate-+l+ neg-sub0 div-sub
10×cube-div exp-neg frac-sub
un-div-inv
pow3 associate-+l-
sqrt-unprod sin-sum rem-cbrt-cube unpow2
rem-square-sqrt neg-log
sqrt-undiv sqrt-pow2 pow-sub expm1-log1p distribute-rgt-neg-out expm1-udef
distribute-lft-neg-out log1p-udef log1p-expm1 unpow-prod-up cos-sum associate--r- distribute-rgt-out-- pow-neg exp-to-pow

sample2.4m (15.4%)

Algorithm
99×intervals
Results
46.3s545636×body80valid
46.3s7869×body10240exit
14.3s27095×body1280valid
6.9s20772×body640valid
3.8s69737×body80nan
2.6s1536×body1280nan
2.4s10732×body320valid
1.8s9537×body160valid
1.5s1232×body640nan
936.0ms49536×pre80true
585.0ms667×body320nan
557.0ms1420×body2560valid
275.0ms4291×body80overflowed
254.0ms610×body5120valid
196.0ms343×body160nan

prune1.2m (7.6%)

series1.1m (7.1%)

Calls

460 calls:

2.2s
(pow (exp (* (cbrt (log (fma (/ (cbrt (pow (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) 1.061405429 -1.453152027) 1.421413741) -0.284496736) 0.254829592) 3)) (exp (* (fabs x) (fabs x)))) (/ (- 1.0) (fma (fabs x) 0.3275911 1.0)) 1.0))) (cbrt (log (fma (/ (cbrt (pow (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) 1.061405429 -1.453152027) 1.421413741) -0.284496736) 0.254829592) 3)) (exp (* (fabs x) (fabs x)))) (/ (- 1.0) (fma (fabs x) 0.3275911 1.0)) 1.0))))) (cbrt (log (fma (/ (cbrt (pow (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) (fma (/ 1.0 (fma (fabs x) 0.3275911 1.0)) 1.061405429 -1.453152027) 1.421413741) -0.284496736) 0.254829592) 3)) (exp (* (fabs x) (fabs x)))) (/ (- 1.0) (fma (fabs x) 0.3275911 1.0)) 1.0))))
2.0s
(fma (/ (/ 1 (/ (/ (* (cbrt (fma i 2.0 (+ alpha beta))) (cbrt (fma i 2.0 (+ alpha beta)))) 1) (* (cbrt (- beta alpha)) (cbrt (- beta alpha))))) (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (/ (/ 1 (/ (/ (cbrt (fma i 2.0 (+ alpha beta))) (+ alpha beta)) (cbrt (- beta alpha)))) (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) 1.0)
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))))
1.9s
(/ (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta))))
1.0s
(* (/ (sqrt i) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (/ (fma i 2.0 (+ alpha beta)) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta))))))

simplify14.3s (1.6%)

Algorithm
251×egg-herbie

regimes5.1s (0.5%)

Accuracy

Total 46.2b remaining (43.6%)

Threshold costs 5.6b (5.3%)

8.8b45.5%Octave 3.8, jcobi/4
8.1b32.2%_divideComplex, imaginary part
7.3b34.2%_divideComplex, real part
6.1b5.1%Octave 3.8, jcobi/2
4.6b20.3%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.3s (0.2%)

bsearch1.1s (0.1%)

end1.0ms (0.0%)