Vec3d.h

#ifndef QNANO_NEW_CLASS_VEC3D_
#define QNANO_NEW_CLASS_VEC3D_

#include <cmath>
#include <ostream>
#include <istream>

class Vec3d{
public:

  double v[3];
 

//Functions:

  inline void set(double x, double y,double z){v[0]=x;v[1]=y;v[2]=z;}
  inline void set(const Vec3d &other){v[0]=other.v[0];v[1]=other.v[1];v[2]=other.v[2];}

  inline double norm2() const{return (v[0]*v[0]+v[1]*v[1]+v[2]*v[2]);}
  inline double norm() const{return sqrt(norm2());}
  inline void normalize(){double n=norm(); v[0]/=n; v[1]/=n; v[2]/=n;}

//Operators:
  double & operator[](const size_t i){return v[i];}
  double operator[](const size_t i) const{return v[i];}
  Vec3d & operator=(const Vec3d &v2);
  Vec3d & operator+=(const Vec3d &other);
  Vec3d & operator-=(const Vec3d &other);
  Vec3d & operator*=(double d);
  Vec3d & operator/=(double d);
  friend std::ostream& operator<<(std::ostream& os, const Vec3d& v);
  
//Constructors:
  Vec3d();
  Vec3d(double x, double y,double z);
  Vec3d(const Vec3d &v2);
};

std::ostream& operator<<(std::ostream& os, const Vec3d& v);
Vec3d operator+(const Vec3d left,const Vec3d right);
Vec3d operator-(const Vec3d left,const Vec3d right);

Vec3d operator*(double left,const Vec3d right);
Vec3d operator*(const Vec3d left, double right);
Vec3d operator*(int left,const Vec3d right);
Vec3d operator*(const Vec3d left, int right);
Vec3d operator/(const Vec3d left, double right);

double operator*(const Vec3d & v1, const Vec3d & v2);
Vec3d cross(const Vec3d &v1, const Vec3d &v2);
double norm(const Vec3d &v);
double dist(const Vec3d &v1, const Vec3d &v2);

std::istream& operator>>(std::istream& is, Vec3d &v);



#endif