Ticket #844: _reloaderFromMcmc.tol

/////////////////////////////////////////////////////////////////////////////
// LIBRARY: StdtLib (Standard Library of TOL)
// FILE: _reloaderFromMcmc.tol
// PURPOSE: 
/////////////////////////////////////////////////////////////////////////////

/////////////////////////////////////////////////////////////////////////////
Class @McmcReloaderFromMcmc : @McmcReloader
/////////////////////////////////////////////////////////////////////////////
{

// (1) por valor:
//  Matrix _.mcmc;
// (2) por referencia:
  Set _.mcmcSet = Copy(Empty);

  Set  _.colNames;

  ////////////////////////////////////////////////////////////////////////////
  Matrix _getMcmc(Real void)
  ////////////////////////////////////////////////////////////////////////////
  {
// (1) por valor:
//    _.mcmc;
// (2) por referencia:
    _.mcmcSet[1]
  };

  ////////////////////////////////////////////////////////////////////////////
  Set selectIndexes(Set selectedNames)
  ////////////////////////////////////////////////////////////////////////////
  {
    If(!Card(selectedNames) | !Card(_.colNames) , Copy(Empty),
    {
      Real If(!HasIndexByName(selectedNames), SetIndexByName(selectedNames));
      SelectIndexByName(selectedNames, _.colNames)
    })
  };

  ////////////////////////////////////////////////////////////////////////////
  VMatrix reloadSimulation(Real numSim, Set indexes)
  ////////////////////////////////////////////////////////////////////////////
  {
    Real numRow = getRowSequence(numSim);
    Matrix row = Sub(_getMcmc(?), numRow, 1, 1, Columns(_getMcmc(?)));
    Mat2VMat(SubCol(row, indexes), TRUE)
  };

  ////////////////////////////////////////////////////////////////////////////
  Static @McmcReloaderFromMcmc New(Matrix mcmc, Set colNames, Real sequential, 
    Real burnin)
  ////////////////////////////////////////////////////////////////////////////
  {
    Real mcmc.rows = Rows(mcmc);
    Real mcmc.cols = Columns(mcmc);
    Case(mcmc.cols!=Card(colNames), {
      WriteLn("[@McmcReloaderFromMcmc] MCMC columns differs from "
        <<"number of parameters","E");
      @McmcReloader::Null(?) 
    }, burnin>=mcmc.rows, {
      WriteLn("[@McmcReloaderFromMcmc] MCMC rows is less or equal "
        <<"than specified burnin","E");
      @McmcReloader::Null(?)
    }, True, {
      @McmcReloaderFromMcmc reloader = [[
        Set  _.colNames  = colNames; 
// (1) por valor:
//        Matrix _.mcmc = mcmc;
// (2) por referencia:
        Set _.mcmcSet = [[ mcmc ]];
        Real _.sampleLength = Rows(mcmc) - burnin;
        Real _.sequential = sequential;
        Real _.burnin = burnin
      ]];
      Real reloader::initialize.sequence(0);
      Real SetIndexByName(reloader::_.colNames);
      reloader
    })
  }

};