for(int i=0; i < ss.numSubsets(); ++i)
for(int k=0; k < ss.numSubsets(); ++k) zero_rep(dest[k]);   const multi1d<int>& lat_color =  ss.latticeColoring();  const int iend = Layout::outerSitesOnNode(); for(int i=0; i < iend; ++i)
sumMulti(const QDPExpr<RHS,OScalar<T> >& s1, const Set& ss){  typename UnaryReturn<OScalar<T>, FnSumMulti>::Type_t  dest(ss.numSubsets());  // lazy - evaluate repeatedly  for(int i=0; i < ss.numSubsets(); ++i)  {    evaluate(dest[i],OpAssign(),s1);  }  return dest;}