source: trunk/test/vector_test.cc @ 341

// $Id: vector_test.cc 341 2005-06-07 14:41:25Z jari $

#include <c++_tools/utility/Exception.h>
#include <c++_tools/utility/FileIO.h>
#include <c++_tools/gslapi/vector.h>
#include <c++_tools/gslapi/utility.h>

#include <fstream>
#include <sstream>

using namespace theplu;

void check_file_access(std::string& str)
{
  if (utility::FileIO().access_rights(str,"r")) {
    std::cerr << "test_nni: Cannot access file " << str << std::endl;
    exit(-1);
  }
}

int main(const int argc,const char* argv[])
{  
  bool print = (argc>1 && argv[1]==std::string("-p"));
  bool ok = true;

  gslapi::vector vec(12);
  for (unsigned int i=0; i<vec.size(); i++) 
    vec(i)=i;

  // checking that shuffle works
  double sum_before = vec.sum();
  shuffle(vec);
  double sum_after = vec.sum();
  if (sum_after != sum_before)
    ok = false;

  // checking that view works
  sum_before=0;
  for (unsigned int i=0; i<vec.size(); i+=2) 
    sum_before+=vec[i];
  gslapi::vector vec_view(vec,0,6,2);
  sum_after=vec_view.sum();
  if (sum_after != sum_before)
    ok = false;
  vec[0]=0;
  vec_view[0]=24;
  if (vec[0]!=vec_view[0])
    ok=false;

  // checking that copy contrutor creates an independent object
  gslapi::vector vec2(vec);
  if (vec.size()!=vec2.size())
    ok=false;
  if (&vec2 == &vec)
    ok=false;
  if (vec2.isview())
    ok=false;

  // checking that copy contrutor creates an independent object when
  // copying a view
  vec2=vec_view;
  if (vec_view.size()!=vec2.size())
    ok=false;
  if ((&vec2 == &vec_view) || (&vec2 == &vec))
    ok=false;
  if (vec2.isview())
    ok=false;
  vec2[0]=0;
  vec_view[0]=24;
  if (vec2[0]==vec_view[0])
    ok=false;

  // checking that reading vectors from properly formatted files works
  try {
    std::string data1("data/vector1.data");
    std::string data2("data/vector2.data");
    std::string data3("data/vector3.data");
    std::string data4("data/vector4.data");
    check_file_access(data1);
    check_file_access(data2);
    check_file_access(data3);
    check_file_access(data4);
    std::ifstream data_stream1(data1.c_str());
    std::ifstream data_stream2(data2.c_str());
    std::ifstream data_stream3(data3.c_str());
    std::ifstream data_stream4(data4.c_str());
    vec=gslapi::vector(data_stream1);
    if (vec.size()!=9)
      ok=false;
    vec=gslapi::vector(data_stream2);
    if (vec.size()!=9)
      ok=false;
    vec=gslapi::vector(data_stream3);
    if (vec.size()!=12)
      ok=false;
    vec=gslapi::vector(data_stream4);
    if (vec.size()!=12)
      ok=false;
  } catch (utility::IO_error& err) {
    if (print)
      std::cerr << err.what() << std::endl;
    ok=false;
  }

  // Checking that vector stream input operator can read whatever
  // vector stream output operator produces.
  std::stringstream s;
  s << vec;
  vec2=gslapi::vector(s);
  if (!(vec==vec2))
    ok=false;

  // Checking that unexpected characters in a stream are skipped. When
  // file read is reimplemented to resemble C/C++ style this should
  // check something else or maybe nothing at all.
  std::string data5("data/vector5.data");
  check_file_access(data5);
  std::ifstream data_stream5(data5.c_str());
  vec=gslapi::vector(data_stream5);
  if (vec.size()!=8)
    ok=false;

  // Checking that badly formatted files are not accepted, or at least
  // throws an exception for the programmer to catch.
  bool this_ok=false;
  try {
    std::string data("data/vector6.data");
    check_file_access(data);
    std::ifstream data_stream(data.c_str());
    vec=gslapi::vector(data_stream); // this will gove an exception
  } catch (utility::IO_error& err) {
    if (print)
      std::cerr << err.what() << std::endl;
    this_ok=true; // good, exceoption thrown, test passed
  }
  if (!this_ok)
    ok=false;
  this_ok=false;
  try {
    std::string data("data/vector7.data");
    check_file_access(data);
    std::ifstream data_stream(data.c_str());
    vec=gslapi::vector(data_stream); // this will gove an exception
  } catch (utility::IO_error& err) {
    if (print)
      std::cerr << err.what() << std::endl;
    this_ok=true; // good, exceoption thrown, test passed
  }
  if (!this_ok)
    ok=false;

  return (ok ? 0 : -1);
}