Details

rewrite10.2m (67.7%)

Algorithm
156×rewrite-expression-head
Calls

469 calls:

779.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))))
477.0ms
(/ (+ (+ (* (* x (* (pow (* x x) 3) (pow x 3))) (+ 0.0008327945 (* (* x x) (* 2.0 0.0001789971)))) (fma (* 0.2909738639 x) (pow x 3) (fma (* 0.7715471019 x) x 1.0))) (* (pow x 6) (+ 0.0694555761 (* (* x x) 0.0140005442)))) (+ (+ (* (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971))) (fma (* 0.1049934947 x) x 1.0)) (* (pow x 4) (+ 0.0424060604 (* (* x x) 0.0072644182)))))
467.0ms
(* (/ x (+ (+ (* (* x (* (pow (* x x) 3) (pow x 3))) (+ 0.0008327945 (* (* x x) (* 2.0 0.0001789971)))) (fma (* 0.2909738639 x) (pow x 3) (fma (* 0.7715471019 x) x 1.0))) (* (pow x 6) (+ 0.0694555761 (* (* x x) 0.0140005442))))) (fma (fma (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971)) (fma (* 0.1049934947 x) x 1.0)) (fma (pow (* x x) 4) (+ 0.0005064034 (* (* x x) 0.0001789971)) (fma (* 0.1049934947 x) x 1.0)) (- (* (* (pow x 4) (+ 0.0424060604 (* (* x x) 0.0072644182))) (* (pow x 4) (+ 0.0424060604 (* (* x x) 0.0072644182)))))))
404.0ms
(- (+ (/ (/ x (+ wj 1)) (exp wj)) wj) (/ wj (+ wj 1)))
358.0ms
(- (/ 1 (/ (+ (+ alpha beta) 2.0) beta)) (- (/ alpha (+ (+ alpha beta) 2.0)) 1.0))
Rules
8413×*-un-lft-identity
7073×add-sqr-sqrt
6819×times-frac
5511×add-cube-cbrt
2780×add-exp-log
2763×prod-diff
1718×add-cbrt-cube
1596×pow1
1270×distribute-lft-out
1124×sqrt-prod
1052×prod-exp
863×unpow-prod-down
833×associate-*l*
796×associate-*r*
764×add-log-exp
759×div-exp
748×cube-prod
533×associate-/l*
483×cbrt-unprod
482×associate-/r*
467×log1p-expm1-u expm1-log1p-u
437×associate-/r/
433×difference-of-squares
390×log-pow
382×fma-neg
374×div-inv
359×cbrt-prod
351×pow-prod-down
328×distribute-lft-out--
326×distribute-rgt-in distribute-lft-in
321×cbrt-undiv
299×flip-+ flip3-+
289×pow-exp
248×1-exp
224×unswap-sqr
203×exp-sum
201×associate--l+
188×rec-exp associate-*r/
187×log-prod
183×sqr-pow
168×swap-sqr
161×distribute-rgt-neg-in
153×distribute-lft-neg-in
150×pow-to-exp
140×fma-def associate--r+ flip3-- flip--
138×unpow3 cube-mult
134×exp-prod
130×*-commutative
112×pow-prod-up
106×sub-neg
99×frac-2neg clear-num diff-log
85×sqrt-div associate-*l/
79×pow-unpow
77×associate-/l/ pow-plus
76×sqrt-pow1
72×pow1/3
61×fma-udef
60×pow-sqr
59×frac-times
55×cbrt-div
54×pow1/2
44×neg-mul-1 exp-to-pow
42×sum-log
38×frac-add pow2
36×pow-pow
32×pow-flip
27×+-commutative
26×inv-pow
25×rem-sqrt-square
23×distribute-neg-in
22×div-sub
20×associate--l-
17×rem-cube-cbrt
16×rem-log-exp
14×associate-+l+ frac-sub
12×hypot-def
11×rem-exp-log hypot-udef associate-+r+
neg-log log-div
sqrt-undiv
cube-div exp-diff
neg-sub0 distribute-rgt-neg-out
pow3 sin-sum cos-sum unpow2
sqrt-unprod rem-square-sqrt associate-+l-
log-rec associate--r-
difference-cubes sqrt-pow2 un-div-inv exp-neg unpow-prod-up pow-sub unsub-neg distribute-rgt-out-- rem-cbrt-cube

sample2.4m (16.0%)

Algorithm
99×intervals
Results
45.9s453510×body80valid
45.0s7720×body10240exit
14.7s26686×body1280valid
7.4s20994×body640valid
5.7s70823×body80nan
2.6s1509×body1280nan
2.6s10495×body320valid
2.1s9561×body160valid
1.6s1249×body640nan
1.1s49536×pre80true
676.0ms1503×body2560valid
540.0ms563×body320nan
359.0ms4180×body80overflowed
295.0ms609×body5120valid
223.0ms333×body160nan

series1.0m (7.0%)

Calls

467 calls:

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))))
1.1s
(- 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (log (exp (/ 1.0 (+ (* 0.3275911 (fabs x)) 1.0)))) (+ 1.421413741 (+ (/ (* -1.453152027 1.0) (fma 0.3275911 (fabs x) 1.0)) (/ (* (/ (* 1.061405429 1.0) (fma 0.3275911 (fabs x) 1.0)) 1.0) (fma 0.3275911 (fabs x) 1.0))))))))) (exp (- (* (fabs x) (fabs x))))))
949.0ms
(exp (log (- 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.453152027 1.0) (fma 0.3275911 (fabs x) 1.0)) (/ (* (/ (* 1.061405429 1.0) (fma 0.3275911 (fabs x) 1.0)) 1.0) (fma 0.3275911 (fabs x) 1.0))))))))) (exp (- (* (fabs x) (fabs x))))))))
946.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.453152027 1.0) (fma 0.3275911 (fabs x) 1.0)) (/ (* (/ (* 1.061405429 1.0) (fma 0.3275911 (fabs x) 1.0)) 1.0) (fma 0.3275911 (fabs x) 1.0))))))))) (exp (- (* (fabs x) (fabs x))))))

prune59.9s (6.6%)

simplify15.8s (1.8%)

Algorithm
254×egg-herbie

regimes5.3s (0.6%)

Accuracy

Total 49.9b remaining (45.4%)

Threshold costs 5.1b (4.6%)

10.0b14.4%_divideComplex, imaginary part
9.9b41.1%Octave 3.8, jcobi/4
9.0b17.7%_divideComplex, real part
4.9b10.7%math.sqrt on complex, imaginary part, im greater than 0 branch
4.2b34.9%math.sqrt on complex, real part

localize2.5s (0.3%)

bsearch967.0ms (0.1%)

end1.0ms (0.0%)