SimParams.cc

/*

Copyright (C) 2003-2009 Kevin Thornton, krthornt[]@[]uci.edu

Remove the brackets to email me.

This file is part of libsequence.

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

libsequence 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 General Public License for more details.

You should have received a copy of the GNU General Public License
long with libsequence.  If not, see <http://www.gnu.org/licenses/>.

*/

#include <Sequence/SimParams.hpp>
#include <iostream>
#include <cctype>
#include <sstream>

namespace Sequence
{
  SimParams::SimParams (void):_command_line(""),_howmany(0),tsam(0)
  {}

  int SimParams::fromfile ( FILE * openfile )
  {
    _command_line.clear();
    _command_line.reserve(100);
    int ch;
    while( (ch=std::fgetc(openfile)) != EOF )
      {
        if (ch == '\n')
          break;
        else
          _command_line += char(ch);
      }
    //some versions of ms output the values with
    //which the RNG was seeded.  This deals with those cases
    while( (ch=std::fgetc(openfile)) != EOF )
      {
        if ((!isdigit(ch)&&!isspace(ch))|| ch=='\n')
          break;
      }
    std::istringstream in(_command_line.c_str());
    std::string ms;//for the program name...
    in >> ms >> tsam >> _howmany;
    return (ch);
  }

  std::istream & SimParams::read (std::istream & s)
  {
    _command_line.clear();
    char ch;
    //read in the command line
    while (s.get(ch))
      {
        if (ch == '\n')
          {
            break;
          }
        else
          {
            _command_line += ch;
          }
      }

    //some versions of ms output the values with
    //which the RNG was seeded.  This deals with those cases
    //    std::string seedline;
    while (1)
      {
        s.get (ch);
        if ((!isdigit(int(ch))&&!isspace(int(ch)))|| ch=='\n')
          break;
        //        else
        //          seedline += ch;
      }

    //open a stringstream to read the sample size and # runs
    //from the command-line args
    std::istringstream in(_command_line.c_str());
    std::string ms;//for the program name...
    in >> ms >> tsam >> _howmany;
    return s;
  }

  std::ostream & operator<< (std::ostream & stream, class SimParams & object)
  {
    stream << object._command_line;
    return (stream);
  }

  std::istream& operator>>(std::istream& s,  SimParams& c)
  {
    return c.read (s);
  }
}