source: trunk/test/matrix_test.cc @ 393

// $Id: matrix_test.cc 393 2005-09-14 18:25:03Z peter $

#include <c++_tools/gslapi/matrix.h>

#include <fstream>
#include <iostream>

int main(const int argc,const char* argv[])
{  
  using namespace theplu;
  std::ostream* error;
  if (argc>1 && argv[1]==std::string("-p"))
    error = &std::cerr;
  else 
    error = new std::ofstream("/dev/null");

  *error << "testing matrix" << std::endl;

  bool ok = true;

  gslapi::matrix m(3,3,9);

  gslapi::matrix m2(m);
  if (m2!=m)
    ok=false;
  std::ofstream my_out("data/tmp_test_matrix.txt");
  my_out << m2;
  my_out.close();

  std::ifstream is("data/tmp_test_matrix.txt");
  gslapi::matrix m3(is);
  is.close();
  if (m3!=m2)
    ok=false;

  gslapi::matrix m4(3,3,1);
  m4 *= 9;
  
  if (m4!=m){
    ok=false;
    *error << "error operator*=(double)" << std::endl;
  }


  if (error!=&std::cerr)
    delete error;

  return (ok ? 0 : -1);

}