Details

simplify11.6m (34.1%)

rewrite10.2m (29.8%)

Algorithm
140×rewrite-expression-head
Calls

403 calls:

5.0m
(* (* (pow (+ (+ (- z 1.0) 7.0) 0.5) (- z 1.0)) (+ (+ 0.9999999999998099 (+ (+ (/ 676.5203681218851 z) (+ (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (+ (/ 771.3234287776531 (+ (- z 1.0) 3.0)) (/ -176.6150291621406 (+ (- z 1.0) 4.0))))) (+ (/ -0.13857109526572012 (+ (- z 1.0) 6.0)) (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0))))) (+ (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)) (/ 12.507343278686905 (+ (- z 1.0) 5.0))))) (/ (pow (+ (+ (- z 1.0) 7.0) 0.5) 0.5) (exp (+ (+ (- z 1.0) 7.0) 0.5))))
3.9m
(* (+ (+ (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (+ 0.9999999999998099 (+ (+ (/ 676.5203681218851 z) (+ (/ -1259.1392167224028 (+ (- z 1.0) 2.0)) (+ (/ 771.3234287776531 (+ (- z 1.0) 3.0)) (/ -176.6150291621406 (+ (- z 1.0) 4.0))))) (+ (/ -0.13857109526572012 (+ (- z 1.0) 6.0)) (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)))))) (/ (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)) (exp (+ (+ (- z 1.0) 7.0) 0.5))))
897.0ms
(/ (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)) (exp (+ (+ (- z 1.0) 7.0) 0.5)))
745.0ms
(* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1))
665.0ms
(+ (* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
Rules
67458×times-frac
42359×*-un-lft-identity
27193×add-sqr-sqrt
25912×add-cube-cbrt
24575×frac-times
24528×sqrt-prod
15994×associate-*r/
9577×associate-*l/
8068×distribute-lft-out
7882×exp-neg
6962×flip-+ flip3-+
6758×exp-diff
6756×neg-sub0
6529×frac-add
2892×pow1
2503×add-exp-log
2101×div-inv
1906×add-cbrt-cube
1575×prod-diff
1152×prod-exp
1130×un-div-inv
883×associate-*r*
774×associate-*l*
751×pow-prod-down
721×cbrt-unprod
713×pow-prod-up
614×add-log-exp
567×associate-/r*
485×div-exp
405×associate-/l*
401×log1p-expm1-u
400×expm1-log1p-u
394×difference-of-squares
348×unpow-prod-down
340×pow-plus
323×cbrt-undiv
280×cube-prod
279×distribute-lft-out--
269×cbrt-prod
266×fma-neg
257×pow-sqr
232×associate-/r/
204×unswap-sqr
168×associate--l+
151×distribute-rgt-in distribute-lft-in
150×pow2
145×flip3-- flip--
140×sqr-pow
135×exp-sum
125×log-pow
121×*-commutative
115×sqrt-div
108×cube-div
104×fma-def
101×pow-exp pow1/3
87×swap-sqr
86×unpow3 cube-mult
83×sub-neg
77×diff-log
72×associate--r+
68×frac-2neg clear-num
66×cbrt-div associate-/l/
65×fma-udef
57×sqrt-pow1
47×log-prod pow-to-exp
43×sum-log
41×pow1/2
35×1-exp
32×rec-exp
31×+-commutative
27×pow-unpow
25×cube-unmult distribute-rgt-out
24×rem-cube-cbrt pow-pow
22×frac-sub rem-sqrt-square
21×exp-prod
15×associate--l-
14×distribute-rgt1-in
10×div-sub hypot-def count-2
unpow-prod-up
associate-+l+ associate-+r+ pow3 associate-+l-
rem-log-exp
pow-flip sin-sum inv-pow cos-sum
sqrt-unprod
sqrt-undiv difference-cubes hypot-udef rem-square-sqrt associate--r-
expm1-log1p expm1-udef pow-sub associate-+r- rem-exp-log rem-cbrt-cube distribute-rgt-out--

prune4.9m (14.4%)

series4.4m (12.8%)

Calls

400 calls:

2.7m
(- (pow (pow 1.0 3) 3) (pow (* (sqrt (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3)) (sqrt (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))) 3))
4.6s
(- (pow 1.0 3) (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
4.5s
(fma 1.0 (* 1.0 1.0) (* (* (/ (- beta alpha) (fma 2.0 i (+ alpha beta))) (/ (+ alpha beta) (+ 2.0 (fma 2.0 i (+ alpha beta))))) (* (* (/ (- beta alpha) (fma 2.0 i (+ alpha beta))) (/ (+ alpha beta) (+ 2.0 (fma 2.0 i (+ alpha beta))))) (* (/ (- beta alpha) (fma 2.0 i (+ alpha beta))) (/ (+ alpha beta) (+ 2.0 (fma 2.0 i (+ alpha beta))))))))
3.9s
(+ (* (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* 1.0 (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))))
3.3s
(- (pow 1.0 3) (* (sqrt (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3)) (sqrt (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))))

sample2.8m (8.1%)

Algorithm
87×intervals
Results
1.2m12017×body10240exit
36.5s342447×body80valid
13.0s13244×body1280valid
6.7s69450×body80nan
6.3s10911×body640valid
4.8s1586×body1280nan
3.1s1315×body640nan
1.9s5852×body320valid
1.9s49536×pre80true
1.8s7259×body160valid
1.3s718×body320nan
964.0ms749×body2560valid
263.0ms327×body160nan
236.0ms335×body5120valid

localize8.5s (0.4%)

regimes7.2s (0.3%)

Accuracy

Total 19.5b remaining (21.2%)

Threshold costs 2.9b (3.2%)

5.7b20.6%Octave 3.8, jcobi/2
5.1b44.1%Octave 3.8, jcobi/4
2.9b77.6%Octave 3.8, jcobi/1
1.7b0%Jmat.Real.lambertw, newton loop step
1.2b55.2%Octave 3.8, jcobi/3

bsearch477.0ms (0.0%)

end1.0ms (0.0%)