//Boost Test


#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <boost/numeric/ublas/vector_sparse.hpp>
#include <boost/numeric/ublas/vector.hpp>

#include <boost/numeric/ublas/operation.hpp>
#include <boost/numeric/ublas/operation_sparse.hpp>

int main()
{
  using namespace boost::numeric::ublas;

  typedef int                                               VectorIndexType;
  typedef unbounded_array<VectorIndexType>                  IndexArrayType;
  typedef unbounded_array<double>                           ValueArrayType;
  
  typedef compressed_matrix<double,row_major,0,IndexArrayType,ValueArrayType>
                                                            CompressedMatrixType;

  typedef compressed_vector<double,0,IndexArrayType,ValueArrayType>
                                                            BoostProbabilityVector;

  
  CompressedMatrixType m (3, 3, 3 * 3);
  for (VectorIndexType i = 0; i < m.size1 (); ++ i)
      for (VectorIndexType j = 0; j < m.size2 (); ++ j)
          m (i, j) = 3 * i + j;
  std::cout << m << std::endl;

  BoostProbabilityVector v (3, 3);
  for (VectorIndexType i = 0; i < v.size (); ++ i)
      v (i) = i;
      
  std::cout << v << std::endl;

  BoostProbabilityVector u;

  u = prod (m,v);
  std::cout << u << std::endl;

  BoostProbabilityVector w(3);
  //w = m*v
  axpy_prod(m, v, w, true);
  std::cout << w << std::endl;

}