int iouter = i >> INNER_LEN; int iinner = i & ((1 << INNER_LEN)-1);
int iouter = i >> INNER_LOG; int iinner = i & ((1 << INNER_LOG)-1);
pokeSite(OLattice<T1>& l, const OScalar<T2>& r, const multi1d<int>& coord){  int i      = Layout::linearSiteIndex(coord);  int iouter = i >> INNER_LEN;  int iinner = i & ((1 << INNER_LEN)-1);  copy_site(l.elem(iouter), iinner, r.elem());  return l;}