% GETSETRANDSTATE gets or sets the random state of the MATLAB functions 'rand' and 'randn'
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DESCRIPTION: GETSETRANDSTATE gets or sets the random state of the MATLAB
%              functions 'rand' and 'randn'. This function looks in the
%              input plist for the key-word 'rand_state' (LTPDA toolbox
%              less than 2.1) or 'rand_stream' (LTPDA toolbox greate equal
%              than 2.1) to set the random state. If these key words
%              doesn't exist in the plist then stores this function the
%              random state in the plist.
%
% CALL:        pl = getSetRandState(pl_in);
%              pl = pl_in.getSetRandState();
%              pl_in.getSetRandState();
%
% <a href="matlab:utils.helper.displayMethodInfo('plist', 'getSetRandState')">Parameters Description</a>
%
% VERSION:     $Id: getSetRandState.m,v 1.10 2011/08/15 10:42:00 hewitson Exp $
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function varargout = getSetRandState(varargin)
  
  %%% Check if this is a call for parameters
  if utils.helper.isinfocall(varargin{:})
    varargout{1} = getInfo(varargin{3});
    return
  end
  
  % Collect all PLISTs
  plin = utils.helper.collect_objects(varargin(:), 'plist');
  
  % Decide on a deep copy or a modify
  plout = copy(plin, nargout);

  if numel(plout) ~= 1
    error('### This method accepts only one input plist object.');
  end
  
  def_struct = struct('Type',          '',    ...
                      'NumStreams',    [],    ...
                      'StreamIndex',   [],    ...
                      'Substream',      1,    ...
                      'Seed',           0,    ...
                      'State',         [],    ...
                      'RandnAlg',      '',    ...
                      'Antithetic',    false, ...
                      'FullPrecision', true);
  
  rand_state  = plout.find('rand_state');
  rand_stream = plout.find('rand_stream');
  
  %%% We have different behavior for the MATLAB version 2008b and less than
  %%% version 2008b because version 2008b is working with random streams.
  
  %%% MATLAB version 2008a and less
  if verLessThan('MATLAB', '7.7')
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %%%                    MATLAB version 2008a and less                    %%%
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    if ~isempty(rand_state)
      %%%    Reset random state    %%%

      %%% It might be that 'rand_state' contains the 'state' or the 'seed'
      %%% of the MATLAB rand or randn methods.
      try
        randn('seed',rand_state);
        rand('seed', rand_state);
      catch
        if numel(rand_state) > 2
          rand('state', rand_state)
        else
          randn('state', rand_state)
        end
      end
      plout.remove('rand_state');
      plout.getSetRandState();
      
    elseif ~isempty(rand_stream)
      %%%    Reset random state    %%%
            
      %%% Legacy mode
      if numel(rand_stream) ~= 1
        %%% Set that of 'randn'
        randn('state', rand_stream(1).State)
        %%% Set that of 'rand'
        rand('state', rand_stream(2).State)
      else
        error('### I have no idea what I should do because the random state is stores as a random stream with MATLAB version 2008b or later.')
      end
      
    else
      %%%    Store random state    %%%

      stream = [def_struct def_struct];
      %%% Store state of 'randn'
      stream(1).Type  = 'randn';
      stream(1).State = randn('state');
      %%% Store state of 'rand'
      stream(2).Type  = 'rand';
      stream(2).State = rand('state');
      plout.append('rand_stream', stream);
    end
    
  else
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %%%                   MATLAB version 2008b and later                    %%%
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    def_stream  = RandStream.getDefaultStream;

    if ~isempty(rand_state)
      %%%    Reset random state    %%%
      
      %%% It might be that 'rand_state' contains the 'state' or the 'seed'
      %%% of the MATLAB rand or randn methods.
      try
        randn('seed',rand_state);
        rand('seed', rand_state);
      catch
        if numel(rand_state) > 2
          rand('state', rand_state)
        else
          randn('state', rand_state)
        end
      end
      plout.remove('rand_state');
      plout.getSetRandState();

    elseif ~isempty(rand_stream)
      %%%    Reset random state    %%%
      
      %%% Legacy mode
      if numel(rand_stream) ~= 1
        
        %%% Set that of 'randn'
        randn('state', rand_stream(1).State)
        %%% Set that of 'rand'
        rand('state', rand_stream(2).State)
      else
        stream = RandStream(rand_stream.Type, 'Seed', rand_stream.Seed, 'RandnAlg', rand_stream.RandnAlg);
        stream.State         = uint32(rand_stream.State);
        stream.Substream     = rand_stream.Substream;
        stream.Antithetic    = rand_stream.Antithetic;
        stream.FullPrecision = rand_stream.FullPrecision;
        RandStream.setDefaultStream(stream);
      end
      
    else
      %%%    Store random state    %%%

      if strcmp(def_stream.Type, 'legacy')
        %%%%%%%%%%   Legacy Mode   %%%%%%%%%%

        stream = [def_struct def_struct];
        %%% Store state of 'randn'
        stream(1).Type  = 'randn';
        stream(1).State = randn('state');
        %%% Store state of 'rand'
        stream(2).Type  = 'rand';
        stream(2).State = rand('state');

      else
        %%%%%%%%%%   Stream Mode   %%%%%%%%%%
        stream = def_struct;
        stream.Type          = def_stream.Type;
        stream.NumStreams    = def_stream.NumStreams;
        stream.StreamIndex   = def_stream.StreamIndex;
        stream.Substream     = def_stream.Substream;
        stream.Seed          = def_stream.Seed;
        stream.State         = double(def_stream.State);
        stream.RandnAlg      = def_stream.RandnAlg;
        stream.Antithetic    = def_stream.Antithetic;
        stream.FullPrecision = def_stream.FullPrecision;
      end
      plout.pset('rand_stream', stream);
    end
    
    
  end
  
  varargout{1} = plout;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% FUNCTION:    getInfo
%
% DESCRIPTION: Get Info Object
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function ii = getInfo(varargin)
  if nargin == 1 && strcmpi(varargin{1}, 'None')
    sets = {};
    pl   = [];
  else
    sets = {'Default'};
    pl   = getDefaultPlist;
  end
  % Build info object
  ii = minfo(mfilename, 'plist', 'ltpda', utils.const.categories.internal, '$Id: getSetRandState.m,v 1.10 2011/08/15 10:42:00 hewitson Exp $', sets, pl);
  ii.setArgsmin(1);
  ii.setArgsmax(1);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% FUNCTION:    getDefaultPlist
%
% DESCRIPTION: Get Default Plist
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function plo = getDefaultPlist()
  plo = plist();
end