Random Numbers, Continued

• One time initialization of an array to a random set of values works as follows ...

  for i := 1 to n do
  
  for j := 1 to m do
  [i,j] A := llrand(seed);
  end;
  end;

• For random arrays, pick a larger separation

for i := 1 to n do

for j := 1 to m do

[i,j] A := llrand(seed);

for k := 1 to 9999 do -- spin generator

temp := llrand(seed); -- to separate

end; -- samples

end;

end;

r0 r10000 r20000 r30000

r40000 r50000 r60000 r70000

r0 r1 r2 r3

r4 r5 r r7

Previous slide

Next slide

Back to first slide

View graphic version