#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <complex>

using namespace boost::numeric::ublas;
//using namespace std;

int main(void) {
  
  typedef std::complex<double> dComplex;
  vector<dComplex> v(4);
  //std::cout << v.size() << std::endl;
  //v.resize(4);

  for (int i=0;i<v.size();++i)
    v[i] = dComplex( i, i+1 );
  
  std::cout << v << std::endl;

  std::cout << "L1 norm = " << norm_1(v) << std::endl;
  std::cout << "L2 norm = " << norm_2(v) << std::endl;
  std::cout << "index of Linf norm " << index_norm_inf(v) << std::endl;
  std::cout << "Linf norm = " << norm_inf(v) << " - this is wrong, should |3+4j| = 5" << std::endl;
  std::cout << "Linf norm = " << abs(v[3]) << " - this is the correct answer" << std::endl;
  
  v *= 3.;

  std::cout << std::endl << v << std::endl;

  std::cout << "L1 norm = " << norm_1(v) << std::endl;
  std::cout << "L2 norm = " << norm_2(v) << std::endl;
  std::cout << "index of Linf norm " << index_norm_inf(v) << std::endl;
  std::cout << "Linf norm = " << norm_inf(v) << " - this is wrong, should |9+12j| = 15" << std::endl;
  std::cout << "Linf norm = " << abs(v[3]) << " - this is the correct answer" << std::endl;
  return 0;
}