operator*(const GammaConst<4,15>&, const PSpinMatrix<T2,4>& r)
operator*(const PSpinMatrix<T2,4>& l, const GammaConst<4,15>&)
operator*(const GammaConst<4,15>&, const PSpinMatrix<T2,4>& r){  typename BinaryReturn<GammaConst<4,15>, PSpinMatrix<T2,4>, OpGammaConstMultiply>::Type_t  d;  for(int i=0; i < 4; ++i)  {    d.elem(0,i) = r.elem(0,i);    d.elem(1,i) = r.elem(1,i);    d.elem(2,i) = -r.elem(2,i);    d.elem(3,i) = -r.elem(3,i);  }    return d;}