gausSmear<LatticePropagator>(u, chi, width, ItrGaus, j_decay);
T psi;  Real ftmp = - (width*width) / Real(4*ItrGaus);   Real ftmpi = Real(1) / ftmp;  for(int n = 0; n < ItrGaus; ++n) { psi = chi * ftmp; klein_gord(u, psi, chi, ftmpi, j_decay); }
void gausSmear(const multi1d<LatticeColorMatrix>& u, 	       LatticePropagator& chi, 	       const Real& width, int ItrGaus, int j_decay){  gausSmear<LatticePropagator>(u, chi, width, ItrGaus, j_decay);}