Details

rewrite10.7m (62.6%)

Algorithm
164×rewrite-expression-head
Calls

520 calls:

11.3s
(* (* (* (* (* (* (/ (sqrt 1.0) 1) (/ (* 1.0 (sqrt 1.0)) (pow (fabs x) 2))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x))) (/ 1.0 (fabs x)))
3.1s
(log (exp (+ (/ (/ (+ alpha beta) 1) (/ (+ (+ (+ alpha beta) (* 2.0 i)) 2.0) (/ (- beta alpha) (+ (+ alpha beta) (* 2.0 i))))) 1.0)))
2.1s
(+ (* (/ (* (cbrt (+ alpha beta)) (cbrt (+ alpha beta))) (* (* (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) (* (cbrt (+ (+ alpha beta) (* 2.0 i))) (cbrt (+ (+ alpha beta) (* 2.0 i))))) (sqrt 1))) (/ (/ (cbrt (+ alpha beta)) (sqrt 1)) (/ (sqrt (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) (/ (- beta alpha) (cbrt (+ (+ alpha beta) (* 2.0 i))))))) 1.0)
2.1s
(* (/ 1 (/ (+ (+ (* (* x (* (pow (* x x) 3) (pow x 3))) (+ 0.0008327945 (* (* x x) (* 2.0 0.0001789971)))) (+ (* (* x x) (+ 0.7715471019 (* 0.2909738639 (* x x)))) 1.0)) (* (pow x 6) (+ 0.0694555761 (* (* x x) 0.0140005442)))) (- (* (+ (* (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971))) (+ 1.0 (* 0.1049934947 (* x x)))) (+ (* (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971))) (+ 1.0 (* 0.1049934947 (* x x))))) (* (* (pow x 4) (+ 0.0424060604 (* (* x x) 0.0072644182))) (* (pow x 4) (+ 0.0424060604 (* (* x x) 0.0072644182))))))) (/ x (- (+ (* (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971))) (+ 1.0 (* 0.1049934947 (* x x)))) (* (pow x 4) (+ 0.0424060604 (* (* x x) 0.0072644182))))))
973.0ms
(+ (/ (/ (+ alpha beta) 1) (/ (+ (+ (+ alpha beta) (* 2.0 i)) 2.0) (/ (- beta alpha) (+ (+ alpha beta) (* 2.0 i))))) 1.0)
Rules
70274×times-frac
40448×*-un-lft-identity
28820×add-sqr-sqrt
26805×add-cube-cbrt
10560×add-exp-log
6590×associate-*r/
6266×distribute-lft-out
5917×flip-+ flip3-+
5504×prod-exp
4775×add-cbrt-cube
3868×associate-/r/
3574×div-exp
3103×frac-times
2524×frac-add
2400×div-inv
2200×cbrt-unprod
2066×associate-*l*
2007×pow1
1950×frac-sub
1612×sqrt-prod
1305×associate-*l/
1246×cbrt-undiv
1210×associate-*r*
975×associate-/r*
873×add-log-exp
675×difference-of-squares
599×distribute-lft-out--
548×associate-/l*
524×pow-prod-down
379×sqrt-div
361×unpow-prod-down
355×1-exp
314×pow-exp
312×cbrt-prod
305×log-pow
283×unswap-sqr
271×pow-to-exp
254×swap-sqr
187×cube-prod
183×sqrt-pow1
173×associate-/l/
162×*-commutative
140×pow-prod-up
136×flip3-- flip--
130×cbrt-div
129×exp-prod
124×pow1/2 log-prod
116×sqr-pow
108×diff-log
107×pow-plus
98×sub-neg
92×clear-num frac-2neg
85×sum-log
77×distribute-rgt-in distribute-lft-in
73×exp-diff pow-sqr rec-exp
68×exp-neg
67×cube-div pow-pow
65×pow1/3
63×pow-unpow rem-sqrt-square
60×neg-sub0
51×associate--l+
50×+-commutative
45×pow2
37×distribute-rgt-neg-in
35×unpow3 cube-mult
34×distribute-lft-neg-in
30×inv-pow pow-flip exp-sum
29×rem-log-exp
25×log-div
24×rem-cube-cbrt
17×un-div-inv
16×associate--l-
15×div-sub
14×rem-exp-log associate--r+
13×distribute-neg-in
12×distribute-neg-frac
11×associate-+r+
10×associate-+l+
neg-mul-1
pow3 distribute-rgt-neg-out
neg-log
exp-to-pow unpow2
sqrt-undiv associate-+l-
distribute-rgt-out-- rem-cbrt-cube
sqrt-unprod sin-sum log-rec
difference-cubes associate-+r- unsub-neg rem-square-sqrt associate--r- sum-cubes
cos-sum sqrt-pow2 unpow-prod-up pow-sub distribute-frac-neg

sample2.4m (14.2%)

Algorithm
99×intervals
Results
45.9s453510×body80valid
45.4s7720×body10240exit
15.0s26686×body1280valid
7.6s20994×body640valid
5.7s70823×body80nan
2.7s1509×body1280nan
2.6s10495×body320valid
2.1s9561×body160valid
1.6s1249×body640nan
1.1s49536×pre80true
697.0ms1503×body2560valid
554.0ms563×body320nan
363.0ms4180×body80overflowed
307.0ms609×body5120valid
228.0ms333×body160nan

prune2.4m (13.8%)

series1.0m (6.1%)

Calls

518 calls:

3.5s
(- (pow 1.0 3) (* (sqrt (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (+ (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) -1.453152027) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3)) (sqrt (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (+ (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) -1.453152027) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))))
3.2s
(- (pow 1.0 3) (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (+ (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) -1.453152027) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
1.5s
(* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (+ (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) -1.453152027) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (+ (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (+ (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) -1.453152027) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 1.0))
1.5s
(* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (+ (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) -1.453152027) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (+ (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (+ (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) -1.453152027) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 1.0))
889.0ms
(pow (* (pow (log u1) 1.0) (pow -2.0 1.0)) 0.5)

simplify21.9s (2.1%)

Algorithm
262×egg-herbie

regimes7.6s (0.7%)

Accuracy

Total 84.2b remaining (30.6%)

Threshold costs 3.9b (1.4%)

10.1b58.5%math.log10 on complex, real part
10.1b58.7%math.log/2 on complex, real part
9.9b59.1%math.log/1 on complex, real part
9.7b59%math.abs on complex
8.9b58.5%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.9s (0.3%)

bsearch2.1s (0.2%)

end1.0ms (0.0%)