gaussian(d,s);
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_exit("gaussianUnorderedSubset not implemented"); #endif
void gaussian(const OSubLattice<T,S>& dd){  OLattice<T>& d = const_cast<OSubLattice<T,S>&>(dd).field();  const S& s = dd.subset();  gaussian(d,s);}