Details

rewrite10.9m (31.3%)

Algorithm
160×rewrite-expression-head
Calls

479 calls:

1.3m
(+ (+ (+ (/ 1.5056327351493116e-07 (- 8 z)) (/ 9.984369578019572e-06 (- 7 z))) (+ (+ (/ -0.13857109526572012 (- 6 z)) (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (/ -176.6150291621406 (+ 1 (+ (- 1 z) 2))))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1 z)) (/ -1259.1392167224028 (- 2 z)))))) (/ 12.507343278686905 (- 5 z)))
1.3m
(+ (+ (+ (/ 1.5056327351493116e-07 (- 8 z)) (/ 9.984369578019572e-06 (- 7 z))) (+ (+ (/ -0.13857109526572012 (- 6 z)) (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (/ -176.6150291621406 (+ 1 (+ (- 1 z) 2))))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1 z)) (/ -1259.1392167224028 (- 2 z)))))) (/ 12.507343278686905 (- 5 z)))
1.1m
(* (* (/ (sqrt (* PI 2)) (exp (+ 0.5 z))) (* (exp -6) (+ (+ (+ (+ (/ 771.3234287776531 (+ z 2)) (/ -1259.1392167224028 (+ 1 z))) (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -176.6150291621406 (+ 3 z)))) (+ (+ (/ -0.13857109526572012 (+ z 5)) (/ 12.507343278686905 (- z -4))) (/ 1.5056327351493116e-07 (+ 7 z)))) (/ 9.984369578019572e-06 (+ 6 z))))) (pow (+ 0.5 (- z -6)) (- z (- 1 0.5))))
51.9s
(+ (+ (+ (/ 1.5056327351493116e-07 (- 8 z)) (/ 9.984369578019572e-06 (- 7 z))) (+ (+ (/ -0.13857109526572012 (- 6 z)) (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (/ -176.6150291621406 (+ 1 (+ (- 1 z) 2))))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1 z)) (/ -1259.1392167224028 (- 2 z)))))) (/ 12.507343278686905 (- 5 z)))
49.6s
(+ (+ (+ (/ 1.5056327351493116e-07 (- 8 z)) (/ 9.984369578019572e-06 (- 7 z))) (+ (+ (/ -0.13857109526572012 (- 6 z)) (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (/ -176.6150291621406 (+ 1 (+ (- 1 z) 2))))) (+ 0.9999999999998099 (+ (/ 676.5203681218851 (- 1 z)) (/ -1259.1392167224028 (- 2 z)))))) (/ 12.507343278686905 (- 5 z)))
Rules
25184×frac-times
17866×associate-*r/
14184×frac-add
10337×associate-*l/
9274×flip-+ flip3-+
8519×*-un-lft-identity
7581×exp-neg
7308×times-frac
6502×exp-diff
6498×neg-sub0
6475×add-sqr-sqrt
5040×add-cube-cbrt
4889×pow1
3657×add-exp-log
2308×add-cbrt-cube
2080×prod-diff
1938×distribute-lft-out
1854×pow-prod-up
1662×prod-exp
1387×sqrt-prod
1140×add-log-exp
1103×un-div-inv
1089×associate-*l*
1019×associate-*r*
1000×pow-plus
942×div-exp
875×pow-prod-down
842×cbrt-unprod
681×unpow-prod-down
605×associate-/l*
501×associate-/r*
498×associate-/r/
481×expm1-log1p-u
479×insert-posit16 log1p-expm1-u
477×cube-prod
476×cbrt-div
469×sum-log
446×pow-sqr
437×sqr-pow
427×difference-of-squares
419×cbrt-undiv
393×associate-/l/
371×div-inv
313×distribute-lft-out--
310×pow2
307×pow-sub
293×cube-mult
292×log-pow
291×unpow3
286×flip3-- flip--
272×cbrt-prod
267×pow-exp
256×unswap-sqr
253×fma-def
241×distribute-rgt-in distribute-lft-in
188×swap-sqr
170×exp-sum
164×pow-to-exp
155×fma-neg
152×*-commutative
151×log-prod
135×unpow-prod-up
133×1-exp rec-exp
119×sqrt-pow1
114×cube-div
104×frac-2neg clear-num
100×pow1/3
78×pow1/2
76×sqrt-div
67×sub-neg
61×fma-udef
52×pow-flip
45×cube-unmult
41×inv-pow
38×diff-log
35×sin-sum
31×+-commutative
29×pow3 rem-sqrt-square
28×distribute-rgt-out
26×cos-sum
25×distribute-rgt1-in pow-unpow
23×div-sub
22×frac-sub
20×exp-prod
16×associate-+l+
11×pow-pow
10×associate--r+
8×pow-div
6×distribute-lft-neg-in distribute-rgt-neg-in hypot-udef rem-log-exp
5×hypot-def
4×rem-cube-cbrt associate--l+ sqrt-unprod associate-+l- associate--r-
3×associate-+r- rem-square-sqrt log-div
2×expm1-log1p expm1-udef associate-+r+ distribute-rgt-neg-out log1p-expm1 sub0-neg neg-mul-1 distribute-rgt-out-- rem-exp-log rem-cbrt-cube distribute-neg-in log1p-udef count-2
1×associate--l- difference-cubes unsub-neg

simplify9.7m (27.8%)

prune4.7m (13.5%)

series3.9m (11.2%)

Calls

478 calls:

18.4s
(fma 1 (* 1 1) (- (* (/ (pow (exp (- (* (fabs x) (fabs x)))) 3) (* (+ 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)))))))))) (/ (pow (* 1 (- (* 0.254829592 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))))))) (* (/ 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)))))))))) 3) (* (* (+ 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))))))))) (* (+ 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))))))))))))))
18.1s
(- (pow 1 3) (/ (pow (* (* 1 (- (* 0.254829592 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))))))) (* (/ 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) (pow (* (+ 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))))))))) 3)))
4.3s
(- (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.8s
(fma (- (/ (pow (exp (- (* (fabs x) (fabs x)))) 3) (* (+ 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))))))))))) (/ (pow (* 1 (- (* 0.254829592 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))))))) (* (/ 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)))))))))) 3) (* (* (+ 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))))))))) (* (+ 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))))))))))) (* (/ (pow (exp (- (* (fabs x) (fabs x)))) 3) (* (+ 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)))))))))) (/ (pow (* 1 (- (* 0.254829592 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))))))) (* (/ 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)))))))))) 3) (* (* (+ 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))))))))) (* (+ 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)))))))))))))
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)))))))

sample3.1m (9.0%)

Algorithm
99×intervals
Results
54.7s12116×body10240exit
53.0s450374×body80valid
26.2s26667×body1280valid
14.1s20987×body640valid
7.7s70285×body80nan
3.8s10690×body320valid
3.0s9616×body160valid
2.8s49536×pre80true
2.8s1506×body1280nan
1.5s1318×body640nan
1.0s1464×body2560valid
907.0ms672×body320nan
486.0ms345×body160nan
380.0ms603×body5120valid

end2.2m (6.3%)

regimes10.3s (0.5%)

Accuracy

Total 28.2b remaining (20.7%)

Threshold costs 6.8b (5%)

7.0b35.4%_divideComplex, real part
6.8b33.2%_divideComplex, imaginary part
2.9b77.4%Octave 3.8, jcobi/1
1.9b29.5%Octave 3.8, jcobi/4
1.5b51.6%Octave 3.8, jcobi/3

localize7.8s (0.4%)

bsearch1.2s (0.1%)