Details

rewrite10.5m (68.2%)

Algorithm
153×rewrite-expression-head
Calls

453 calls:

9.9m
(/ (* (* (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)))
7.2s
(* (/ (+ (+ (+ (+ (+ 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)))) (+ (+ (+ (+ (+ (+ 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) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x)
5.1s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
4.2s
(/ (+ (+ (+ (+ (+ 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)))) (+ (+ (+ (+ (+ (+ 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) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))
849.0ms
(* (+ (/ 1.0 (/ 2.0 (pow (/ 1.0 (fabs x)) 3))) (/ (* (+ (fma 1.875 (/ 1 (pow (fabs x) 6)) (* 0.75 (/ 1 (pow (fabs x) 4)))) 1) 1.0) (fabs x))) (* (/ 1.0 (sqrt PI)) (pow (exp (fabs x)) (/ (fabs x) 2))))
Rules
14126×times-frac
13598×*-un-lft-identity
9295×add-sqr-sqrt
7284×add-cube-cbrt
2660×distribute-lft-out
2612×add-exp-log
2515×cbrt-prod
1740×add-cbrt-cube
1504×prod-diff
1415×sqrt-prod
1349×pow1
1043×associate-*r*
931×prod-exp
798×associate-*l*
756×associate-/r/
697×div-exp
665×add-log-exp difference-of-squares
630×distribute-lft-out--
576×associate-/l*
537×div-inv
523×cbrt-unprod
487×fma-def
454×log1p-expm1-u
450×expm1-log1p-u
415×associate-/r*
343×pow-prod-down
339×cbrt-undiv
322×flip-+ flip3-+
282×flip3-- flip--
250×fma-neg
230×unpow-prod-down
229×unswap-sqr
195×log-pow
167×associate-*r/
165×exp-sum
164×1-exp
157×distribute-rgt-in distribute-lft-in
142×sqrt-pow1
136×associate--r+
129×associate--l+
126×associate-*l/
123×sqrt-div
120×cbrt-div
115×*-commutative
114×pow-exp
113×log-prod
109×swap-sqr
104×pow1/3
89×clear-num frac-2neg
88×sub-neg
84×associate-/l/
83×cube-prod
82×rec-exp
81×diff-log
78×frac-times
74×pow-unpow sqr-pow
72×pow-to-exp
66×pow1/2
60×fma-udef
47×pow-prod-up
44×pow-flip
38×rem-sqrt-square frac-add sum-log
37×pow-sqr
33×distribute-rgt-neg-in distribute-lft-neg-in
29×exp-prod
25×neg-mul-1
24×+-commutative
23×inv-pow
21×pow-pow
17×pow-plus associate-+l+
16×pow2
14×hypot-def associate--l-
13×unpow3 cube-mult
12×rem-cube-cbrt div-sub
11×log-div
10×frac-sub
exp-diff
rem-log-exp
rem-exp-log
un-div-inv hypot-udef associate-+l-
cube-div unpow2 expm1-log1p expm1-udef
sum-cubes log1p-udef sin-sum log1p-expm1
distribute-neg-in neg-sub0 neg-log
sqrt-unprod rem-square-sqrt cos-sum associate--r- distribute-rgt-neg-out
sqrt-pow2 distribute-lft-neg-out sqrt-undiv unpow-prod-up pow-sub associate-+r+

sample1.8m (11.4%)

Algorithm
99×intervals
Results
42.7s11830×body10240exit
22.7s482608×body80valid
11.5s27380×body1280valid
5.2s20759×body640valid
2.7s69699×body80nan
1.7s10743×body320valid
1.6s1552×body1280nan
1.1s9447×body160valid
926.0ms1326×body640nan
677.0ms49536×pre80true
400.0ms1536×body2560valid
309.0ms616×body320nan
158.0ms598×body5120valid
108.0ms322×body160nan

series1.6m (10.4%)

Calls

450 calls:

17.7s
(log (fma (* (cbrt (/ (/ (+ alpha beta) (/ (fma i 2.0 (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (cbrt (/ (/ (+ alpha beta) (/ (fma i 2.0 (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)))) (cbrt (/ (/ (+ alpha beta) (/ (fma i 2.0 (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) 1.0))
13.9s
(fma (* (cbrt (/ (/ (+ alpha beta) (/ (fma i 2.0 (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) (cbrt (/ (/ (+ alpha beta) (/ (fma i 2.0 (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)))) (cbrt (/ (/ (+ alpha beta) (/ (fma i 2.0 (+ alpha beta)) (- beta alpha))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))) 1.0)
1.9s
(/ (* (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (fma i 2.0 (+ alpha beta))) (fma beta alpha (* i (+ (+ alpha beta) i))))
1.5s
(/ (+ (pow (* (* 1.0 (/ (pow (* -2.0 (log u1)) 0.5) 6.0)) (cos (* (* 2.0 PI) u2))) 3) (pow 0.5 3)) (fma (* (* 1.0 (/ (pow (* -2.0 (log u1)) 0.5) 6.0)) (* 1.0 (/ (pow (* -2.0 (log u1)) 0.5) 6.0))) (* (cos (* (* 2.0 PI) u2)) (cos (* (* 2.0 PI) u2))) (* 0.5 (- 0.5 (* (* 1.0 (/ (pow (* -2.0 (log u1)) 0.5) 6.0)) (cos (* (* 2.0 PI) u2)))))))
1.1s
(/ (/ (* i (/ (+ (+ alpha beta) i) (fma i 2.0 (+ alpha beta)))) (+ (pow (+ (+ alpha beta) (* 2.0 i)) 3) (* (fma i 2.0 (+ alpha beta)) (- 1.0)))) (/ 1 (fma beta alpha (* i (+ (+ alpha beta) i)))))

prune1.2m (7.5%)

simplify14.7s (1.6%)

Algorithm
251×egg-herbie

regimes5.4s (0.6%)

Accuracy

Total 44.6b remaining (30%)

Threshold costs 6.3b (4.3%)

8.5b31.7%_divideComplex, imaginary part
7.6b31%_divideComplex, real part
5.5b3.3%math.sqrt on complex, real part
5.5b18.5%Octave 3.8, jcobi/2
4.1b29.8%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.6s (0.3%)

bsearch850.0ms (0.1%)

end1.0ms (0.0%)