REAL max=0; REAL tmp=0;
T max=0; T tmp=0;
template<class T> T MaxDifference(const CrystVector<T> &a,const CrystVector<T> &b){   const T *p1=a.data();   const T *p2=b.data();   REAL max=0;   REAL tmp=0;   for(long i=0;i<a.numElements();i++)   {      tmp=(T)fabs((REAL) (*p1++ - *p2++));      if(tmp>max) max=tmp;   }   return (T)max;}