newmat9.cpp

Go to the documentation of this file.

//$$ newmat9.cpp         Input and output

// Copyright (C) 1991,2,3,4: R B Davies


#ifndef WANT_STREAM
#define WANT_STREAM
#endif

#include "include.h"

#include "newmat.h"
#include "newmatio.h"
#include "newmatrc.h"
#include <iomanip>

using namespace std;

#ifdef use_namespace
namespace NEWMAT {
#endif



#ifdef DO_REPORT
#define REPORT { static ExeCounter ExeCount(__LINE__,9); ++ExeCount; }
#else
#define REPORT {}
#endif

// for G++ 3.01
#ifndef ios_format_flags
#define ios_format_flags long
#endif

ostream& operator<<(ostream& s, const BaseMatrix& X)
{
   GeneralMatrix* gm = (const_cast<BaseMatrix&>(X)).Evaluate(); operator<<(s, *gm);
   gm->tDelete(); return s;
}


ostream& operator<<(ostream& s, const GeneralMatrix& X)
{
   MatrixRow mr(const_cast<GeneralMatrix*>(&X), LoadOnEntry);
   int w = s.width();  int nr = X.Nrows();  ios_format_flags f = s.flags();
   s.setf(ios::fixed, ios::floatfield);
   for (int i=1; i<=nr; i++)
   {
      int skip = mr.skip;  int storage = mr.storage;
      Real* store = mr.data;  skip *= w+1;
      while (skip--) s << " ";
      while (storage--) { s.width(w); s << *store++ << " "; }
//      while (storage--) s << setw(w) << *store++ << " ";
      mr.Next();  s << "\n";
   }
   s << flush;  s.flags(f); return s;
}

ostream& operator<<(ostream& s, const printmat& p) {
  const int printwidth = s.width();
  s.width(0);
  s << "[";
  for (int i = 1; i <= p.mat.Nrows(); i++) {
    if (i > 1) s << "; ";
    for (int j = 1; j <= p.mat.Ncols(); j++) {
      if ( j > 1 ) s << ", ";
      s << setw(printwidth) << p.mat(i,j);
    }
  }
  s << "]";
  return s;
}


// include this stuff if you are using an old version of G++
// with an incomplete io library

/*

ostream& operator<<(ostream& os, Omanip_precision i)
   { os.precision(i.x); return os; }

Omanip_precision setprecision(int i) { return Omanip_precision(i); }

ostream& operator<<(ostream& os, Omanip_width i)
   { os.width(i.x); return os; }

Omanip_width setw(int i) { return Omanip_width(i); }

*/

#ifdef use_namespace
}
#endif