template<class T> ostream& operator<<(ostream &os, CrystVector<T> &vect)
template<class T> ostream& operator<<(ostream &os, const CrystMatrix<T> &vect)
template<class T> ostream& operator<<(ostream &os, CrystVector<T> &vect){   //return os << FormatHorizVector(vect);   for(long i=0;i<vect.numElements();i++) os << vect(i) <<endl;   return os;}