Details

simplify3.3m (42.8%)

rewrite3.1m (39.8%)

Algorithm
74×rewrite-expression-head
Calls

182 calls:

1.4m
(* (/ (* (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.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)
24.2s
(+ (/ -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))))))))
14.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)))))
14.2s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
Rules
3008×*-un-lft-identity
2499×add-sqr-sqrt
2272×times-frac
2040×add-exp-log
1960×add-cube-cbrt
1943×add-cbrt-cube
963×prod-exp
914×cbrt-unprod
837×frac-add
811×prod-diff
761×pow1
631×associate-*l*
515×div-exp
505×distribute-lft-out
485×cbrt-undiv
394×add-log-exp
390×associate-*r*
347×flip-+ flip3-+
333×frac-times
241×associate-*l/
234×associate-*r/
214×pow-prod-down
187×associate-/l*
182×log1p-expm1-u expm1-log1p-u
175×cbrt-prod
174×sqrt-prod
159×distribute-rgt-in distribute-lft-in
147×unpow-prod-down
140×div-inv
129×exp-sum
128×difference-of-squares
121×fma-def
118×associate-/r/
116×distribute-lft-out--
115×swap-sqr
113×unswap-sqr
112×associate-/r*
109×sum-log
97×cube-prod
93×associate--l+
92×fma-neg
79×pow-prod-up
77×pow-plus
73×log-pow
70×flip3-- flip--
63×*-commutative
61×sub-neg
53×sqr-pow
40×cbrt-div
37×diff-log
35×sqrt-div
33×frac-2neg clear-num
32×associate-/l/
25×unpow3 cube-mult log-prod
23×pow-sqr
21×pow2 pow1/3
18×+-commutative sqrt-pow1
17×pow-exp
15×fma-udef
14×pow1/2
13×pow-to-exp rec-exp 1-exp
12×distribute-rgt-neg-in
pow-unpow
rem-sqrt-square associate--l-
cube-div hypot-udef associate-+r+ div-sub associate-+l-
pow-flip associate--r+
pow3
inv-pow pow-pow
exp-diff rem-cube-cbrt associate-+l+ exp-prod
pow-sub rem-log-exp
cos-sum difference-cubes sin-sum unpow-prod-up associate--r- rem-cbrt-cube log-div

sample31.9s (6.9%)

Algorithm
64×intervals
Results
16.4s160398×body80valid
3.5s532×body10240exit
1.6s26271×body80nan
1.5s1070×body1280valid
1.3s4172×body160valid
476.0ms705×body640valid
260.0ms121×body1280nan
185.0ms364×body320valid
113.0ms92×body640nan
40.0ms1536×pre80true
36.0ms43×body320nan
20.0ms46×body2560valid
12.0ms20×body160nan
7.0ms18×body5120valid

prune28.5s (6.2%)

series16.5s (3.6%)

Calls

93 calls:

1.4s
(- (pow 1.0 3) (pow (* (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (* (fma (fma (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (fma (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (fma 1.061405429 (* (/ (/ 1.0 (- 1.0 (* 0.3275911 (fabs x)))) (fma 0.3275911 (fabs x) 1.0)) (- 1.0 (* 0.3275911 (fabs x)))) -1.453152027) 1.421413741) -0.284496736) (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) 0.254829592) (exp (- (pow (fabs x) 2))))) 3))
1.3s
(- 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))))))
1.3s
(- 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 (* (fma 0.3275911 (fabs x) 1.0) (- 1.0 (* 0.3275911 (fabs x))))) (- 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))))
1.1s
(* (/ 1.0 6.0) (pow (* -2.0 (log u1)) 0.5))
679.0ms
(log (exp (- 1.0 (* (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (* (fma (fma (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (fma (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) (fma 1.061405429 (* (/ (/ 1.0 (- 1.0 (* 0.3275911 (fabs x)))) (fma 0.3275911 (fabs x) 1.0)) (- 1.0 (* 0.3275911 (fabs x)))) -1.453152027) 1.421413741) -0.284496736) (/ 1.0 (fma 0.3275911 (fabs x) 1.0)) 0.254829592) (exp (- (pow (fabs x) 2))))))))

localize2.9s (0.6%)

regimes568.0ms (0.1%)

Accuracy

Total 0.3b remaining (2%)

Threshold costs 0.3b (2%)

0.1b0%Jmat.Real.erfi, branch x less than or equal to 0.5
0.1b0%math.log10 on complex, imaginary part
0.0b0%Jmat.Real.erf
0.0b0%FastMath dist
0.0b0%_multiplyComplex, imaginary part

bsearch0.0ms (0.0%)

end0.0ms (0.0%)