int numSubsets() const {return Layout::lattSize()[dir_decay];}
TimeSliceFunc::numSubsets() const { if ((dir_decay<0)||(dir_decay>=Nd)) { return 1 ; } else { return Layout::lattSize()[dir_decay] ; } }
  int numSubsets() const {return Layout::lattSize()[dir_decay];}