More sophisticated plateau finder

[carstenbauer]

Maybe take a look at the ALPS source code. But don't expect too much :)

[carstenbauer]

Regarding ALPS, alps/src/alps/alea/simplebinning.h might be the file of interest. It contains:

template <class T>
typename SimpleBinning<T>::convergence_type SimpleBinning<T>::converged_errors() const
{
  convergence_type conv;
  result_type err=error();
  resize_same_as(conv,err);
  const unsigned int range=4;
  typename slice_index<convergence_type>::type it;
  if (binning_depth()<range) {
    for (it= slices(conv).first; it!= slices(conv).second; ++it)
      slice_value(conv,it) = MAYBE_CONVERGED;
  }
  else {
    for (it= slices(conv).first; it!= slices(conv).second; ++it)
      slice_value(conv,it) = CONVERGED;

    for (unsigned int i=binning_depth()-range;i<binning_depth()-1;++i) {
      result_type this_err(error(i));
      for (it= slices(conv).first; it!= slices(conv).second; ++it)
        if (std::abs(slice_value(this_err,it)) >= std::abs(slice_value(err,it)))
          slice_value(conv,it)=CONVERGED;
        else if (std::abs(slice_value(this_err,it)) < 0.824 * std::abs(slice_value(err,it)))
          slice_value(conv,it)=NOT_CONVERGED;
        else if (std::abs(slice_value(this_err,it)) <0.9* std::abs(slice_value(err,it))  &&
            slice_value(conv,it)!=NOT_CONVERGED)
          slice_value(conv,it)=MAYBE_CONVERGED;
    }
  }
  return conv;
}