Details

simplify41.9m

Calls

16044 calls:

Slowest
2.9s
(cbrt (* (* (+ (+ (+ (+ (+ 1 (* 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.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.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))))))
2.2s
(/ (* (* (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))) (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))) (* (* (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) (* (* (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))))
2.0s
(* (/ (/ (* (* (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592)) (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592)) (* (* (+ 1 (* (fabs x) 0.3275911)) (+ 1 (* (fabs x) 0.3275911))) (+ 1 (* (fabs x) 0.3275911)))) (* (* (exp (* (fabs x) (fabs x))) (exp (* (fabs x) (fabs x)))) (exp (* (fabs x) (fabs x))))) (/ (/ (* (* (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592)) (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592)) (* (* (+ 1 (* (fabs x) 0.3275911)) (+ 1 (* (fabs x) 0.3275911))) (+ 1 (* (fabs x) 0.3275911)))) (* (* (exp (* (fabs x) (fabs x))) (exp (* (fabs x) (fabs x)))) (exp (* (fabs x) (fabs x))))))
2.0s
(* (exp (* (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))) (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))))) (exp (* 1 (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))))))
1.8s
(/ (* (* (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) (* (* (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))))

sample9.2m

Algorithm
87×intervals
Results
6.6m100396×body10240exit
1.2m387441×body80valid
25.0s22297×body1280valid
12.9s17421×body640valid
11.9s70528×body80nan
5.9s1573×body1280nan
3.8s8948×body320valid
3.1s1291×body640nan
2.2s33536×pre80true
2.0s4736×body160valid
1.4s1422×body2560valid
1.3s648×body320nan
563.0ms582×body5120valid
435.0ms330×body160nan

rewrite5.1m

Algorithm
150×rewrite-expression-head
Calls

475 calls:

Slowest
33.1s
(cbrt (* (* (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))) (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))))
26.8s
(* (+ (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7))) (+ (+ (/ -0.13857109526572012 (- z -5)) (/ 12.507343278686905 (+ z 4))) (+ (+ (/ 771.3234287776531 (+ z 2)) (+ (/ 676.5203681218851 z) (+ 0.9999999999998099 (/ -1259.1392167224028 (- z -1))))) (/ -176.6150291621406 (+ z 3))))) (exp (- (+ (- z -6) 0.5))))
19.9s
(* (* (pow (+ (- z -6) 0.5) (- z (- 1 0.5))) (sqrt (* PI 2))) (* (+ (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7))) (+ (+ (/ -0.13857109526572012 (- z -5)) (/ 12.507343278686905 (+ z 4))) (+ (+ (/ 771.3234287776531 (+ z 2)) (+ (/ 676.5203681218851 z) (+ 0.9999999999998099 (/ -1259.1392167224028 (- z -1))))) (/ -176.6150291621406 (+ z 3))))) (exp (- (+ (- z -6) 0.5)))))
19.6s
(* (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))) (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))))
19.1s
(* (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))) (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))))
Rules
22195×*-un-lft-identity
18752×times-frac
9811×add-sqr-sqrt
7370×add-cube-cbrt
7057×distribute-lft-out
3198×frac-add
2435×add-exp-log
2261×pow1
1947×add-cbrt-cube
1513×flip-+ flip3-+
1459×associate-*l*
1439×associate-*r*
1332×swap-sqr
1118×frac-times
1094×sqrt-prod
897×prod-exp
884×associate-/r/
816×difference-of-squares
814×distribute-lft-out--
792×associate-*l/
766×associate-/l*
714×add-log-exp
663×associate-*r/
587×sqrt-div
560×cbrt-unprod
527×div-inv
514×div-exp
498×pow-prod-up
494×pow-prod-down
475×insert-posit16
427×cbrt-undiv
425×associate-/r*
394×unpow-prod-down
354×cbrt-prod
336×flip3-- flip--
332×cube-prod
329×unswap-sqr
300×sqrt-pow1
250×exp-diff
246×exp-neg
245×pow-sub
243×neg-sub0
241×pow-plus
157×*-commutative
156×log-pow
155×cbrt-div
153×pow-exp
146×pow-sqr
137×pow1/3
133×pow1/2
97×sqr-pow associate-/l/
95×pow2
87×frac-2neg clear-num
82×pow-to-exp sub-neg
78×log-prod
69×diff-log
66×distribute-rgt-in distribute-lft-in
63×sum-log
54×rem-sqrt-square
46×exp-prod
44×+-commutative
35×associate--l+
33×1-exp rec-exp
30×frac-sub exp-sum
24×log-div
23×pow-flip cube-unmult
21×pow-unpow
17×div-sub
13×rem-cube-cbrt
11×associate--r+ inv-pow un-div-inv
pow3 pow-pow
pow-div associate--l- rem-log-exp
difference-of-sqr-1 sqrt-unprod
cube-mult distribute-rgt-out
unpow3 difference-cubes rem-exp-log cube-div distribute-rgt1-in
associate-+l+ sin-sum rem-square-sqrt
unpow-prod-up distribute-rgt-neg-in associate-+l- rem-cbrt-cube associate--r- distribute-lft-neg-in
remove-posit16 sqrt-undiv associate-+r+ neg-mul-1 sum-cubes cos-sum distribute-neg-in distribute-rgt-out--

prune3.8m

series1.8m

Calls

475 calls:

Slowest
8.3s
(+ (pow (/ (/ (+ alpha beta) (/ (+ (+ alpha beta) (* 2 i)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2 i)) 2.0)) 3) (pow 1.0 3))
4.5s
(- (pow 1 3) (pow (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))) 3))
3.9s
(- (pow 1 3) (* (sqrt (pow (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))) 3)) (sqrt (pow (/ (/ (+ (/ (+ -0.284496736 (/ (+ (/ (+ (/ (/ 1.061405429 (sqrt (+ 1 (* (fabs x) 0.3275911)))) (sqrt (+ 1 (* (fabs x) 0.3275911)))) -1.453152027) (+ 1 (* (fabs x) 0.3275911))) 1.421413741) (+ 1 (* (fabs x) 0.3275911)))) (+ 1 (* (fabs x) 0.3275911))) 0.254829592) (+ 1 (* (fabs x) 0.3275911))) (exp (* (fabs x) (fabs x)))) 3))))
3.4s
(* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))))
3.2s
(/ (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))

regimes10.8s

Accuracy

Total 83.4b remaining (37.2%)

15.6b44.5%math.sqrt on complex, imaginary part, im greater than 0 branch
11.5b58%math.sqrt on complex, real part
10.1b54%math.log/2 on complex, real part
9.7b59.2%math.log10 on complex, real part
9.4b60.1%math.log/1 on complex, real part

bsearch10.0s

localize6.9s

end2.0ms