gaussian(d(all));
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)); }
void gaussian(OLattice<T>& d){  gaussian(d(all));}