source: trunk/test/vector_test.cc @ 790

// $Id: vector_test.cc 790 2007-03-10 20:49:31Z jari $

/*
  Copyright (C) The authors contributing to this file.

  This file is part of the yat library, http://lev.thep.lu.se/trac/yat

  The yat library is free software; you can redistribute it and/or
  modify it under the terms of the GNU General Public License as
  published by the Free Software Foundation; either version 2 of the
  License, or (at your option) any later version.

  The yat library is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with this program; if not, write to the Free Software
  Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA
  02111-1307, USA.
*/

#include "yat/utility/Exception.h"
#include "yat/utility/FileUtil.h"
#include "yat/utility/utility.h"
#include "yat/utility/vector.h"

#include <fstream>
#include <sstream>

using namespace theplu::yat;

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

int main(const int argc,const char* argv[])
{  
  std::ostream* message;
  if (argc>1 && argv[1]==std::string("-v"))
    message = &std::cerr;
  else {
    message = new std::ofstream("/dev/null");
    if (argc>1)
      std::cout << argv[0] << " -v : for printing extra "  << "information\n";
  }
  *message << "testing vector" << std::endl;
  bool ok = true;

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

  // checking that shuffle works
  *message << "shuffle" << std::endl;
  double sum_before = utility::sum(vec);
  shuffle(vec);
  double sum_after = utility::sum(vec);
  ok &= (sum_after==sum_before);

  // checking that view works
  *message << "view" << std::endl;
  sum_before=0;
  for (unsigned int i=0; i<vec.size(); i+=2) 
    sum_before+=vec[i];
  utility::vector vec_view(vec,0,6,2);
  sum_after=utility::sum(vec_view);
  ok &= (sum_after==sum_before);
  vec[0]=0;
  vec_view[0]=24;
  ok &= (vec[0]==vec_view[0]);

  // Test of const view implementation (make sure no zero pointer is
  // used). Here we use the bad style of not making the view const!
  {
    *message << "const view implementation" << std::endl;
    const utility::vector vv(10,3.0);
    utility::vector vview(vv,0,5,1);
    // const utility::vector vview(vv,0,5,1); // this is the proper line
    utility::vector vv2(5,2.0);
    vv2.mul(vview); // should work even without const since const arg passing
    vv2.div(vview); // should work even without const since const arg passing
    ok &= (vview*vview == 3.0*3.0*vview.size());
  }

  // checking that copy constructor creates an independent object when
  // a non-view vector is copied
  {
    *message << "copy constructor" << std::endl;
    utility::vector vec2(vec);
    ok &= (vec.size()==vec2.size());
    ok &= (vec2==vec);
    ok &= (&vec2 != &vec);
    ok &= !vec2.isview();
  }

  // checking that copy constructor creates an independent object when
  // a view vector is copied
  {
    *message << "copy contructor on view" << std::endl;
    utility::vector vec3(vec_view);
    ok &= (vec_view.size()==vec3.size());
    ok &= (vec3 == vec_view);
    ok &= (&vec3 != &vec_view);
    ok &= !vec3.isview();
  }

  // checking that assignment operator throws an exception if vectors
  // differ in size
  {
    *message << "assignment operator" << std::endl;
    // GSL will catch the error in this test there for the GSL error
    // handling must be disabled until after the exception is
    // catched. The GSL error handler is reinstated after the
    // try-catch construct.
    gsl_error_handler_t* err_handler=gsl_set_error_handler_off();
    bool exception_happens=false;
    try {
      utility::vector v(vec_view.size()+1,0.0);
      v=vec_view;
    } catch (utility::GSL_error& err) {
      exception_happens=true;
    }
    if (!exception_happens) {
      *message << "Vector assignment operator did not throw expected exception"
               << std::endl;
      ok=false;
    }
    gsl_set_error_handler(err_handler);
  }

  // checking that assignment operator changes the underlying object when
  // a view is changed.
  {
    *message << "assignment operator on view" << std::endl;
    vec[3]=vec[4]=vec[5]=13;
    utility::vector vec_view(vec,3,3,1);
    utility::vector vec2(3,123.0);
    vec_view=vec2;
    if (vec[3]!=vec_view[0] || vec[4]!=vec_view[1] || vec[5]!=vec_view[2])
      ok=false;
  }

  // checking clone functionality
  {
    *message << "clone functionality" << std::endl;
    bool this_ok=true;
    *message << "\tcloning normal vector" << std::endl;
    utility::vector vec2(3,123.0);
    vec2.clone(vec);
    if (vec.size()!=vec2.size())
      this_ok=false;
    else
      for (unsigned int i=0; i<vec.size(); ++i)
        if (vec(i)!=vec2(i))
          this_ok=false;
    if (vec.gsl_vector_p()==vec2.gsl_vector_p())
      this_ok=false;
    *message << "\tcloning vector view" << std::endl;
    utility::vector* vec_view=new utility::vector(vec,3,3,1);
    utility::vector vec_view2;
    vec_view2.clone(*vec_view);
    if (!vec_view2.isview())
      this_ok=false;
    if (vec_view->size()!=vec_view2.size())
      this_ok=false;
    else
      for (u_int i=0; i<vec_view2.size(); ++i)
        if ((*vec_view)(i)!=vec_view2(i))
          this_ok=false;
    *message << "\tcloned vector view independence" << std::endl;
    delete vec_view;
    if (vec[3]!=vec_view2[0] || vec[4]!=vec_view2[1] || vec[5]!=vec_view2[2])
      this_ok=false;

    if (!this_ok) {
      *message << "FAIL: clone test" << std::endl;
      ok=false;
    }
  }

  // checking that reading vectors from properly formatted files works
  try {
    *message << "stream" << std::endl;
    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());
    utility::vector vec1(data_stream1);
    ok &= (vec1.size()==9);
    vec1=utility::vector(data_stream2);
    ok &= (vec1.size()==9);
    utility::vector vec2(data_stream3);
    ok &= (vec2.size()==12);
    vec2=utility::vector(data_stream4);
    ok &= (vec2.size()==12);
  } catch (utility::IO_error& err) {
    *message << err.what() << std::endl;
    ok=false;
  }

  // Checking that vector stream input operator can read whatever
  // vector stream output operator produces.
  {
    *message << "checking that output stream is valid as an input stream"
             << std::endl;
    std::stringstream s;
    s << vec;
    utility::vector vec2(s);
    ok &= (vec==vec2);
  }

  // Checking that badly formatted files are not accepted, or at least
  // throws an exception for the programmer to catch.
  *message << "checking that bad stream are rejected" << std::endl;
  bool this_ok=false;
  // Checking that unexpected characters in a stream gives error. 
  try {
    std::string data5("data/vector5.data");
    check_file_access(data5);
    std::ifstream data_stream5(data5.c_str());
    utility::vector dummy(data_stream5); // this will give an exception
  } catch (utility::IO_error& err) {
    *message << err.what() << std::endl;
    this_ok=true; // good, exceoption thrown, test passed
  }
  try {
    std::string data("data/vector6.data");
    check_file_access(data);
    std::ifstream data_stream(data.c_str());
    utility::vector dummy(data_stream); // this will give an exception
  } catch (utility::IO_error& err) {
    *message << err.what() << std::endl;
    this_ok=true; // good, exceoption thrown, test passed
  }
  ok &= this_ok;

  this_ok=false;
  try {
    std::string data("data/vector7.data");
    check_file_access(data);
    std::ifstream data_stream(data.c_str());
    vec=utility::vector(data_stream); // this will give an exception
  } catch (utility::IO_error& err) {
    *message << err.what() << std::endl;
    this_ok=true; // good, exceoption thrown, test passed
  }
  ok &= this_ok;

  if (!ok)
    *message << "vector test failed" << std::endl;

  return (ok ? 0 : -1);
}