Details

simplify8.5m (43.2%)

sample4.0m (20.4%)

Algorithm
100×intervals
Results
1.3m15785×body10240exit
55.9s456200×body80valid
33.7s30262×body1280valid
17.3s22978×body640valid
7.3s71216×body80nan
7.2s3120×body1280nan
5.6s12176×body320valid
3.6s2495×body640nan
3.3s9978×body160valid
2.0s49536×pre80true
1.3s1365×body320nan
957.0ms1481×body2560valid
931.0ms611×body5120valid
401.0ms681×body160nan

rewrite3.3m (16.9%)

Algorithm
151×rewrite-expression-head
Calls

465 calls:

1.6m
(* (/ (* (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)) (sqrt (* PI 2.0))) (exp (+ (+ (- z 1.0) 7.0) 0.5))) (+ (+ (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0))) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (+ (/ -176.6150291621406 (+ (- z 1.0) 4.0)) (+ (+ (/ -0.13857109526572012 (+ (- z 1.0) 6.0)) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (+ (/ 676.5203681218851 z) (+ (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) 0.9999999999998099))))))
25.3s
(+ (/ -1259.1392167224028 (+ (- z) 2.0)) (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (/ -176.6150291621406 (+ 4.0 (- z)))) (/ 771.3234287776531 (+ 3.0 (- z)))) (+ (/ 12.507343278686905 (+ 5.0 (- z))) (+ (/ 1.5056327351493116e-07 (+ 8.0 (- z))) (+ (/ 9.984369578019572e-06 (+ (- z) 7.0)) (/ -0.13857109526572012 (+ 6.0 (- z))))))))
24.9s
(* (/ (+ (+ (+ (+ (+ 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)
14.0s
(/ (+ (+ (+ (+ (+ 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)))))
5.3s
(* (- (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0))) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (fma -176.6150291621406 (* (* (- (/ -0.13857109526572012 (+ (- z 1.0) 6.0)) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) z) (fma (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (* 0.9999999999998099 (- 0.9999999999998099 (/ -1259.1392167224028 (+ (- z 1.0) 2.0)))))) (* (+ (- z 1.0) 4.0) (fma (- (* (/ -0.13857109526572012 (+ (- z 1.0) 6.0)) (/ -0.13857109526572012 (+ (- z 1.0) 6.0))) (* (/ 771.3234287776531 (+ (- z 1.0) 3.0)) (/ 771.3234287776531 (+ (- z 1.0) 3.0)))) (* (fma (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (* 0.9999999999998099 (- 0.9999999999998099 (/ -1259.1392167224028 (+ (- z 1.0) 2.0))))) z) (* (- (/ -0.13857109526572012 (+ (- z 1.0) 6.0)) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (fma 676.5203681218851 (fma (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (* 0.9999999999998099 (- 0.9999999999998099 (/ -1259.1392167224028 (+ (- z 1.0) 2.0))))) (* z (+ (pow (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) 3) (pow 0.9999999999998099 3)))))))))
Rules
15650×*-un-lft-identity
12488×times-frac
9174×add-sqr-sqrt
6944×add-cube-cbrt
3555×distribute-lft-out
3030×prod-diff
1963×add-exp-log
1737×sqrt-prod
1600×associate-*r*
1564×add-cbrt-cube
1225×distribute-lft-out--
1204×pow1
1043×distribute-rgt-in distribute-lft-in
932×associate-*l*
906×frac-add
725×difference-of-squares
715×add-log-exp
632×flip-+ flip3-+
600×unpow-prod-down
579×cbrt-prod
572×prod-exp
568×associate-/l*
560×associate-/r/
530×div-inv
489×associate-/r*
468×log1p-expm1-u
467×expm1-log1p-u
433×div-exp
417×cbrt-unprod
313×cbrt-undiv
295×cube-prod
290×pow-prod-down
262×associate-*r/
256×unswap-sqr
237×swap-sqr
211×frac-times
200×sqr-pow
190×fma-def
181×fma-neg
143×distribute-rgt-neg-in
135×exp-sum
126×*-commutative
124×sum-log
118×log-pow
115×associate-*l/
98×frac-2neg clear-num
94×flip3-- flip--
86×pow1/3
82×pow-unpow
80×cbrt-div
78×associate-/l/ pow-exp
77×fma-udef
73×associate--r+
66×sqrt-div
62×sqrt-pow1
61×sub-neg
59×unpow3 cube-mult
52×rec-exp 1-exp
49×associate--l+
47×pow-to-exp
46×pow1/2
44×diff-log
41×log-prod
29×pow-prod-up
27×pow-pow
25×rem-cube-cbrt rem-sqrt-square
23×pow-sqr +-commutative
18×hypot-udef pow-plus
17×div-sub
14×exp-diff pow2
12×cube-div un-div-inv pow-flip associate-+l-
10×exp-neg neg-sub0
inv-pow frac-sub
associate-+r+
associate-+l+ log1p-udef distribute-lft-neg-in
cos-sum log1p-expm1 pow-sub expm1-log1p expm1-udef rem-log-exp associate--l-
sin-sum unpow-prod-up
exp-prod associate--r-
unpow1/3 distribute-neg-in neg-mul-1
difference-cubes sqrt-pow2 sqrt-unprod pow3 sub-div rem-square-sqrt distribute-rgt-out-- rem-cbrt-cube distribute-rgt-neg-out log-div unpow2

prune1.9m (9.6%)

series1.7m (8.7%)

Calls

465 calls:

4.1s
(* (+ (sqrt (* d1 (+ (- d2 d3) d4))) d1) (- (sqrt (* d1 (+ (- d2 d3) d4))) d1))
3.0s
(* (/ i (* (+ (sqrt 1.0) (fma 2.0 i (+ alpha beta))) (/ (fma 2.0 i (+ alpha beta)) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i))))))) (/ (/ (+ (+ alpha beta) i) (/ (fma 2.0 i (+ alpha beta)) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i)))))) (- (fma 2.0 i (+ alpha beta)) (sqrt 1.0))))
2.9s
(fma -1259.1392167224028 (* (* (- 4.0 z) (- 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))) (* (- 3.0 z) (- (- (/ 12.507343278686905 (- 5.0 z)) (/ 1.5056327351493116e-07 (- 8.0 z))) (+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0)))))) (* (+ (- z) 2.0) (fma (fma (fma (- (* 0.9999999999998099 0.9999999999998099) (* (/ 676.5203681218851 (- 1.0 z)) (/ 676.5203681218851 (- 1.0 z)))) (- 4.0 z) (* (- 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) -176.6150291621406)) (- 3.0 z) (* (+ (sqrt 0.9999999999998099) (sqrt (/ 676.5203681218851 (- 1.0 z)))) (* (- (sqrt 0.9999999999998099) (sqrt (/ 676.5203681218851 (- 1.0 z)))) (* (- 4.0 z) 771.3234287776531)))) (- (- (/ 12.507343278686905 (- 5.0 z)) (/ 1.5056327351493116e-07 (- 8.0 z))) (+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0)))) (* (* (- 4.0 z) (- 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))) (* (- 3.0 z) (fma (/ 12.507343278686905 (- 5.0 z)) (/ 12.507343278686905 (- 5.0 z)) (- (* (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0)))) (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))))))))))))
1.5s
(/ (pow (* -2.0 (log u1)) 0.5) 6.0)
1.4s
(- (pow 1.0 3) (log1p (expm1 (pow (* (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (/ (fma (fma (fma (fma 1.061405429 (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) -1.453152027) (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) 1.421413741) (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) -0.284496736) (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) 0.254829592) (exp (pow (fabs x) 2)))) 3))))

regimes8.1s (0.7%)

Accuracy

Total 21.6b remaining (17.3%)

Threshold costs 7.2b (5.8%)

6.7b37.4%_divideComplex, real part
2.9b77.1%Octave 3.8, jcobi/1
2.4b36.1%Octave 3.8, jcobi/4
1.2b51.5%Octave 3.8, jcobi/3
1.2b0%_divideComplex, imaginary part

localize5.0s (0.4%)

bsearch1.0s (0.1%)

end1.0ms (0.0%)