Details

simplify46.3m

Calls

17341 calls:

Slowest
3.2s
(* (* (* (/ 1 (sqrt PI)) (/ 1 (sqrt PI))) (/ 1 (sqrt PI))) (* (* (+ (+ (+ (* 2 (fabs x)) (* (* (* (fabs x) (fabs x)) (fabs x)) 2/3)) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (+ (+ (+ (* 2 (fabs x)) (* (* (* (fabs x) (fabs x)) (fabs x)) 2/3)) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x))))) (+ (+ (+ (* 2 (fabs x)) (* (* (* (fabs x) (fabs x)) (fabs x)) 2/3)) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x))))))
2.7s
(* (* (* (pow (* -2 (log u1)) 0.5) (pow (* -2 (log u1)) 0.5)) (pow (* -2 (log u1)) 0.5)) (/ (* (* (* (* (cos (* (* PI 2) u2)) (cos (* (* PI 2) u2))) (cos (* (* PI 2) u2))) (* (* (pow (* -2 (log u1)) 0.5) (pow (* -2 (log u1)) 0.5)) (pow (* -2 (log u1)) 0.5))) (* (* (* (cos (* (* PI 2) u2)) (cos (* (* PI 2) u2))) (cos (* (* PI 2) u2))) (* (* (pow (* -2 (log u1)) 0.5) (pow (* -2 (log u1)) 0.5)) (pow (* -2 (log u1)) 0.5)))) (* (* (/ 216 (cos (* (* PI 2) u2))) (/ 216 (cos (* (* PI 2) u2)))) (/ 216 (cos (* (* PI 2) u2))))))
2.6s
(* 1 (* (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)))))) x))
2.4s
(* (* (* (pow (* -2 (log u1)) 0.5) (pow (* -2 (log u1)) 0.5)) (pow (* -2 (log u1)) 0.5)) (/ (* (* (* (* (cos (* (* PI 2) u2)) (cos (* (* PI 2) u2))) (cos (* (* PI 2) u2))) (* (* (pow (* -2 (log u1)) 0.5) (pow (* -2 (log u1)) 0.5)) (pow (* -2 (log u1)) 0.5))) (* (* (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5)) (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5))) (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5)))) (/ (* (* 216 216) 216) (* (* (cos (* (* PI 2) u2)) (cos (* (* PI 2) u2))) (cos (* (* PI 2) u2))))))
2.3s
(* (* (* (pow (* -2 (log u1)) 0.5) (pow (* -2 (log u1)) 0.5)) (pow (* -2 (log u1)) 0.5)) (/ (* (* (* (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5)) (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5))) (* (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5)) (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5)))) (* (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5)) (* (cos (* (* PI 2) u2)) (pow (* -2 (log u1)) 0.5)))) (* (* (/ 216 (cos (* (* PI 2) u2))) (/ 216 (cos (* (* PI 2) u2)))) (/ 216 (cos (* (* PI 2) u2))))))

sample5.8m

Algorithm
87×intervals
Results
3.2m52940×body10240exit
1.1m385967×body80valid
26.7s22572×body1280valid
13.5s18090×body640valid
11.1s69111×body80nan
6.0s1521×body1280nan
4.9s9575×body320valid
4.1s1435×body640nan
2.3s5197×body160valid
2.0s41536×pre80true
1.7s655×body320nan
1.3s1424×body2560valid
551.0ms624×body5120valid
378.0ms328×body160nan

prune4.2m

rewrite3.8m

Algorithm
144×rewrite-expression-head
Calls

451 calls:

Slowest
25.5s
(* (+ (+ (/ 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))))
24.4s
(* (/ 1 (sqrt (+ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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)))))) x))
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)))))
16.1s
(* (/ (+ (+ (+ (+ (+ 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))))) x)
12.8s
(* (/ 1 (sqrt (+ (+ (+ (+ (+ (+ 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)))))) (* (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)))) (sqrt (+ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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)))))))) x))
Rules
26697×frac-times
18166×associate-*r/
13517×*-un-lft-identity
10856×times-frac
10545×associate-*l/
9466×frac-add
8731×flip-+ flip3-+
7842×exp-neg
6972×add-sqr-sqrt
6759×exp-diff
6741×neg-sub0
4945×add-cube-cbrt
4791×add-exp-log
4783×distribute-lft-out
4300×pow1
2645×prod-exp
2602×add-cbrt-cube
2134×cbrt-div
1869×sqrt-prod
1594×associate-*r*
1539×pow-prod-up
1169×cbrt-unprod
1104×un-div-inv
1037×associate-*l*
1019×pow-prod-down
942×div-exp
920×pow-plus
894×flip3-- flip--
667×add-log-exp
643×associate-/l*
614×associate-/r/
593×difference-of-squares
454×associate-/r*
451×insert-posit16
395×pow-exp
382×cbrt-undiv
368×div-inv
348×distribute-lft-out--
310×sqrt-div
291×unswap-sqr
283×pow-sqr
274×cbrt-prod
260×log-pow
245×pow-sub
223×pow2
222×sqrt-pow1
202×pow-to-exp
157×*-commutative
143×pow1/3
141×swap-sqr
131×pow1/2
104×log-prod
100×unpow-prod-down
97×associate-/l/
88×distribute-rgt-in distribute-lft-in
78×frac-2neg clear-num
74×1-exp sub-neg
72×rec-exp
69×rem-sqrt-square
68×cube-unmult pow-div
62×sum-log frac-sub
58×exp-prod
56×diff-log
52×sqr-pow
45×+-commutative
44×div-sub
42×log-div
35×associate--l+
33×cube-prod
29×pow3
25×pow-unpow
20×pow-flip inv-pow cube-div
17×pow-pow
16×exp-sum
15×sub0-neg
11×associate-+r+
10×associate--r+
rem-log-exp
sub-div sum-cubes
difference-of-sqr-1 associate--l- distribute-rgt-out
rem-cube-cbrt distribute-rgt-neg-in distribute-lft-neg-in
sqrt-unprod distribute-rgt1-in
sin-sum cube-mult rem-square-sqrt
associate-+l+ unpow-prod-up associate-+r- sqrt-undiv unpow3 distribute-rgt-out-- associate--r-
neg-mul-1
remove-posit16 difference-cubes exp-to-pow rem-exp-log associate-+l- cos-sum distribute-neg-in

series1.4m

Calls

451 calls:

Slowest
4.0s
(- (pow 1 3) (pow (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
3.5s
(+ (* (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1 (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
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.0s
(* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ 1 (/ (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))))
1.2s
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0))

regimes10.9s

Accuracy

Total 103.1b remaining (38.9%)

20.3b33.1%math.log10 on complex, real part
15.2b45.6%math.sqrt on complex, real part
13.7b49.5%math.sqrt on complex, imaginary part, im greater than 0 branch
9.8b57.9%math.log/1 on complex, real part
9.5b56%math.abs on complex

bsearch9.3s

localize6.5s

end1.0ms