OLattice<T>  r1, r2;
OLattice<T>& d = const_cast<OSubLattice<T,S>&>(dd).field(); const S& s = dd.subset();
void gaussian(OLattice<T>& d, const UnorderedSubset& s){  OLattice<T>  r1, r2;  random(r1,s);  random(r2,s);#if ! defined(QDP_NOT_IMPLEMENTED)  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));  }#else  QDP_error("gaussianUnorderedSubset not implemented");#endif}