Details

simplify34.1m

Calls

11023 calls:

Slowest
2.5s
(/ (/ (* (* (* (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))))) (* (* (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0) (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0)) (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0)))
2.3s
(/ (/ (* (* (* (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))))) (* (* (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0) (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0)) (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0)))
1.8s
(* 1 (* (* (pow (+ (- z -6) 0.5) (+ (- z 1) 0.5)) (sqrt (* PI 2))) (+ (pow (+ (+ (+ (/ 771.3234287776531 (+ z 2)) (/ -1259.1392167224028 (+ z 1))) (+ 0.9999999999998099 (/ 676.5203681218851 z))) (+ (/ -176.6150291621406 (+ z 3)) (+ (/ -0.13857109526572012 (- z -5)) (/ 12.507343278686905 (+ z 4))))) 3) (pow (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7))) 3))))
1.6s
(/ (* (* (* (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.5s
(/ (/ (* (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (fma (+ (+ alpha beta) i) i (* beta alpha))) (fma (+ (+ alpha beta) i) i (* beta alpha))) (* (* (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.0) (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0)) (- (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) 1.0)) (/ (* (* (* (+ (+ 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))))))

sample8.3m

Algorithm
87×intervals
Results
6.0m96143×body10240exit
1.0m385877×body80valid
26.3s22572×body1280valid
13.7s18090×body640valid
10.0s68798×body80nan
4.9s9575×body320valid
2.2s41536×pre80true
1.8s5197×body160valid
1.8s1424×body2560valid
597.0ms624×body5120valid
495.0ms41×body320nan
239.0ms96×body1280nan
194.0ms82×body640nan
20.0ms22×body160nan

prune3.4m

series1.3m

Calls

414 calls:

Slowest
4.5s
(- (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 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
3.5s
(+ (* (* (* (/ 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 (+ 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 (+ 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 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
1.3s
(* 1/6 (pow (* -2 (log u1)) 0.5))
1.1s
(log (* 1/6 (pow (* -2 (log u1)) 0.5)))
1.1s
(- 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 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))))

rewrite37.9s

Algorithm
139×rewrite-expression-head
Calls

414 calls:

Slowest
2.7s
(* (/ 1 (sqrt (+ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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))
2.5s
(/ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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))))))
2.5s
(/ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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))))))
2.4s
(/ (+ (+ (+ (+ (+ 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)))) (sqrt (+ (+ (+ (+ (+ (+ 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))))))
2.4s
(* (/ (+ (+ (+ (+ (+ 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
27622×frac-times
18602×associate-*r/
12586×frac-add
10948×associate-*l/
9716×flip-+ flip3-+
8099×exp-neg
7020×exp-diff
7015×neg-sub0
6216×times-frac
3867×*-un-lft-identity
3768×add-sqr-sqrt
3725×add-cube-cbrt
2093×cbrt-div
1175×add-exp-log
1099×un-div-inv
1020×add-cbrt-cube
920×pow1
879×flip3-- flip--
777×associate-+l-
776×pow-sub
687×sqrt-prod
611×div-inv
486×add-log-exp
485×associate-*l*
423×associate-*r*
417×associate-/r*
416×log1p-expm1-u associate-/l*
414×expm1-log1p-u insert-posit16
312×prod-exp
304×associate-/r/
224×cbrt-unprod
220×pow-prod-down
210×cbrt-prod
190×div-exp
165×cbrt-undiv
136×sqrt-div
133×prod-diff
132×*-commutative
96×fma-def
93×distribute-lft-out associate-/l/
84×frac-2neg clear-num
80×fma-neg
79×pow1/3
68×sub-neg
63×distribute-rgt-in distribute-lft-in
62×log-pow
60×rec-exp
59×pow1/2
57×pow-prod-up
56×log-prod
53×fma-udef
42×pow-div
41×pow-plus
34×exp-sum
32×unpow-prod-down
29×+-commutative
28×pow-flip pow2
27×inv-pow exp-prod
26×frac-sub rem-sqrt-square
23×sum-log
19×difference-of-squares
18×cube-div
17×div-sub
15×cube-mult diff-log
14×unpow3
12×associate--l+
11×sin-sum
10×pow-exp
associate-+r+ distribute-lft-out-- associate--r+ log-div
associate-+l+ rem-log-exp distribute-rgt-out
pow-to-exp
cube-prod hypot-def
distribute-rgt1-in pow3 rem-exp-log
cube-unmult difference-of-sqr-1 sqrt-unprod associate--l- pow-pow expm1-def
rem-cube-cbrt sqrt-undiv rem-square-sqrt associate--r-
unpow-prod-up expm1-log1p expm1-udef log1p-udef
distribute-lft1-in log1p-expm1 associate-+r- difference-cubes distribute-rgt-out-- cos-sum rem-cbrt-cube distribute-neg-in sin-mult

regimes7.1s

Accuracy

Total 36.1b remaining (33%)

10.0b25.7%_divideComplex, real part
8.6b18.6%_divideComplex, imaginary part
5.7b8.4%Octave 3.8, jcobi/2
2.9b77.4%Octave 3.8, jcobi/1
2.2b23.2%math.cos on complex, imaginary part

localize3.8s

bsearch1.4s

end1.0ms