mcsstreammanager.cpp

Free C++ toolkit to facilitate Monte-Carlo simulation. This is a library covered under the LGPL. "MC

/*
    MCS-libre, a library to facilite Monte-Carlo simulation.
    Copyright (C) 2002  D. E. Williams

    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation; either
    version 2.1 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

    The full text of the GNU Lesser General Public License is
    contained in the files COPYING and LICENSE.

    See the 'MCS-libre-HOWTO' in the 'doc' directory for author's
    contact information.  

    Have a nice day!
*/

#include <iostream>
#include <fstream>

#include "mcsStreamManager.h"

streamManager::streamManager()
{
  streamlist.push_back( *(new stream(1636651392)));
  streamlist.push_back( *(new stream(202148464)));
  streamlist.push_back( *(new stream(1320382080)));
  streamlist.push_back( *(new stream(752575680)));
  streamlist.push_back( *(new stream(652689536)));
  streamlist.push_back( *(new stream(1792670848)));
  streamlist.push_back( *(new stream(2124451328)));
  streamlist.push_back( *(new stream(10470686)));
  streamlist.push_back( *(new stream(301214048)));
  streamlist.push_back( *(new stream(404617120)));
  streamlist.push_back( *(new stream(1988574336)));
  streamlist.push_back( *(new stream(715252800)));
  streamlist.push_back( *(new stream(894628032)));
  streamlist.push_back( *(new stream(1607214080)));
  streamlist.push_back( *(new stream(1149578368)));
  streamlist.push_back( *(new stream(1798325632)));
  streamlist.push_back( *(new stream(1424393856)));
  streamlist.push_back( *(new stream(1620518528)));
  streamlist.push_back( *(new stream(1197292672)));
  streamlist.push_back( *(new stream(1985659648)));
  streamlist.push_back( *(new stream(1433123712)));
  streamlist.push_back( *(new stream(425081408)));
  streamlist.push_back( *(new stream(1784672384)));
  streamlist.push_back( *(new stream(962164160)));
  streamlist.push_back( *(new stream(1389965056)));
  streamlist.push_back( *(new stream(615296064)));
  streamlist.push_back( *(new stream(183282288)));
  streamlist.push_back( *(new stream(1310163200)));
  streamlist.push_back( *(new stream(1329164544)));
  streamlist.push_back( *(new stream(748796736)));
  streamlist.push_back( *(new stream(1490840320)));
  streamlist.push_back( *(new stream(174655824)));
  streamlist.push_back( *(new stream(1921753984)));
  streamlist.push_back( *(new stream(90691488)));
  streamlist.push_back( *(new stream(1844555776)));
  streamlist.push_back( *(new stream(2118883328)));
  streamlist.push_back( *(new stream(599019776)));
  streamlist.push_back( *(new stream(622696448)));
  streamlist.push_back( *(new stream(2025160704)));
  streamlist.push_back( *(new stream(668767808)));
  streamlist.push_back( *(new stream(1803721344)));
  streamlist.push_back( *(new stream(205118112)));
  streamlist.push_back( *(new stream(407310336)));
  streamlist.push_back( *(new stream(570436992)));
  streamlist.push_back( *(new stream(996164352)));
  streamlist.push_back( *(new stream(715042432)));
  streamlist.push_back( *(new stream(702159168)));
  streamlist.push_back( *(new stream(560048256)));
  streamlist.push_back( *(new stream(1032647488)));
  streamlist.push_back( *(new stream(1075215488)));
  streamlist.push_back( *(new stream(2134198016)));
  streamlist.push_back( *(new stream(7679534)));
  streamlist.push_back( *(new stream(1373992320)));
  streamlist.push_back( *(new stream(717163840)));
  streamlist.push_back( *(new stream(1923909120)));
  streamlist.push_back( *(new stream(360446496)));
  streamlist.push_back( *(new stream(1354610816)));
  streamlist.push_back( *(new stream(410985568)));
  streamlist.push_back( *(new stream(670076928)));
  streamlist.push_back( *(new stream(620837376)));
  streamlist.push_back( *(new stream(1909764224)));
  streamlist.push_back( *(new stream(776868608)));
  streamlist.push_back( *(new stream(124742400)));
  streamlist.push_back( *(new stream(1990720640)));
  streamlist.push_back( *(new stream(1930928000)));
  streamlist.push_back( *(new stream(396497184)));
  streamlist.push_back( *(new stream(1278990848)));
  streamlist.push_back( *(new stream(1624754048)));
  streamlist.push_back( *(new stream(707575040)));
  streamlist.push_back( *(new stream(2024914048)));
  streamlist.push_back( *(new stream(903738624)));
  streamlist.push_back( *(new stream(764248192)));
  streamlist.push_back( *(new stream(884995072)));
  streamlist.push_back( *(new stream(20226652)));
  streamlist.push_back( *(new stream(956673408)));
  streamlist.push_back( *(new stream(1182072320)));
  streamlist.push_back( *(new stream(196679984)));
  streamlist.push_back( *(new stream(1716410112)));
  streamlist.push_back( *(new stream(1184083840)));
  streamlist.push_back( *(new stream(464176384)));
  streamlist.push_back( *(new stream(476144576)));
  streamlist.push_back( *(new stream(1611158400)));
  streamlist.push_back( *(new stream(539038208)));
  streamlist.push_back( *(new stream(2069039232)));
  streamlist.push_back( *(new stream(2045549056)));
  streamlist.push_back( *(new stream(144749280)));
  streamlist.push_back( *(new stream(910373056)));
  streamlist.push_back( *(new stream(1416116096)));
  streamlist.push_back( *(new stream(1371423872)));
  streamlist.push_back( *(new stream(115334704)));
  streamlist.push_back( *(new stream(874800640)));
  streamlist.push_back( *(new stream(31421040)));
  streamlist.push_back( *(new stream(1575471616)));
  streamlist.push_back( *(new stream(680241408)));
  streamlist.push_back( *(new stream(188537536)));
  streamlist.push_back( *(new stream(343140224)));
  streamlist.push_back( *(new stream(1275366272)));
  streamlist.push_back( *(new stream(1352005376)));
  streamlist.push_back( *(new stream(1629645952)));
  streamlist.push_back( *(new stream(1026296704)));
}

double streamManager::get_random_number_given_stream(int i)
{
  if ( (i < 0) or (i >= streamlist.size()))
    {
      cerr << "ERROR:  You have requested stream " << i << "." << endl;
      cerr << "mcsLibre (with the given initial conditions) only knows " << endl;
      cerr << "about streams 0 through " << streamlist.size() - 1 << "." << endl;
      exit(1);
    }
  else
    {
  return streamlist[i].get_new_random_value();
    }
}

int streamManager::get_number_of_streams()
{
  return streamlist.size();
}

void streamManager::loadNewSeeds(string s)
{
  streamlist.clear();
  ifstream input_file(s.data());
  double seed;
  while (! input_file.eof())
    {
      input_file >> seed;
      streamlist.push_back( *(new stream(seed)));      
    }
  input_file.close();
}