Details

rewrite21.2m (54.9%)

Algorithm
155×rewrite-expression-head
Calls

459 calls:

6.6m
(* (* (* (+ (+ (/ 1.5056327351493116e-07 (+ 8.0 (- z))) (/ 12.507343278686905 (+ (- z) 5.0))) (+ (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ -176.6150291621406 (+ 4.0 (- z)))) (+ (/ 771.3234287776531 (+ 3.0 (- z))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1.0 z)) (/ -1259.1392167224028 (+ 2.0 (- z))))))) (/ 9.984369578019572e-06 (+ 7.0 (- z))))) (/ PI (sin (* PI z)))) (pow (+ (+ 7.0 (- z)) 0.5) (+ (- z) 0.5))) (/ (sqrt (* 2.0 PI)) (exp (+ (+ 7.0 (- z)) 0.5))))
28.2s
(+ (+ (/ 1.5056327351493116e-07 (+ 8.0 (- z))) (/ 12.507343278686905 (+ (- z) 5.0))) (+ (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ -176.6150291621406 (+ 4.0 (- z)))) (+ (/ 771.3234287776531 (+ 3.0 (- z))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1.0 z)) (/ -1259.1392167224028 (+ 2.0 (- z))))))) (/ 9.984369578019572e-06 (+ 7.0 (- z)))))
21.6s
(* (/ (+ (* (* x x) (+ (* (* (* x x) (* (* x x) (* x x))) (+ 0.0005064034 (* (* x x) 0.0001789971))) 0.1049934947)) (+ (+ 1.0 (* (* (* x x) (* x x)) 0.0424060604)) (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))))) (+ (* (* (* (* x x) (* (* x x) (* x x))) (* (* x x) (* x x))) (+ (* (* (* x x) 0.0001789971) 2.0) 0.0008327945)) (+ (* (* x x) (+ (* (* (* x x) (* x x)) (+ (* (* x x) 0.0140005442) 0.0694555761)) (+ 0.7715471019 (* (* x x) 0.2909738639)))) 1.0))) x)
20.9s
(* (/ (+ (* (* x x) (+ (* (* (* x x) (* (* x x) (* x x))) (+ 0.0005064034 (* (* (* x x) (sqrt 0.0001789971)) (sqrt 0.0001789971)))) 0.1049934947)) (+ (+ 1.0 (* (* (* x x) (* x x)) 0.0424060604)) (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))))) (+ (* (* (* (* x x) (* (* x x) (* x x))) (* (* x x) (* x x))) (+ (* (* (* x x) 0.0001789971) 2.0) 0.0008327945)) (+ (* (* x x) (+ (* (* (* x x) (* x x)) (+ (* (* x x) 0.0140005442) 0.0694555761)) (+ 0.7715471019 (* (* x x) 0.2909738639)))) 1.0))) x)
11.8s
(/ (+ (* (* x x) (+ (* (* (* x x) (* (* x x) (* x x))) (+ 0.0005064034 (* (* x x) 0.0001789971))) 0.1049934947)) (+ (+ 1.0 (* (* (* x x) (* x x)) 0.0424060604)) (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))))) (+ (* (* (* (* x x) (* (* x x) (* x x))) (* (* x x) (* x x))) (+ (* (* (* x x) 0.0001789971) 2.0) 0.0008327945)) (+ (* (* x x) (+ (* (* (* x x) (* x x)) (+ (* (* x x) 0.0140005442) 0.0694555761)) (+ 0.7715471019 (* (* x x) 0.2909738639)))) 1.0)))
Rules
25463×frac-times
24730×times-frac
17120×*-un-lft-identity
16199×associate-*r/
12281×add-sqr-sqrt
10881×add-cube-cbrt
10417×associate-*l/
8264×flip-+ flip3-+
8081×frac-add
7847×exp-neg
6932×neg-sub0
6747×exp-diff
5288×sqrt-prod
3276×add-exp-log
3029×distribute-lft-out
2942×pow1
2310×add-cbrt-cube
1533×prod-exp
1528×associate-/l*
1499×associate-*l*
1219×associate-/r/
1202×associate-*r*
1137×un-div-inv
899×cbrt-unprod
803×add-log-exp
792×pow-prod-up
671×pow-prod-down
658×div-inv
639×div-exp
444×associate-/r*
433×cbrt-undiv
369×difference-of-squares
365×pow-plus
348×cbrt-div
293×sqrt-div
265×unswap-sqr
250×cbrt-prod
241×log-pow
234×distribute-lft-out--
214×flip3-- flip--
213×associate-+l-
204×pow-sub
185×pow-sqr
169×sum-log
163×*-commutative
159×log-prod
151×sqrt-pow1
129×pow2
128×pow1/2
123×pow-exp
110×associate-/l/
95×swap-sqr
84×distribute-rgt-in distribute-lft-in
80×frac-2neg clear-num
77×pow1/3
67×pow-to-exp
58×sqr-pow
57×cube-unmult +-commutative
54×sub-neg
53×1-exp
52×diff-log
49×rem-sqrt-square
46×exp-prod
38×frac-sub
36×rec-exp unpow-prod-down
31×log-div
27×exp-sum
23×pow3
21×div-sub
18×associate--r+ cube-div
16×associate-+r+ associate--l+
12×pow-unpow rem-log-exp
10×associate-+l+
pow-pow
pow-flip cube-prod distribute-rgt-out
pow-div
sin-sum difference-cubes inv-pow
associate--l- associate-+r- sum-cubes distribute-rgt1-in
rem-cube-cbrt sub-div rem-exp-log
distribute-rgt-out-- unsub-neg
distribute-lft-neg-out unpow3 sqrt-unprod neg-mul-1 cube-mult distribute-rgt-neg-in rem-square-sqrt cos-sum count-2 rem-cbrt-cube associate--r- distribute-lft-neg-in

simplify8.5m (22.0%)

prune4.5m (11.6%)

sample2.1m (5.5%)

Algorithm
97×intervals
Results
45.4s429062×body80valid
25.2s27121×body1280valid
13.5s21329×body640valid
7.7s8306×body10240exit
6.2s67697×body80nan
4.0s10658×body320valid
2.8s9643×body160valid
1.7s49536×pre80true
1.3s1496×body2560valid
361.0ms625×body5120valid
170.0ms101×body1280nan
111.0ms86×body640nan
58.0ms28×body160nan
36.0ms39×body320nan

series1.9m (4.9%)

Calls

456 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.5s
(- (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))
3.9s
(+ (* (* (* (/ 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.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.3s
(* (/ (/ i (/ (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) 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))))

regimes15.5s (0.7%)

Accuracy

Total 101.1b remaining (39.9%)

Threshold costs 3.8b (1.5%)

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

localize7.9s (0.3%)

bsearch3.0s (0.1%)

end1.0ms (0.0%)