Details

rewrite11.7m (44.5%)

Algorithm
166×rewrite-expression-head
Calls

484 calls:

9.6m
(/ (* (* (sqrt (* PI 2.0)) (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -1259.1392167224028 (+ (- z 1.0) 2.0))) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (/ -176.6150291621406 (+ (- z 1.0) 4.0))) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (/ -0.13857109526572012 (+ (- z 1.0) 6.0))) (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0))) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)))) (exp (+ (+ (- z 1.0) 7.0) 0.5)))
18.9s
(* (/ (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))) (/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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)
17.7s
(/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))))
13.9s
(* (/ (- (- (+ (+ (+ (+ 1.0 (* 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 (pow x 6)) (pow x 4))) (- (- (+ (+ (+ (+ (+ 1.0 (* 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 0.0001789971) (pow x 8)) (pow x 4)))) x)
10.9s
(/ (sqrt (+ (+ (+ (+ (+ 1.0 (* 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))))) (/ (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 0.0001789971)) (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2)))))) (* (pow x 2) (+ (* 0.0140005442 (* (* (* x x) x) (* (* x x) x))) (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))) (sqrt (+ (+ (+ (+ (+ 1.0 (* 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)))))))
Rules
27649×times-frac
17700×*-un-lft-identity
14588×add-sqr-sqrt
12166×add-cube-cbrt
2826×add-exp-log
2661×sqrt-prod
2385×distribute-lft-out
2297×add-cbrt-cube
1998×pow1
1630×associate-*r*
1421×associate-/l*
1125×associate-*l*
1100×associate-/r/
1058×div-inv
972×prod-exp
874×flip-+ flip3-+
793×cbrt-unprod
775×difference-of-squares
719×add-log-exp
693×div-exp
544×unpow-prod-down frac-times
507×associate-/r*
499×pow-prod-down
471×cbrt-prod
466×cbrt-undiv
449×distribute-lft-out--
438×cube-prod
422×associate-*r/
351×unswap-sqr
347×sqr-pow
330×swap-sqr
321×cbrt-div
292×associate-*l/
272×sqrt-div
260×flip3-- flip--
212×sqrt-pow1
211×pow-prod-up
207×log-pow
162×*-commutative
138×pow1/2
133×pow-unpow
130×exp-prod
120×pow-sqr
116×log-prod
114×pow-plus
110×distribute-rgt-in distribute-lft-in
103×associate-/l/
102×frac-add
98×pow1/3
88×pow-exp
85×frac-2neg clear-num
80×sub-neg
75×diff-log
65×1-exp
62×rem-sqrt-square
61×pow2 pow-pow
60×pow-to-exp
57×sum-log
50×exp-neg
49×exp-diff
47×rec-exp
46×distribute-rgt-neg-in
44×neg-sub0
42×distribute-lft-neg-in
34×associate--l+
33×+-commutative
32×cube-div
26×unpow3 cube-mult
24×distribute-neg-frac
18×div-sub
17×rem-cube-cbrt
16×pow3
15×associate-+r+
14×associate-+l+ neg-mul-1 log-div
13×unpow2
12×associate--r+
11×un-div-inv associate--l- exp-sum
inv-pow sin-sum pow-flip
frac-sub rem-log-exp
difference-cubes
distribute-neg-in rem-cbrt-cube associate-+l-
sqrt-pow2 rem-exp-log
sqrt-undiv sqrt-unprod rem-square-sqrt distribute-rgt-neg-out
cos-sum distribute-lft-neg-out
unpow1/3 distribute-rgt-out-- unpow-prod-up pow-sub unsub-neg distribute-rgt-out associate--r- exp-to-pow

series10.8m (41.0%)

Calls

477 calls:

19.9s
(- (pow 1.0 3) (/ (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (- (* 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 (- 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)))
3.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 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
1.5s
(* (* (* (/ 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))
829.0ms
(- 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))))))
788.0ms
(pow (* (/ 1.0 6.0) (pow (* -2.0 (log u1)) 0.5)) 3)

sample2.1m (8.2%)

Algorithm
98×intervals
Results
45.2s7687×body10240exit
35.9s458085×body80valid
14.1s27267×body1280valid
6.9s21054×body640valid
3.8s68309×body80nan
2.7s1644×body1280nan
2.4s10373×body320valid
1.6s1308×body640nan
955.0ms49536×pre80true
816.0ms5539×body160valid
576.0ms1427×body2560valid
521.0ms601×body320nan
277.0ms4306×body80overflowed
259.0ms598×body5120valid
199.0ms344×body160nan

prune1.2m (4.5%)

simplify16.3s (1.0%)

Algorithm
262×egg-herbie

regimes6.3s (0.4%)

Accuracy

Total 101.7b remaining (41.3%)

Threshold costs 3.6b (1.5%)

11.6b53.5%math.log10 on complex, real part
11.6b53.7%math.log/1 on complex, real part
11.2b53.6%math.abs on complex
10.3b58.7%math.log/2 on complex, real part
10.2b21.4%Octave 3.8, jcobi/4

localize2.7s (0.2%)

bsearch2.4s (0.1%)

end1.0ms (0.0%)