Details

rewrite17.5m (51.1%)

Algorithm
158×rewrite-expression-head
Calls

477 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.7m
(* (+ (+ (+ (+ (+ (/ 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)))
38.9s
(* (/ (sqrt (+ (* (* (* 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 (* 0.0003579942 (* x x)))) 0.0140005442)))) (sqrt (+ (* (* (* 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)))))) x)
23.7s
(/ (sqrt (+ (* (* (* 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 (* 0.0003579942 (* x x)))) 0.0140005442)))) (sqrt (+ (* (* (* 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))))))
21.7s
(* (/ (+ (* (* (* 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
24755×frac-times
16185×associate-*r/
13149×times-frac
11405×*-un-lft-identity
9706×associate-*l/
8672×add-sqr-sqrt
7847×exp-neg
7645×flip-+ flip3-+
6751×frac-add
6735×exp-diff
6730×neg-sub0
6489×add-cube-cbrt
4017×sqrt-prod
3143×add-exp-log
2682×pow1
2535×associate-*r*
2352×add-cbrt-cube
1994×distribute-lft-out
1419×prod-exp
1354×associate-*l*
1128×un-div-inv
847×cbrt-unprod associate-/r/
762×associate-/l*
740×difference-of-squares
738×add-log-exp
644×pow-prod-up
601×distribute-lft-out--
595×pow-prod-down
586×div-exp
453×associate-/r*
441×cbrt-undiv
439×div-inv
336×sqrt-div
321×pow-plus
320×unswap-sqr
294×cbrt-prod
249×log-pow
206×swap-sqr
177×flip3-- flip--
176×cbrt-div
174×*-commutative
162×unpow-prod-down
161×log-prod sqrt-pow1
144×pow-sqr
139×associate-/l/
136×pow1/2
124×pow-exp
109×pow2
103×sum-log
102×distribute-rgt-in distribute-lft-in
94×distribute-rgt-neg-in
89×frac-2neg clear-num
85×sqr-pow
79×pow1/3
66×cube-prod pow-to-exp
58×sub-neg
57×+-commutative
52×cube-unmult
49×rem-sqrt-square
47×diff-log
36×1-exp
32×exp-sum div-sub
31×rec-exp
25×exp-prod
24×cube-div
23×pow3
22×frac-sub
21×pow-unpow
20×associate--l+
19×log-div
16×associate--r+
13×difference-cubes pow-pow
12×associate-+r+ rem-log-exp
11×associate-+r-
10×unpow-prod-up
rem-cube-cbrt
distribute-rgt-out
associate-+l+ pow-div
sin-sum pow-sub
distribute-rgt1-in
pow-flip unpow3 cube-mult sum-cubes inv-pow unsub-neg
associate-+l- associate--r-
sqrt-unprod rem-square-sqrt rem-cbrt-cube
associate--l- distribute-rgt-out-- rem-exp-log cos-sum count-2

simplify7.0m (20.4%)

prune4.3m (12.6%)

sample2.9m (8.6%)

Algorithm
98×intervals
Results
52.3s12154×body10240exit
46.0s437147×body80valid
25.9s27121×body1280valid
11.8s21329×body640valid
6.8s67962×body80nan
5.3s10658×body320valid
3.0s1648×body1280nan
2.0s9643×body160valid
1.8s1340×body640nan
1.5s49536×pre80true
865.0ms1496×body2560valid
833.0ms625×body5120valid
524.0ms669×body320nan
175.0ms339×body160nan

series2.1m (6.0%)

Calls

474 calls:

21.2s
(- (pow 1.0 3) (/ (pow (* (* 1.0 (- (* 0.254829592 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.061405429))))))) (* (/ 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.061405429)))))))))) (exp (- (* (fabs x) (fabs x))))) 3) (pow (* (+ 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.061405429))))))))) 3)))
4.8s
(- (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.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
4.0s
(+ (* (* (* (/ 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.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.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1.0 (* (* (/ 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.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
3.3s
(* (/ (* 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.2s
(* (/ (/ (* 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))))

regimes16.7s (0.8%)

Accuracy

Total 99.9b remaining (39.4%)

Threshold costs 3.0b (1.2%)

10.4b58.8%math.log/2 on complex, real part
10.3b58.7%math.log10 on complex, real part
10.1b59.3%math.log/1 on complex, real part
10.0b58.8%math.abs on complex
9.4b38.4%Octave 3.8, jcobi/4

localize8.7s (0.4%)

bsearch2.3s (0.1%)

end1.0ms (0.0%)