// vec.cpp

#include <iostream>
#include <cmath>
#include "vec.h"
using namespace std;

template <typename T>
vec<T>::vec(const T& x, const T& y, const T& z)
      : _x(x), _y(y), _z(z) {}

template <typename T>
vec<T>::vec(const vec<T>& v)
      : _x(v._x), _y(v._y), _z(v._z) {}

template <typename T>
vec<T> operator+(const vec<T>& v1, const vec<T>& v2)
{
  vec<T> temp(v1);
  temp._x += v2._x;
  temp._y += v2._y;
  temp._z += v2._z;
  return temp;
}

template <typename T>
T vec<T>::size()
{
  return sqrt(_x*_x + _y*_y + _z*_z);
}

template <typename T>
void vec<T>::show()
{
  cout << "\n  vec: (" << _x << ", " << _y << ", " << _z << ")";
}

#ifdef TEST_VEC

int main()
{
  vec<double> v1(1,1,1);
  vec<double> v2(2,2,2);

  vec<double> v3 = v1 + v2;
  v3.show();
  cout << "\n  size = " << v3.size();

  vec<int> v4(1,1,1);
  vec<int> v5(-2,-2,-2);
  vec<int> v6 = v4 + v5;
  v6.show();
  cout << "\n  size = " << v6.size();

  cout << "\n\n";
}

#endif