#ifndef TEST_UBLAS_H
#define TEST_UBLAS_H
/***************************************************************************
    test_ublas.h : Test uBLAS
                             -------------------
    begin                : Sun 2005-03-14
    copyright            : (C) 2005 by Paul C. Leopardi
    email                : leopardi@bigpond.net.au
 ***************************************************************************
 *   This library is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU Lesser General Public License as        *
 *   published by the Free Software Foundation; either version 2.1 of the  *
 *   License, or (at your option) any later version.                       *
 *   See http://www.fsf.org/copyleft/lesser.html for details               *
 ***************************************************************************/

// Declarations
#define HAVE_STD_ABS
#include <boost/limits.hpp>

#include <exception>
#include <stdexcept>

// Include <utility> to define pair<>
#include <utility>
#include <vector>
#include <fstream>
#include <string>
#include <complex>

// Add real equivalents to complex functions
#include <math.h>

#include <boost/numeric/ublas/config.hpp>
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/triangular.hpp>
#include <boost/numeric/ublas/vector_sparse.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/operation_sparse.hpp>
#include <boost/numeric/ublas/io.hpp>
#ifdef   _GLUCAT_HAVE_UBLAS_LU_H
#include <boost/numeric/ublas/operation.hpp>
#include <boost/numeric/ublas/lu.hpp>
#endif

#include <stdio.h>
#include <iomanip>

namespace ublas_test
{

  using namespace boost::numeric::ublas;
  using namespace std;
  typedef double Scalar_T;
  typedef row_major                            orientation_t;
  typedef compressed_matrix<Scalar_T, orientation_t > matrix_t;
  typedef matrix_t::size_type         matrix_index_t;

  const double MS_PER_S = 1000.0;

  inline
  double
  elapsed( clock_t cpu_time )
  { return ((clock() - cpu_time)*MS_PER_S) / CLOCKS_PER_SEC; }
}
#endif // TEST_UBLAS_H