Details

prune6.4m (28.0%)

rewrite6.1m (26.7%)

Algorithm
153×rewrite-expression-head
Calls

453 calls:

2.4m
(* (/ (* (pow (+ (- z -6) 0.5) (+ (- z 1) 0.5)) (sqrt (* PI 2))) (exp (+ (- z -6) 0.5))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7)))))
1.3m
(* (* (/ (pow (+ (- z -6) 0.5) (+ (- z 1) 0.5)) (sqrt (exp (+ (- z -6) 0.5)))) (/ (sqrt (* PI 2)) (sqrt (exp (+ (- z -6) 0.5))))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7)))))
1.3m
(* (* (pow (+ (- z -6) 0.5) (+ (- z 1) 0.5)) (sqrt (* PI 2))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7)))))
47.7s
(/ (* (* (pow (+ (- z -6) 0.5) (+ (- z 1) 0.5)) (sqrt (* PI 2))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7))))) (* (sqrt (exp (+ (- z -6) 0.5))) (sqrt (exp (+ (- z -6) 0.5)))))
2.5s
(* (* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1)) d1)
Rules
35167×*-un-lft-identity
30537×times-frac
13989×distribute-lft-out
13378×add-sqr-sqrt
12088×add-cube-cbrt
4047×pow1
3229×sqrt-prod
3007×add-exp-log
1857×add-cbrt-cube
1724×associate-/r*
1540×pow-prod-up
1432×frac-add
1160×prod-diff
1064×prod-exp
1022×associate-/r/
963×associate-*r*
904×flip-+ flip3-+
879×div-exp
861×associate-*l*
810×pow-plus
713×pow-prod-down
702×div-inv
650×cbrt-prod
622×cbrt-unprod
597×add-log-exp
595×associate-/l*
538×sqrt-pow1
454×log1p-expm1-u
453×expm1-log1p-u insert-posit16
419×sqrt-div
370×frac-times pow-sqr
347×cbrt-undiv associate-*r/
331×unswap-sqr
324×associate-*l/
269×log-pow
258×pow2
242×difference-of-squares
231×distribute-rgt-in distribute-lft-in
205×swap-sqr
188×1-exp rec-exp
185×distribute-lft-out--
174×fma-neg
173×pow1/3
161×exp-sum
157×associate-/l/
150×*-commutative
143×flip3-- flip--
134×unpow-prod-down
129×cbrt-div
119×log-prod
115×pow1/2
96×frac-2neg clear-num
93×pow-flip
92×associate--l+ associate-+l-
81×pow-sub fma-def
75×associate--r+ sub-neg
65×sqr-pow
60×cube-unmult inv-pow
56×pow-exp
53×fma-udef
46×diff-log
41×associate-+r+
38×pow-div pow3 pow-to-exp
35×distribute-rgt-out
31×sum-log
26×distribute-rgt1-in rem-sqrt-square
25×+-commutative exp-prod
24×frac-sub
13×exp-diff div-sub
11×rem-log-exp
10×un-div-inv
unpow-prod-up pow-unpow count-2
hypot-def
pow-pow
associate-+l+
log-div
associate--l- sin-sum difference-of-sqr-1 hypot-udef cos-sum expm1-def
sqrt-unprod rem-exp-log rem-square-sqrt
sqrt-undiv rem-cbrt-cube associate--r-
expm1-log1p expm1-udef associate-+r- log1p-expm1 sqrt-pow2 exp-to-pow unpow1/3 log1p-udef

sample5.0m (21.8%)

Algorithm
100×intervals
Results
1.6m16272×body10240exit
1.3m448603×body80valid
41.2s30076×body1280valid
20.9s23829×body640valid
10.0s3119×body1280nan
9.6s69965×body80nan
7.7s12099×body320valid
5.8s2645×body640nan
5.0s10144×body160valid
2.6s49536×pre80true
2.0s1328×body320nan
1.3s1443×body2560valid
657.0ms657×body160nan
631.0ms618×body5120valid

simplify3.8m (16.5%)

series1.3m (5.7%)

Calls

453 calls:

1.0s
(log (sqrt (exp (- 1 (* (fma (fma (+ (/ (+ -1.453152027 (/ 1.061405429 (fma 0.3275911 (fabs x) 1))) (fma 0.3275911 (fabs x) 1)) 1.421413741) (/ 1 (fma 0.3275911 (fabs x) 1)) -0.284496736) (/ 1 (fma 0.3275911 (fabs x) 1)) 0.254829592) (/ (/ 1 (fma 0.3275911 (fabs x) 1)) (exp (* (fabs x) (fabs x)))))))))
1.0s
(/ (pow (* -2 (log u1)) 0.5) 6)
1.0s
(* (/ 1 6) (pow (* -2 (log u1)) 0.5))
986.0ms
(* (/ (pow (* -2 (log u1)) (/ 0.5 2)) 6) (pow (* -2 (log u1)) (/ 0.5 2)))
969.0ms
(log (sqrt (exp (- 1 (* (fma (fma (+ (/ (+ -1.453152027 (/ 1.061405429 (fma 0.3275911 (fabs x) 1))) (fma 0.3275911 (fabs x) 1)) 1.421413741) (/ 1 (fma 0.3275911 (fabs x) 1)) -0.284496736) (/ 1 (fma 0.3275911 (fabs x) 1)) 0.254829592) (/ (/ 1 (fma 0.3275911 (fabs x) 1)) (exp (* (fabs x) (fabs x)))))))))

regimes9.2s (0.7%)

Accuracy

Total 26.2b remaining (17.5%)

Threshold costs 6.1b (4.1%)

7.9b31.5%_divideComplex, real part
3.4b73.1%Octave 3.8, jcobi/1
2.8b0%Octave 3.8, jcobi/4
2.6b0%_divideComplex, imaginary part
1.4b49.4%Octave 3.8, jcobi/3

localize7.6s (0.6%)

bsearch1.2s (0.1%)

end1.0ms (0.0%)