evaluate(dest[i],OpAssign(),s1);
int j = lat_color[i]; dest[j].elem() += forEach(s1, EvalLeaf1(i), OpCombine());
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;}