#include <iostream>
#include "vector.hh"
#include "vector.hh"


void test_vec_as_const(const Vector& u)
{
  for (unsigned int i = 0; i < u.size(); ++i) {
    std::cout << "u.get_entry(" << i << ") = " << u.get_entry(i);
  }
  std::cout << "u = " << u << std::endl;
  std::cout << "||u||_2 = " << u.l2_norm() << std::endl;

  const Vector v(3, 1.0);

  std::cout << "u.add(v) = " << u.add(v) << std::endl;

  std::cout << "u + v = " << u + v << std::endl;

  for (unsigned int i = 0; i < u.size(); ++i) {
    std::cout << "u[" << i << "] = " << u[i];
  }

  //std::cout << "u - v = " << u - v << std::endl;
  //std::cout << "u*v = " << u*v << std::endl;
}


int main()
{
  Vector u(3, 0.0);
  for (unsigned int i = 0; i < u.size(); ++i) {
    u.set_entry(i, i);
  }

  test_vec_as_const(u);

  for (unsigned int i = 0; i < u.size(); ++i) {
    u[i] = i*i;
  }

  const Vector v(3, 1.0);

  u += v;
  std::cout << u << std::endl;
  (u += v) += v;

  //Vector empty;
  //std::cout << "empty = " << empty << std::endl;
  
  //u.scal(2.5);
  //std::cout << "u = " << u << std::endl;
  //u *= 0.75;
  //std::cout << "u = " << u << std::endl;

  //u.axpy(v, 0.5);
  //std::cout << "u = " << u << std::endl;

  return 0;
}