normal_multivar.h

#ifndef _NORMAL_MULTIVAR_H_
#define _NORMAL_MULTIVAR_H_

// Boost random generator
#include <boost/random/mersenne_twister.hpp>
#include <boost/random/normal_distribution.hpp>
// Eigen dense matrices
#include <Eigen/Dense>

#include "numeric_utils.h"

namespace numeric_utils {
class NormalMultiVar : public RandomGenerator {
 public:
  NormalMultiVar();

  NormalMultiVar(int seed);

  virtual ~NormalMultiVar(){};

  bool generate(
      Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic>& random_numbers,
      const Eigen::VectorXd& means, const Eigen::MatrixXd& cov,
      unsigned int cases = 1) override;

  std::string name() const override;

 private:
  boost::random::mt19937 generator_; 
  boost::random::normal_distribution<double> distribution_; 
};
}  // namespace numeric_utils

#endif  // _NORMAL_MULTIVAR_H_