Details

rewrite12.5m (39.4%)

Algorithm
151×rewrite-expression-head
Calls

442 calls:

3.3m
(* (/ (pow (sqrt (+ (- z -6) 0.5)) (+ (- z 1) 0.5)) (exp z)) (* (/ (pow (sqrt (+ (- z -6) 0.5)) (+ (- z 1) 0.5)) (exp (- -6))) (* (/ (sqrt (* PI 2)) (exp 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)))))))
2.3m
(* (/ (* (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.6m
(* (/ (sqrt (* PI 2)) (exp 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)) (exp (- z -6))) (* (/ (sqrt (* PI 2)) (exp 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))))))
17.2s
(* (* (/ 1 (sqrt PI)) (exp (* (fabs x) (fabs x)))) (+ (+ (+ (/ 1 (fabs x)) (* (/ 1 2) (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))))) (* (/ 3 4) (* (* (* (* (/ 1 (fabs x)) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))))) (* (/ 15 8) (* (* (* (* (* (pow (fabs x) -2) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))) (/ 1 (fabs x))))))
Rules
2907757×frac-times
2178468×associate-*l/
1608045×associate-*r/
1559390×frac-add
577526×un-div-inv
24089×flip-+ flip3-+
18905×times-frac
13003×*-un-lft-identity
11120×add-sqr-sqrt
9315×add-cube-cbrt
5972×sqrt-prod
3419×add-exp-log
2801×pow1
2305×add-cbrt-cube
1647×associate-*l*
1613×distribute-lft-out
1442×prod-exp
1171×associate-*r*
1035×prod-diff
922×cbrt-unprod
852×div-exp
806×pow-prod-down
599×associate-/l*
587×add-log-exp
584×div-inv
552×pow-prod-up
551×associate-/r/
513×associate-/r*
439×expm1-log1p-u insert-posit16 log1p-expm1-u
412×cbrt-undiv
399×exp-sum
365×difference-of-squares
354×unswap-sqr
310×log-pow
260×distribute-lft-out--
250×swap-sqr
249×pow-plus
220×sqrt-pow1
204×flip3-- flip--
189×1-exp cbrt-prod
185×rec-exp
175×unpow-prod-down
167×distribute-rgt-in distribute-lft-in
164×fma-neg
160×sqrt-div
154×pow-sqr
150×*-commutative
132×cbrt-div
127×log-prod
115×pow1/2
94×pow2
93×sub-neg
86×pow1/3
84×frac-2neg clear-num
81×pow-flip
78×inv-pow
75×associate-/l/
67×fma-udef
64×fma-def
51×diff-log
50×associate--l+ sqr-pow
45×pow-exp
40×rem-sqrt-square
30×unpow-prod-up
29×exp-prod
28×sum-log pow-to-exp
26×cube-unmult
24×associate-+l+
22×frac-sub
20×div-sub
18×distribute-rgt-out
16×+-commutative
15×hypot-def
14×distribute-rgt1-in associate--r+
12×exp-diff
11×rem-log-exp
pow3 pow-unpow log-div
associate-+l-
sqrt-unprod
associate--l- pow-pow
distribute-rgt-out-- rem-square-sqrt cos-sum
sin-sum difference-of-sqr-1 exp-to-pow rem-exp-log expm1-def
distribute-lft-neg-in pow-sub distribute-rgt-neg-in count-2
rem-cbrt-cube
expm1-log1p expm1-udef distribute-lft-neg-out distribute-lft1-in neg-mul-1 associate--r-

simplify11.3m (35.6%)

sample4.2m (13.3%)

Algorithm
98×intervals
Results
1.3m11831×body10240exit
1.2m437275×body80valid
34.7s26913×body1280valid
18.3s20824×body640valid
10.1s68853×body80nan
7.3s1536×body1280nan
6.2s10368×body320valid
3.6s9624×body160valid
3.3s1282×body640nan
2.5s49536×pre80true
1.7s1480×body2560valid
1.1s613×body320nan
543.0ms610×body5120valid
512.0ms359×body160nan

prune2.1m (6.5%)

series1.4m (4.4%)

Calls

439 calls:

1.2s
(* (cbrt (* (* 0.5 (cos re)) (- (* (* (* im im) im) -1/3) (fma (pow im 5) 1/60 (+ im im))))) (cbrt (* (* 0.5 (cos re)) (- (* (* (* im im) im) -1/3) (fma (pow im 5) 1/60 (+ im im))))))
1.1s
(* (* (pow (* -2 (log u1)) 0.5) (sqrt 1/6)) (sqrt 1/6))
1.1s
(cbrt (fma 1/60 (pow im 5) (+ im im)))
1.1s
(exp (* (pow (* -2 (log u1)) (/ 0.5 2)) (* (pow (* -2 (log u1)) (/ 0.5 2)) 1/6)))
1.1s
(* (pow (* -2 (log u1)) 0.5) 1/6)

regimes7.9s (0.4%)

Accuracy

Total 17.8b remaining (10.9%)

Threshold costs 4.0b (2.5%)

3.0b76.3%Octave 3.8, jcobi/1
2.6b0%_divideComplex, real part
2.3b0%_divideComplex, imaginary part
1.6b43.7%Octave 3.8, jcobi/4
1.6b0%Jmat.Real.lambertw, newton loop step

localize6.4s (0.3%)

bsearch1.5s (0.1%)

end1.0ms (0.0%)