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
gaussian(d,s);
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}