Details

rewrite4.4m (50.4%)

Algorithm
71×rewrite-expression-head
Calls

181 calls:

1.9m
(+ (/ -176.6150291621406 (+ 4.0 (- z))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (+ (/ 771.3234287776531 (+ (- z) 3.0)) (+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z))))) (/ -1259.1392167224028 (+ 2.0 (- z)))) (/ 12.507343278686905 (+ 5.0 (- z)))) (+ (/ 1.5056327351493116e-07 (+ 8.0 (- z))) (/ 9.984369578019572e-06 (+ 7.0 (- z))))))
1.1m
(* (/ (* (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))) (+ (+ (/ -176.6150291621406 (+ (- z 1.0) 4.0)) (+ (+ (/ 676.5203681218851 z) 0.9999999999998099) (/ -1259.1392167224028 (+ (- z 1.0) 2.0)))) (+ (/ 771.3234287776531 (+ (- z 1.0) 3.0)) (+ (+ (/ 12.507343278686905 (+ (- z 1.0) 5.0)) (/ -0.13857109526572012 (+ (- z 1.0) 6.0))) (+ (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0)) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)))))))
24.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)
14.2s
(* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))
13.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)))))
Rules
2017×*-un-lft-identity
1680×times-frac
1634×add-sqr-sqrt
1236×add-cube-cbrt
1011×add-exp-log
851×add-cbrt-cube
810×frac-add
631×pow1
484×associate-*l*
436×distribute-lft-out
427×flip3-+ flip-+
418×prod-exp
390×associate-*r*
381×add-log-exp
324×cbrt-unprod
266×sqrt-prod
216×associate-*r/
206×unpow-prod-down
191×frac-times
188×pow-prod-down
170×div-exp
168×associate-/l*
148×associate-/r/
129×cbrt-undiv
122×cbrt-prod
121×unswap-sqr
118×sum-log
114×cube-prod
108×div-inv
92×associate-*l/
91×associate-/r*
79×log-pow
69×sqrt-div
65×difference-of-squares
64×*-commutative
58×cbrt-div pow-exp
52×flip-- flip3--
45×pow-unpow
44×sqr-pow
42×swap-sqr
40×distribute-lft-out-- associate-/l/
36×distribute-rgt-in distribute-lft-in
35×log-prod
34×pow-to-exp sqrt-pow1
30×pow1/2
28×pow1/3
26×+-commutative
25×frac-2neg clear-num
22×rem-cube-cbrt
20×cube-div rem-sqrt-square sub-neg
18×diff-log
16×pow-pow pow-prod-up
15×rec-exp 1-exp
14×cube-mult unpow3
13×pow-plus
12×pow-sqr
associate-+r+ pow2
associate-+l- div-sub
associate-+l+
associate--l+
associate--r+
un-div-inv exp-sum
exp-diff log-div pow-sub
cos-sum difference-cubes associate--r- unpow-prod-up pow3 rem-cbrt-cube sin-sum rem-log-exp associate-+r-

simplify2.9m (33.1%)

series36.1s (6.9%)

Calls

78 calls:

6.5s
(- (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 (pow (fabs x) 2))) 3)) (sqrt (* (/ (pow (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 3) (exp (* 2 (pow (fabs x) 2)))) (/ (pow (+ 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)))))))) 3) (exp (pow (fabs x) 2)))))))
4.5s
(- (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 (pow (fabs x) 2))) 3))
4.5s
(- (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 (pow (fabs x) 2))) 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 (pow (fabs x) 2))) 3))))
2.0s
(* (/ (+ 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 (pow (fabs x) 2))) (+ 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 (pow (fabs x) 2)))))
1.9s
(* (/ (+ 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 (pow (fabs x) 2))) (+ 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 (pow (fabs x) 2)))))

sample29.9s (5.7%)

Algorithm
64×intervals
Results
16.1s164147×body80valid
3.0s532×body10240exit
1.9s4172×body160valid
1.6s26271×body80nan
1.0s1070×body1280valid
469.0ms705×body640valid
229.0ms121×body1280nan
165.0ms92×body640nan
132.0ms364×body320valid
45.0ms18×body5120valid
38.0ms43×body320nan
31.0ms1536×pre80true
29.0ms46×body2560valid
19.0ms20×body160nan

prune17.0s (3.2%)

localize2.3s (0.4%)

regimes904.0ms (0.2%)

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

bsearch0.0ms (0.0%)

end0.0ms (0.0%)