#include <iostream>
#include <valarray>
#include <boost/timer.hpp>
#include <boost/numeric/ublas/vector.hpp>

const size_t size = 1024;
const int iterations = 1000000;
typedef double RealType;

template<int iterations, class vector_t>
void arithmetic_benchmark(vector_t& v)
{
  {
    vector_t tmp = v;
    for ( int i = 0; i < iterations; ++i)
    {
      tmp += v;
      tmp -= v;
    }
  }
}

int main()
{
  {
    std::valarray<RealType> v(size);
    boost::timer timer;
    arithmetic_benchmark<iterations>(v);
    std::cout << timer.elapsed() << std::endl;
  }

  {
    boost::numeric::ublas::vector<RealType> v;
    boost::timer timer;
    arithmetic_benchmark<iterations>(v);
    std::cout << timer.elapsed() << std::endl;
  }

  return 0;
}