dynare/dynare++/sylv/testing/MMMatrix.cc

/* $Header: /var/lib/cvs/dynare_cpp/sylv/testing/MMMatrix.cpp,v 1.1.1.1 2004/06/04 13:01:13 kamenik Exp $ */

/* Tag $Name:  $ */

#include "MMMatrix.hh"

#include <fstream>
#include <iomanip>

MMMatrixIn::MMMatrixIn(const std::string &fname)
{
  std::ifstream fd{fname};
  if (fd.fail())
    throw MMException("Cannot open file "+fname+" for reading\n");

  // jump over initial comments
  while (fd.peek() == '%')
    fd.ignore(std::numeric_limits<std::streamsize>::max(), '\n');

  // read in number of rows and cols
  fd >> rows >> cols;
  if (fd.fail())
    throw MMException("Couldn't parse rows and cols\n");
  // read in data
  data = std::shared_ptr<double>(static_cast<double *>(operator new[](rows*cols*sizeof(double))), [](double *arr) { operator delete[](static_cast<void *>(arr)); });
  int len = rows*cols;
  int i = 0;
  while (!fd.eof() && i < len)
    {
      fd >> data.get()[i];
      if (fd.fail())
        throw MMException("Couldn't parse float number\n");
      i++;
    }
  if (i < len)
    throw MMException("Couldn't read all " + std::to_string(len) + " elements, read "
                      + std::to_string(i) + " so far\n");
  fd.close();
}

void
MMMatrixOut::write(const std::string &fname, const GeneralMatrix &m)
{
  std::ofstream fd{fname, std::ios::out | std::ios::trunc};
  if (fd.fail())
    throw MMException("Cannot open file "+fname+" for writing\n");

  fd << "%%%%MatrixMarket matrix array real general" << std::endl
     << m.numRows() << ' ' << m.numCols() << std::endl
     << std::setprecision(35);
  for (int i = 0; i < m.numCols(); i++)
    for (int j = 0; j < m.numRows(); j++)
      fd << std::setw(40) << m.get(i, j);
  fd.close();
}