//

#include "global.hpp"
#include "Point.hpp"
/*
template<int DIM>
   Distance<DIM>::Distance (DOUBLE x, DOUBLE y, DOUBLE z) 
     { assert (DIM > 0 && DIM < 4);           //  Unsupported dimensionality
       value[0] = x;
       if (DIM == 1) {assert(z==0.0 && y == 0.0); return;}
       value[1] = y;
       if (DIM == 2) {assert(z==0.0);             return;}
       value[2] = z;
     }
 */
  template<int DIM>
   Distance<DIM>::Distance<DIM>()  {for (int dm = 0; dm < DIM; ++dm) value[dm]=0;} 

   Distance<1>::Distance<1> (DOUBLE x)         
                {value[0] = x;}
   Distance<2>::Distance<2> (DOUBLE x, DOUBLE y)
                {value[0] = x; value[1] = y;}
   Distance<3>::Distance<3> (DOUBLE x, DOUBLE y, DOUBLE z)
                {value[0] = x; value[1] = y; value[2] = z;}


   DISTANCE<3> Distance<3>::ZERO;
   DISTANCE<2> Distance<2>::ZERO; 
   DISTANCE<1> Distance<1>::ZERO;

   DISTANCE<3> Distance<3>::ONE(1,1,1);
   DISTANCE<2> Distance<2>::ONE(1,1); 
   DISTANCE<1> Distance<1>::ONE(1);