int j = lat_color[i]; dest[j].elem() += forEach(s1, EvalLeaf1(i), OpCombine());
evaluate(dest[i],OpAssign(),s1);
sumMulti(const QDPExpr<RHS,OLattice<T> >& s1, const Set& ss){  typename UnaryReturn<OLattice<T>, FnSumMulti>::Type_t  dest(ss.NumSubsets());  // Initialize result with zero  for(int k=0; k < ss.NumSubsets(); ++k)    zero_rep(dest[k]);  // Loop over all sites and accumulate based on the coloring   const multi1d<int>& lat_color =  ss.LatticeColoring();  for(int i=0; i < Layout::vol(); ++i)   {    int j = lat_color[i];    dest[j].elem() += forEach(s1, EvalLeaf1(i), OpCombine());   // SINGLE NODE VERSION FOR NOW  }  return dest;}