9a11cfaf-b87f-4d6b-9321-84972669fce9 56 Johannes Piel Johannes Piel 4 5 false false 0 decimal

&[DATE] &[TIME] - &[FILENAME]

N μ σ rnorm1 R N Random : μ σ 2 sqrt 2 R * 1 - erfinv vectorize * (* + 21 line :

N μ σ rnorm2 R N Random : i 1 N range R i el R i el μ σ ICDF.Normal : for R 31 line :

N μ σ rnorm3 N' N 20 / : R N' 1 matrix : i 1 N' range R i el 20 μ σ rnorm1 : for R R mat2sys sys2mat : 41 line :

σ 13 / :

μ 10 :

N 1000 :

N 5000 :

N 10000 :

N 100000 :

t.0 0 time :

r1a N μ σ rnorm1 eval :

0 time t.0 - 0.6 s *

0 time t.0 - 13.4 s *

0 time t.0 - 0.6 s *

t.0 0 time :

r2a N μ σ rnorm2 eval :

0 time t.0 - 0.2 s *

0 time t.0 - 0.8 s *

0 time t.0 - 1.6 s *

0 time t.0 - 70.3 s *

t.0 0 time :

r3a N μ σ rnorm3 eval :

0 time t.0 - 0.1 s *

0 time t.0 - 0.6 s *

0 time t.0 - 1.0 s *

0 time t.0 - 12.9 s *

I have no idea why this different definitions of rnorm execute with this different speeds. And even the fastest solution here is still slow...