Details

simplify8.1m (39.2%)

rewrite4.9m (23.9%)

Algorithm
152×rewrite-expression-head
Calls

460 calls:

1.5m
(* (/ (* (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))))))
1.0m
(* (exp (- (- 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))))))
43.3s
(* (/ (* (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)) (sqrt (* PI 2.0))) (exp z)) (* (exp (- (- 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.6s
(* (/ (+ (+ (+ (+ (+ 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)
23.1s
(+ (/ -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))))))))
Rules
11422×times-frac
10413×*-un-lft-identity
8521×add-sqr-sqrt
6523×add-cube-cbrt
3601×sqrt-prod
2803×frac-add
2425×add-exp-log
1967×add-cbrt-cube
1535×distribute-lft-out
1487×pow1
1112×flip-+ flip3-+
1092×associate-*l*
888×associate-*r*
887×prod-exp
829×associate-/l*
819×prod-diff
708×associate-*r/
697×add-log-exp
673×frac-times
666×cbrt-unprod
591×difference-of-squares
573×distribute-lft-out--
532×associate-/r*
507×div-exp
463×log1p-expm1-u
461×expm1-log1p-u
458×unpow-prod-down
414×div-inv
406×pow-prod-down
383×cbrt-undiv
361×cbrt-prod
318×unswap-sqr cube-prod
267×associate-/r/
257×exp-diff
221×log-pow
218×cbrt-div
209×associate-*l/
202×sqrt-div
167×fma-def
147×distribute-rgt-in distribute-lft-in
145×fma-neg
142×pow1/3
126×sqr-pow
121×*-commutative
120×pow-exp
119×sum-log
117×flip-- flip3--
111×exp-sum
107×pow-prod-up
101×frac-2neg clear-num
88×associate-/l/ pow-to-exp
85×sqrt-pow1
81×pow-unpow
73×log-prod
68×associate--l+
66×fma-udef
63×swap-sqr
58×pow-plus
57×pow1/2 pow-sqr
52×sub-neg
49×1-exp rec-exp pow-pow
48×unpow3 cube-mult
39×diff-log
32×pow2
26×associate-+l+
24×+-commutative
23×cube-div
22×rem-cube-cbrt
20×frac-sub distribute-rgt-neg-in
19×rem-sqrt-square unpow2
17×div-sub
14×associate-+l-
13×pow-flip
11×inv-pow hypot-udef associate--r+
10×exp-prod
un-div-inv pow-sub
sub-div cos-sum
pow3 unpow-prod-up hypot-def rem-log-exp
sqrt-pow2 log1p-udef associate--l- log1p-expm1 expm1-log1p expm1-udef log-div
sqrt-undiv sin-sum associate-+r+
sqrt-unprod rem-square-sqrt rem-cbrt-cube
rem-exp-log distribute-rgt-out-- distribute-rgt-out associate--r- sum-cubes

sample3.9m (18.9%)

Algorithm
99×intervals
Results
1.3m15544×body10240exit
49.2s437767×body80valid
29.5s30263×body1280valid
17.7s23358×body640valid
7.9s69959×body80nan
6.6s3043×body1280nan
4.5s9993×body160valid
4.4s12027×body320valid
3.9s2442×body640nan
2.1s49536×pre80true
1.2s1468×body2560valid
1.2s1320×body320nan
355.0ms630×body160nan
354.0ms569×body5120valid

prune1.9m (9.4%)

series1.5m (7.4%)

Calls

460 calls:

1.9s
(fma -1259.1392167224028 (* (fma (/ 12.507343278686905 (- 5.0 z)) (/ 12.507343278686905 (- 5.0 z)) (* (+ (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))) (/ -0.13857109526572012 (- 6.0 z))) (- (+ (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))) (/ -0.13857109526572012 (- 6.0 z))) (/ 12.507343278686905 (- 5.0 z))))) (* (- 3.0 z) (- (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (/ -176.6150291621406 (- 4.0 z))))) (* (fma (fma (- (* (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))) (* (/ -176.6150291621406 (- 4.0 z)) (/ -176.6150291621406 (- 4.0 z)))) (- 3.0 z) (* (- (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (/ -176.6150291621406 (- 4.0 z))) 771.3234287776531)) (fma (/ 12.507343278686905 (- 5.0 z)) (/ 12.507343278686905 (- 5.0 z)) (* (+ (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))) (/ -0.13857109526572012 (- 6.0 z))) (- (+ (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))) (/ -0.13857109526572012 (- 6.0 z))) (/ 12.507343278686905 (- 5.0 z))))) (* (* (- 3.0 z) (- (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))) (/ -176.6150291621406 (- 4.0 z)))) (+ (pow (+ (+ (/ 1.5056327351493116e-07 (- 8.0 z)) (/ 9.984369578019572e-06 (+ (- z) 7.0))) (/ -0.13857109526572012 (- 6.0 z))) 3) (pow (/ 12.507343278686905 (- 5.0 z)) 3)))) (+ (- z) 2.0)))
1.7s
(* (pow (* (pow (log u1) 1.0) (pow -2.0 1.0)) 0.5) (- 0.16666666666666666 (* 0.3333333333333333 (* (pow u2 2) (pow PI 2)))))
1.5s
(cbrt (- (fma (exp (- (pow (fabs x) 2))) (+ (/ 0.284496736 (pow (fma (fabs x) 0.3275911 1.0) 2)) (/ 1.453152027 (pow (fma (fabs x) 0.3275911 1.0) 4))) 1.0) (fma (exp (- (pow (fabs x) 2))) (+ (/ 1.061405429 (pow (fma (fabs x) 0.3275911 1.0) 5)) (/ 0.254829592 (fma 0.3275911 (fabs x) 1.0))) (/ 1.421413741 (* (exp (pow (fabs x) 2)) (pow (fma (fabs x) 0.3275911 1.0) 3))))))
1.5s
(* (* (* (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5)) (sqrt (* PI 2.0))) (exp (- (- 1.0 7.0) 0.5))) (fma (+ (pow (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0))) 3) (pow (/ 12.507343278686905 (+ (- z 1.0) 5.0)) 3)) (* (* (+ (- z 1.0) 4.0) (* (+ (- z 1.0) 6.0) (+ (- 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))))))) (* (fma (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0))) (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0))) (* (/ 12.507343278686905 (+ (- z 1.0) 5.0)) (- (/ 12.507343278686905 (+ (- z 1.0) 5.0)) (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)))))) (fma -176.6150291621406 (* (* (+ (- z 1.0) 6.0) (+ (- 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 (fma -0.13857109526572012 (+ (- z 1.0) 3.0) (* (+ (- z 1.0) 6.0) 771.3234287776531)) (* 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) 6.0) (+ (- 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)))))))))))
1.4s
(pow (* (pow (log u1) 1.0) (pow -2.0 1.0)) 0.5)

regimes7.5s (0.6%)

Accuracy

Total 26.6b remaining (19.8%)

Threshold costs 7.5b (5.6%)

7.0b35.7%_divideComplex, real part
3.7b72.8%Octave 3.8, jcobi/1
2.8b0%Octave 3.8, jcobi/3
2.6b33.4%Octave 3.8, jcobi/4
1.9b0%math.sqrt on complex, real part

localize5.4s (0.4%)

bsearch1.3s (0.1%)

end0.0ms (0.0%)