Details

rewrite26.1m (64.8%)

Algorithm
157×rewrite-expression-head
Calls

480 calls:

2.3m
(* (+ (+ (+ (+ (+ (/ 676.5203681218851 (- 1.0 z)) (/ -1259.1392167224028 (+ 2.0 (- z)))) (/ 771.3234287776531 (+ 3.0 (- z)))) 0.9999999999998099) (/ 12.507343278686905 (+ (- z) 5.0))) (+ (/ -176.6150291621406 (+ (- z) 4.0)) (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ 1.5056327351493116e-07 (+ 8.0 (- z)))) (/ 9.984369578019572e-06 (+ 7.0 (- z)))))) (sqrt (* 2.0 PI)))
1.8m
(* (+ (+ (+ (+ (+ (/ 676.5203681218851 (- 1.0 z)) (/ -1259.1392167224028 (+ 2.0 (- z)))) (/ 771.3234287776531 (+ 3.0 (- z)))) 0.9999999999998099) (/ 12.507343278686905 (+ (- z) 5.0))) (+ (/ -176.6150291621406 (+ (- z) 4.0)) (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ 1.5056327351493116e-07 (+ 8.0 (- z)))) (/ 9.984369578019572e-06 (+ 7.0 (- z)))))) (sqrt (* 2.0 PI)))
21.8s
(* (* (/ 1 (sqrt (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* (* (* x x) 2.0) 0.0001789971))) 0.0140005442)))))) (/ (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0))) (sqrt (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* (* (* x x) 2.0) 0.0001789971))) 0.0140005442))))))) x)
20.9s
(* (/ 1 (sqrt (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* (* (* x x) 2.0) 0.0001789971))) 0.0140005442)))))) (/ (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0))) (sqrt (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* (* (* x x) 2.0) 0.0001789971))) 0.0140005442)))))))
20.9s
(* (/ (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0))) (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* (* (* x x) 2.0) 0.0001789971))) 0.0140005442))))) x)
Rules
31963×times-frac
19728×*-un-lft-identity
16485×add-sqr-sqrt
12426×add-cube-cbrt
4850×sqrt-prod
3644×add-exp-log
3233×frac-times
3103×add-cbrt-cube
3099×pow1
2939×cbrt-div
2876×distribute-lft-out
2808×flip-+ flip3-+
2392×associate-*l/
2249×associate-/r/
1824×associate-*r/
1619×difference-of-squares
1510×prod-exp
1259×associate-/l*
1224×div-inv
1220×associate-*r*
1196×cbrt-unprod
1192×flip3-- flip--
1179×pow-prod-up
1104×associate-*l*
895×div-exp
876×distribute-lft-out--
849×frac-add
840×associate-/l/
788×cbrt-undiv
679×add-log-exp
654×pow-plus
548×sqrt-div
494×associate-/r*
461×pow-prod-down
395×cbrt-prod
312×unswap-sqr
308×log-pow
303×swap-sqr
202×pow-sqr
187×sqrt-pow1
176×pow2
148×*-commutative
147×unpow-prod-down
146×log-prod
140×pow1/2
139×pow-exp
127×exp-prod
124×pow1/3
85×sum-log
84×frac-2neg clear-num
79×rem-sqrt-square
75×distribute-rgt-in distribute-lft-in
74×pow-to-exp
72×sqr-pow
58×cube-unmult sub-neg
49×+-commutative
48×frac-sub
43×1-exp rec-exp
42×exp-sum
33×log-div
31×associate--l+
29×diff-log
27×pow3
24×distribute-rgt-neg-in
22×associate-+l+
20×div-sub
18×rem-log-exp
14×unpow-prod-up associate--r+
13×exp-diff
12×pow-unpow
11×associate-+r-
pow-pow pow-sub
pow-flip cube-prod associate-+r+
inv-pow un-div-inv distribute-rgt-out
neg-sub0 unsub-neg distribute-rgt1-in
sin-sum sub-div rem-exp-log
sqrt-unprod cos-sum
sum-cubes rem-square-sqrt rem-cbrt-cube associate--r-
distribute-rgt-out-- sqrt-pow2 difference-cubes associate-+l- unpow1/3

simplify6.0m (14.9%)

sample3.2m (8.0%)

Algorithm
98×intervals
Results
55.0s12082×body10240exit
51.4s439271×body80valid
24.1s26848×body1280valid
12.6s20636×body640valid
7.0s70273×body80nan
4.8s10752×body320valid
3.2s1658×body1280nan
3.2s9427×body160valid
1.6s1463×body2560valid
1.6s1329×body640nan
1.2s41536×pre80true
839.0ms667×body320nan
432.0ms615×body5120valid
190.0ms339×body160nan

prune2.9m (7.1%)

series1.6m (4.0%)

Calls

475 calls:

5.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 (* (log (exp (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
3.5s
(/ (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))
3.4s
(/ (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))))
3.1s
(* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))))
1.5s
(/ 1.0 (/ 6.0 (pow (* -2.0 (log u1)) 0.5)))

regimes15.9s (0.7%)

Accuracy

Total 97.5b remaining (34.1%)

Threshold costs 4.2b (1.5%)

10.8b56.6%math.log10 on complex, real part
10.5b57.3%math.log/1 on complex, real part
10.4b56.3%math.abs on complex
10.3b51.9%math.sqrt on complex, real part
10.1b58%math.log/2 on complex, real part

localize9.0s (0.4%)

bsearch3.2s (0.1%)

end1.0ms (0.0%)