void gaussian(OSubLattice<T> dd)
void gaussian(OLattice<T>& d)
void gaussian(OSubLattice<T> dd){  OLattice<T>& d = dd.field();  const Subset& s = dd.subset();  OLattice<T>  r1, r2;  random(r1(s));  random(r2(s));  const int *tab = s.SiteTable()->slice();  for(int j=0; j < s.NumSiteTable(); ++j)   {    int i = tab[j];    fill_gaussian(d.elem(i), r1.elem(i), r2.elem(i));  }}