systematic.hpp

Go to the documentation of this file.

#ifndef SSMPACK_FILTER_RESAMPLER_SYSTEMATIC
#define SSMPACK_FILTER_RESAMPLER_SYSTEMATIC

#include "ssmkit/filter/resampler/base.hpp"
#include "ssmkit/random/generator.hpp"

#include <armadillo>

#include <random>

namespace ssmkit {
namespace filter {
namespace resampler {

template <class Criterion>
class Systematic : public BaseResampler<Systematic<Criterion>> {
  friend class BaseResampler<Systematic<Criterion>>;

 private:
  std::uniform_real_distribution<double> uniform_;

 protected:
  arma::vec generateOrderedNumbers(const int &num_par) {
    double u0 = uniform_(random::Generator::get().getGenerator());
    arma::vec u(num_par);
    int k = 0;
    u.imbue([&u0, &num_par, &k]() { return (k++ + u0) / num_par; });
    return u;
  }

 public:
  Systematic(Criterion criterion)
      : BaseResampler<Systematic<Criterion>>(criterion) {}
};

template<class Criterion>
Systematic<Criterion> makeSystematic(Criterion criterion){
  return Systematic<Criterion>(criterion);
}

} // namespace resampler
} // namespace filter
} // namespace ssmkit
#endif // SSMPACK_FILTER_RESAMPLER_SYSTEMATIC