Seed seed; Seed skewed_seed;
OLattice<T>& d = const_cast<OSubLattice<T>&>(dd).field(); const Subset& s = dd.subset();
random(OLattice<T>& d, const Subset& s){  Seed seed;  Seed skewed_seed;  const int *tab = s.SiteTable()->slice();  for(int j=0; j < s.NumSiteTable(); ++j)   {    int i = tab[j];    seed = RNG::ran_seed;    skewed_seed.elem() = RNG::ran_seed.elem() * RNG::lattice_ran_mult->elem(i);    fill_random(d.elem(i), seed, skewed_seed, RNG::ran_mult_n);  }  RNG::ran_seed = seed;  // The seed from any site is the same as the new global seed}