Details

rewrite20.0m (65.6%)

Algorithm
159×rewrite-expression-head
Calls

456 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)))
3.4m
(/ (/ (* (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))))) (hypot (pow (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))) 3/2) (pow (+ (+ (+ (+ 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)))) 3/2))) (fma (* 2.0 0.0001789971) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (fma 0.0008327945 (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3))))) (fma 0.0140005442 (* (pow x 2) (* (pow x 2) (* x (pow x 3)))) (fma 0.0694555761 (pow x 6) (fma 0.2909738639 (pow x 4) (fma (* 0.7715471019 x) x 1.0))))))) (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)))) (+ (+ (+ (+ 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))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (+ (+ (+ (+ 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))))))))
34.6s
(- (* (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (+ (+ (+ (+ 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)))))
7.3s
(* (/ (+ (+ (+ (+ (+ 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)
4.1s
(/ (+ (+ (+ (+ (+ 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)))))
Rules
13211×times-frac
11447×*-un-lft-identity
10060×add-sqr-sqrt
7890×add-cube-cbrt
2617×sqrt-prod
2521×add-exp-log
2266×prod-diff
1940×add-cbrt-cube
1471×associate-*r*
1430×pow1
1366×distribute-lft-out
1098×associate-*l*
827×prod-exp
734×associate-/l*
711×add-log-exp
709×unpow-prod-down
663×div-exp
651×associate-/r/
647×cbrt-prod
616×difference-of-squares
571×cbrt-unprod
563×associate-/r*
534×fma-def
522×cube-prod
515×div-inv
476×distribute-lft-out--
451×log1p-expm1-u
450×expm1-log1p-u
439×cbrt-undiv
357×pow-prod-down
356×flip-+ flip3-+
350×log-pow
337×distribute-rgt-in distribute-lft-in
315×fma-neg
288×unswap-sqr
215×sqr-pow
207×distribute-rgt-neg-in
198×associate--l+
197×exp-prod
195×distribute-lft-neg-in
184×flip3-- flip--
183×associate-*r/
180×sqrt-div
174×sqrt-pow1
151×frac-times
131×*-commutative
129×associate-*l/
122×exp-sum
116×unpow3 cube-mult
113×1-exp log-prod
108×pow-exp
100×cube-div
94×sub-neg
93×associate-/l/ clear-num frac-2neg
91×pow-unpow
88×diff-log
87×pow1/2
76×swap-sqr
75×rec-exp neg-mul-1
73×associate--r+
62×pow1/3
59×pow-to-exp cbrt-div
57×exp-neg
55×pow-prod-up
54×neg-sub0 pow-pow
51×exp-diff
45×rem-sqrt-square
43×fma-udef
42×exp-to-pow
41×pow-sqr
36×sum-log
34×rem-cube-cbrt
25×un-div-inv
24×pow-plus
22×frac-add +-commutative
21×associate--l-
20×hypot-def
18×distribute-neg-in
17×associate-+l+
16×pow2
12×associate-+r+ div-sub
10×inv-pow pow-flip
unpow2
rem-log-exp
rem-exp-log associate-+l- neg-log
pow-div sin-sum expm1-log1p frac-sub expm1-udef
sqrt-unprod distribute-lft-neg-out rem-cbrt-cube
sqrt-undiv difference-cubes sqrt-pow2 log1p-udef pow3 hypot-udef rem-square-sqrt distribute-rgt-neg-out log-div
log1p-expm1 unpow-prod-up pow-sub unsub-neg distribute-rgt-out cos-sum associate--r- distribute-rgt-out--

series6.7m (22.1%)

Calls

450 calls:

5.7m
(- (pow (pow 1.0 3) 3) (pow (log1p (expm1 (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))
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))
2.0s
(/ (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta))))
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))))
987.0ms
(- (pow 1.0 3) (log1p (expm1 (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))))

sample2.1m (6.9%)

Algorithm
98×intervals
Results
44.8s7687×body10240exit
33.5s454135×body80valid
14.1s27267×body1280valid
6.9s21054×body640valid
3.8s68309×body80nan
2.7s1644×body1280nan
2.3s10373×body320valid
1.7s9489×body160valid
1.6s1308×body640nan
969.0ms49536×pre80true
545.0ms1427×body2560valid
521.0ms601×body320nan
259.0ms4306×body80overflowed
247.0ms598×body5120valid
197.0ms344×body160nan

prune1.2m (3.9%)

simplify16.4s (0.9%)

Algorithm
255×egg-herbie

regimes5.6s (0.3%)

Accuracy

Total 44.3b remaining (40.4%)

Threshold costs 5.5b (5.1%)

8.2b34.5%_divideComplex, imaginary part
7.0b36.3%_divideComplex, real part
5.5b14.9%Octave 3.8, jcobi/2
5.1b56.7%Octave 3.8, jcobi/4
4.4b23.1%math.sqrt on complex, imaginary part, im greater than 0 branch

localize2.6s (0.1%)

bsearch640.0ms (0.0%)

end5.0ms (0.0%)