#include "Rn.hpp"

std::ostream& operator<<(std::ostream& os, const Rn& r)
{ 
	os<<"< "<<r.vect[0];
	for(int i=1;i<r.dim;++i)
		os<<" , "<<r.vect[i];
	os<<">"<<std::endl; 
	return os;
}	

Rn operator+(const Rn& u,const Rn& v)  // friend function
{
	Rn w(u.dim);
	for(int i=0;i<u.dim;++i)
	    w[i]=u[i]+v[i];
	return w;
}

Rn operator-(const Rn& u,const Rn& v)  // friend function
{
	Rn w(u.dim);
	for(int i=0;i<u.dim;++i)
	    w[i]=u[i]-v[i];
	return w;
} 

Rn Rn::operator+(const Rn& v)  // somme
{
	for(int i=0;i<dim;++i)
	    vect[i]=vect[i]+v[i];
	return *this;
}

Rn Rn::operator-(const Rn& v)  // difference
{
	Rn r(*this);
	r-=v;
	return r;	
}

Rn& Rn::operator=(const Rn& v)
{
	if(&v != this)
	{
		for(int i=0;i<dim;++i)
		    vect[i]=v[i];
	}
	return *this;
}

Rn& Rn::operator+=(const Rn& v)
{
	for(int i=0;i<dim;++i)
	    	vect[i]+=v[i];
	return *this;
}

Rn& Rn::operator-=(const Rn& v)
{
	for(int i=0;i<dim;++i)
	    	vect[i]-= v[i];
	return *this;
}	

double Rn::operator*(const Rn& v)const
{
	double prod=0;
	for(int i=0;i<dim;++i)
   		prod=prod+vect[i]*v[i];
	return prod;
}

Rn& Rn::operator*(const double& d)	  // produit par un scalaire	
{
	for(int i=0;i<dim;++i)
   		vect[i]*=d;
        return *this;
}

Rn& Rn::operator*=(const double& d)	  // produit par un scalaire
{
	for(int i=0;i<dim;++i)
   		vect[i]*=d;
        return *this;
}

Rn& operator*(const Rn& v, const double& d)
{
	return v*d;
}

Rn Rn::operator-() const
{
	for(int i=0;i<dim;++i)
   		vect[i]=-vect[i];
        return *this;
}