Details

simplify38.1m

Calls

14865 calls:

Slowest
2.0s
(* (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.0s
(* (* (* (* (sqrt (/ 1 PI)) (sqrt (/ 1 PI))) (sqrt (/ 1 PI))) (* (* (exp (* (fabs x) (fabs x))) (exp (* (fabs x) (fabs x)))) (exp (* (fabs x) (fabs x))))) (* (* (+ (/ 15/8 (pow (fabs x) 7)) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x)))) (+ (/ 15/8 (pow (fabs x) 7)) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))) (+ (/ 15/8 (pow (fabs x) 7)) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))))
1.8s
(* (* (* (* (sqrt (/ 1 PI)) (exp (* (fabs x) (fabs x)))) (* (sqrt (/ 1 PI)) (exp (* (fabs x) (fabs x))))) (* (sqrt (/ 1 PI)) (exp (* (fabs x) (fabs x))))) (* (* (+ (/ 15/8 (pow (fabs x) 7)) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x)))) (+ (/ 15/8 (pow (fabs x) 7)) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))) (+ (/ 15/8 (pow (fabs x) 7)) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))))
1.7s
(* (* (* (* (sqrt (/ 1 PI)) (sqrt (/ 1 PI))) (sqrt (/ 1 PI))) (* (* (exp (* (fabs x) (fabs x))) (exp (* (fabs x) (fabs x)))) (exp (* (fabs x) (fabs x))))) (* (* (+ (* 15/8 (/ 1 (pow (fabs x) 7))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x)))) (+ (* 15/8 (/ 1 (pow (fabs x) 7))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))) (+ (* 15/8 (/ 1 (pow (fabs x) 7))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))))
1.6s
(* (* (* (* (sqrt (/ 1 PI)) (exp (* (fabs x) (fabs x)))) (* (sqrt (/ 1 PI)) (exp (* (fabs x) (fabs x))))) (* (sqrt (/ 1 PI)) (exp (* (fabs x) (fabs x))))) (* (* (+ (* 15/8 (/ 1 (pow (fabs x) 7))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x)))) (+ (* 15/8 (/ 1 (pow (fabs x) 7))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))) (+ (* 15/8 (/ 1 (pow (fabs x) 7))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))))

prune3.9m

sample3.7m

Algorithm
86×intervals
Results
1.1m22969×body10240exit
1.1m378440×body80valid
24.4s22569×body1280valid
12.8s18278×body640valid
11.3s70229×body80nan
7.3s1534×body1280nan
5.8s9699×body320valid
3.8s1394×body640nan
1.9s5058×body160valid
1.8s683×body320nan
1.7s41536×pre80true
1.3s1394×body2560valid
883.0ms586×body5120valid
846.0ms316×body160nan

rewrite3.1m

Algorithm
149×rewrite-expression-head
Calls

455 calls:

Slowest
26.1s
(* (+ (+ (/ 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))))
21.0s
(* (* (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)))))
15.7s
(* (/ (+ (+ (+ (+ (+ 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)
15.6s
(* (/ (+ (+ (+ (+ (+ 1 (* 0.1049934947 (* x x))) (* (* (* x x) (* x x)) 0.0424060604)) (* 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)
9.0s
(/ (+ (+ (+ (+ (+ 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)))))
Rules
24572×frac-times
16220×associate-*r/
16074×*-un-lft-identity
11690×times-frac
9840×associate-*l/
9338×frac-add
8091×flip-+ flip3-+
7827×exp-neg
7095×add-sqr-sqrt
6745×exp-diff
6741×neg-sub0
6293×distribute-lft-out
5039×add-cube-cbrt
4079×sqrt-prod
2608×pow1
2446×add-exp-log
2042×add-cbrt-cube
1184×associate-*r*
1099×un-div-inv
1037×associate-*l*
1021×prod-exp
733×cbrt-unprod
684×associate-/r/
678×add-log-exp
632×pow-prod-down
562×pow-prod-up
548×difference-of-squares
529×associate-/l*
455×insert-posit16
439×div-exp
438×swap-sqr
394×div-inv
371×associate-/r*
364×distribute-lft-out--
345×cbrt-undiv
340×sqrt-div
294×cbrt-prod
289×pow-plus
282×unswap-sqr
265×flip3-- flip--
251×pow-sub
227×cbrt-div
221×log-pow
174×pow-sqr
163×*-commutative
153×unpow-prod-down
146×sqrt-pow1
122×exp-prod
117×associate-/l/
112×pow2
111×pow1/2
108×sub-neg
105×log-prod
104×frac-sub
101×pow1/3
98×distribute-rgt-in distribute-lft-in
87×frac-2neg clear-num
82×pow-unpow
77×diff-log
69×pow-exp
68×pow-pow
63×sqr-pow
54×rem-sqrt-square
51×1-exp rec-exp
44×sum-log
42×pow-to-exp
39×associate--l+
35×+-commutative
30×distribute-rgt-out--
29×div-sub log-div
24×cube-unmult
21×cube-prod
19×associate--r+ rem-log-exp cube-div
15×associate--l-
14×pow-flip difference-of-sqr-1
12×unpow-prod-up
11×inv-pow sum-cubes
10×rem-cube-cbrt pow3
pow-div distribute-rgt-out
sqrt-undiv
distribute-rgt1-in
sin-sum sqrt-unprod cube-mult rem-exp-log rem-square-sqrt
associate-+r+ exp-sum unpow3 sub-div exp-to-pow
associate-+l+ remove-posit16 rem-cbrt-cube associate--r- e-exp-1
difference-cubes associate-+l- cos-sum distribute-neg-in

series1.4m

Calls

455 calls:

Slowest
4.5s
(- (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.6s
(+ (* (* (* (/ 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.1s
(* (/ (/ (* 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))) (/ (* (/ (+ alpha i) (sqrt (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))) (/ (+ i beta) (+ (+ alpha beta) (* 2 i)))) (sqrt (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))))
1.2s
(- 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))))))

regimes10.7s

Accuracy

Total 81.8b remaining (41.4%)

12.4b55.7%math.sqrt on complex, real part
11.3b48%math.log/2 on complex, real part
10.4b55.4%math.log10 on complex, real part
10.3b55.9%math.log/1 on complex, real part
9.4b56.7%math.abs on complex

localize7.6s

bsearch4.9s

end2.0ms