//  NumVec class implementation

//  See NumVec.hpp for details

#include "global.hpp"
#include "NumVec.hpp"

NumVec<double> dummy(1);                //  Force instantiation

template<class T>
  void NumVec<T>::copy(const NumVec<T>& rs) 
  {assert(dim() == rs.dim());
   for (sizet  count = 0; count <= dim(); count++)
      *(data+count) = *(rs.data+count);
   return;
  }

template<class T>                   //  
  void NumVec<T>::copy(T rs) 
  {for (sizet count  = 0; count <= dim(); count++)
     *(data+count) = rs;
   return;
  }

template<class T>                   //  
  ostream& NumVec<T>::put(ostream& s) const
  {sizet  ctr = lbound();
   s << data[ctr++];

   while (ctr <= hbound())
     s << ", " << data[ctr++];
   return s;
  }