Details

simplify48.8m

Calls

21501 calls:

Slowest
2.1s
(/ (* (* (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i)) (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i))) (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i))) (* (* (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta)))) (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta)))))
1.7s
(/ (* (* (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i)) (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i))) (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i))) (* (* (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) (fma 2 i (+ alpha beta))) (* (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) (fma 2 i (+ alpha beta)))))
1.7s
(/ (* (* (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma (+ (+ alpha beta) i) i (* beta alpha))) (fma (+ (+ alpha beta) i) i (* beta alpha))) (* (* (* (+ (+ alpha beta) i) i) (* (+ (+ alpha beta) i) i)) (* (+ (+ alpha beta) i) i))) (* (* (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) (fma 2 i (+ alpha beta))) (* (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) (fma 2 i (+ alpha beta)))))
1.7s
(/ (* (* (* (/ (/ (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma 2 i (+ alpha beta))) (+ (fma 2 i (+ alpha beta)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (/ (fma 2 i (+ alpha beta)) i))) (* (/ (/ (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma 2 i (+ alpha beta))) (+ (fma 2 i (+ alpha beta)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (/ (fma 2 i (+ alpha beta)) i)))) (* (/ (/ (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma 2 i (+ alpha beta))) (+ (fma 2 i (+ alpha beta)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (/ (fma 2 i (+ alpha beta)) i)))) (* (* (- (fma 2 i (+ alpha beta)) (sqrt 1.0)) (- (fma 2 i (+ alpha beta)) (sqrt 1.0))) (- (fma 2 i (+ alpha beta)) (sqrt 1.0))))
1.6s
(/ (* (* (* (/ (/ (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma 2 i (+ alpha beta))) (+ (fma 2 i (+ alpha beta)) (sqrt 1.0))) (/ (* (+ (+ alpha beta) i) i) (fma 2 i (+ alpha beta)))) (* (/ (/ (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma 2 i (+ alpha beta))) (+ (fma 2 i (+ alpha beta)) (sqrt 1.0))) (/ (* (+ (+ alpha beta) i) i) (fma 2 i (+ alpha beta))))) (* (/ (/ (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma 2 i (+ alpha beta))) (+ (fma 2 i (+ alpha beta)) (sqrt 1.0))) (/ (* (+ (+ alpha beta) i) i) (fma 2 i (+ alpha beta))))) (* (* (- (fma 2 i (+ alpha beta)) (sqrt 1.0)) (- (fma 2 i (+ alpha beta)) (sqrt 1.0))) (- (fma 2 i (+ alpha beta)) (sqrt 1.0))))

sample9.3m

Algorithm
87×intervals
Results
6.4m100267×body10240exit
1.0m379688×body80valid
37.7s25612×body1280valid
19.1s19903×body640valid
12.1s70550×body80nan
7.1s10243×body320valid
6.5s1655×body1280nan
3.0s1358×body640nan
2.5s5383×body160valid
1.7s1422×body2560valid
1.5s33536×pre80true
1.1s650×body320nan
710.0ms582×body5120valid
331.0ms303×body160nan

prune5.3m

rewrite3.9m

Algorithm
144×rewrite-expression-head
Calls

455 calls:

Slowest
1.3m
(* (/ 1 (exp (+ (- (- 1 z) -6) 0.5))) (+ (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (+ (+ (+ (/ -1259.1392167224028 (+ (- 1 z) 1)) 0.9999999999998099) (/ 676.5203681218851 (- 1 z))) (/ -176.6150291621406 (- 5 (+ z 1))))) (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6))))))
33.8s
(cbrt (* (* (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))) (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))))
17.9s
(* (* (+ (+ (+ (+ (+ 1 (* 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 (+ (+ (+ (+ (+ (+ 1 (* 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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))) x)
16.9s
(* (+ (+ (+ (+ (+ 1 (* 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 (+ (+ (+ (+ (+ (+ 1 (* 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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))))
16.0s
(* (/ (+ (+ (+ (+ (+ 1 (* 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.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.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x)
Rules
37580×times-frac
36263×*-un-lft-identity
24306×frac-times
18202×add-sqr-sqrt
16003×associate-*r/
15600×add-cube-cbrt
12386×sqrt-prod
11367×distribute-lft-out
10039×associate-*l/
7684×frac-add
7674×flip-+ flip3-+
7581×exp-neg
6499×exp-diff
6498×neg-sub0
2601×pow1
2541×add-exp-log
1896×add-cbrt-cube
1790×prod-diff
1625×associate-*l*
1090×un-div-inv
1049×div-inv
1037×associate-*r*
930×associate-/r/
891×associate-/l*
839×prod-exp
792×associate-/r*
789×pow-prod-up
732×div-exp
695×difference-of-squares
688×sqrt-pow1
653×unpow-prod-down
638×add-log-exp
614×frac-sub
540×pow-prod-down
529×cbrt-unprod
519×cube-prod
460×log1p-expm1-u
455×expm1-log1p-u insert-posit16 distribute-lft-out--
440×cbrt-undiv
406×unswap-sqr
383×flip3-- flip--
378×fma-def
362×pow1/2
322×log-pow
318×cbrt-prod
293×pow-plus
278×swap-sqr
266×exp-sum
245×pow-sqr
189×fma-neg
182×log-prod
179×pow1/3
164×sqrt-div
150×sqr-pow
135×*-commutative
116×distribute-rgt-in distribute-lft-in
106×1-exp pow2
97×frac-2neg clear-num
96×rec-exp
86×cube-mult
85×unpow3
84×associate-/l/
71×unpow-prod-up
69×associate--r+
67×cbrt-div
63×sub-neg
61×pow-flip
59×fma-udef distribute-rgt-out
58×sum-log distribute-rgt1-in
55×pow-exp
46×associate-+l+
45×pow-to-exp
43×diff-log
38×+-commutative
35×inv-pow
34×rem-sqrt-square exp-prod
32×pow-div
24×cube-div
23×cube-unmult
19×sqrt-unprod div-sub
17×pow3 hypot-def
16×sub-div
15×associate-+r+
13×hypot-udef
12×sin-sum pow-unpow
10×associate--l+ associate-+l-
pow-pow pow-sub
expm1-log1p expm1-udef rem-square-sqrt
rem-cube-cbrt
sum-cubes
associate--r-
difference-cubes rem-log-exp rem-cbrt-cube
associate--l- associate-+r- rem-exp-log cos-sum

series1.6m

Calls

455 calls:

Slowest
4.7s
(- (pow 1 3) (pow (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (* (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x))))) (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x)))))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
3.6s
(+ (* (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (* (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x))))) (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x)))))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (* (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x))))) (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x)))))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1 (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (* (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x))))) (/ 1 (sqrt (+ 1 (* 0.3275911 (fabs x)))))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
1.5s
(/ (pow (* -2 (log u1)) 0.5) 6)
1.4s
(expm1 (expm1 (log1p (/ (pow (* -2 (log u1)) 0.5) 6))))
1.3s
(cbrt (expm1 (/ (pow (* -2 (log u1)) 0.5) 6)))

localize8.1s

regimes7.3s

Accuracy

Total 29.8b remaining (29.9%)

8.5b45.5%Octave 3.8, jcobi/4
7.1b1.6%_divideComplex, real part
7.0b0.9%_divideComplex, imaginary part
3.1b75.4%Octave 3.8, jcobi/1
0.6b0%math.cos on complex, imaginary part

bsearch1.3s

end1.0ms