OLattice<T>& d = const_cast<OSubLattice<T,S>&>(dd).field(); const S& s = dd.subset();
OLattice<T>  r1, r2;
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);}