#ifndef _vec_arithmetic_H_
#define _vec_arithmetic_H

/* Define all operators after definition of  basic operators */

#define vec_arithmetic1(V) \
friend V operator-( const V &a ){ return V(0) - a; } \
friend V operator+( const V &a ){ return a; } 

#define vec_arithmetic(V,S) \
friend V operator+( const V &a, const S &b ){ return a + V( b ); } \
friend V operator-( const V &a, const S &b ){ return a - V( b ); } \
friend V operator*( const V &a, const S &b ){ return a * V( b ); } \
friend V operator/( const V &a, const S &b ){ return a / V( b ); } \
friend V operator+( const S &a, const V &b ){ return V(a)+b; } \
friend V operator-( const S &a, const V &b ){ return V(a)-b; } \
friend V operator*( const S &a, const V &b ){ return V(a)*b; } \
friend V operator/( const S &a, const V &b ){ return V(a)/b; } 
#endif