source: trunk/test/vector_test.cc @ 328

// $Id: vector_test.cc 328 2005-05-30 23:04:20Z jari $

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

#include <fstream>

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 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 file works
  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;

  std::cout << "ok: " << ok << std::endl;

  if (ok) 
    return 0;
  return -1;
}