Re: sampling from two instances of a random number generator

[Søren Hauberg]

søn, 25 05 2008 kl. 22:18 +0200, skrev Søren Hauberg:
>   function retval = myrand (sz, state)
>     rand ("state", state);
>     retval = rand (sz);
>   endfunction

Aarrghh, posted the wrong function. It should of course only set the
state once. So, something like this should do:

  function retval = myrand (sz, state),
    persistent static_state = [];
    if (isempty (static_state))
      static_state = state;
      rand ("state", state);
    endif
    retval = rand (sz);
  endfunction

Søren

P.S. This code isn't really tested, so you should check it before using
it :-)