//  Velocity class template
//  -----------------------
#ifndef VELOCITY
#define VELOCITY const Velocity
template<int DIM> class Distance;

template<int DIM>
    class  Velocity {

friend class Distance<DIM>;
      
    double value;

static char*  units;
public:
    Velocity(DOUBLE v) : value(v) {}


#define Class Velocity
#include "Additive.hpp"
 
   
};

template<int DIM>
 inline Velocity<DIM>& Velocity<DIM>::operator+=(VELOCITY<DIM> rs)
        {value+=rs.value; return *this;}
template<int DIM>
 inline Velocity<DIM>& Velocity<DIM>::operator-=(VELOCITY<DIM> rs)
        {value-=rs.value; return *this;}

template<int DIM>
 inline Velocity<DIM>& Velocity<DIM>::operator*=(DOUBLE        rs)
        {value*=rs;       return *this;}
template<int DIM>
 inline Velocity<DIM>& Velocity<DIM>::operator/=(DOUBLE        rs)
        {value/=rs;       return *this;}
     

#undef Class
#endif