numeric_utils.h

#ifndef _NUMERIC_UTILS_H_
#define _NUMERIC_UTILS_H_

#include <ctime>
#include <utility>
#include <vector>
#include <Eigen/Dense>

namespace numeric_utils {

Eigen::MatrixXd corr_to_cov(const Eigen::MatrixXd& corr,
                const Eigen::VectorXd& std_dev);

bool convolve_1d(const std::vector<double>& input_x,
                 const std::vector<double>& input_y,
                 std::vector<double>& response);

double trapazoid_rule(const std::vector<double>& input_vector, double spacing);

double trapazoid_rule(const Eigen::VectorXd& input_vector, double spacing);

class RandomGenerator {
 public:
  RandomGenerator() = default;

  virtual ~RandomGenerator() {};

  RandomGenerator(const RandomGenerator&) = delete;

  RandomGenerator& operator=(const RandomGenerator&) = delete;

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

   virtual std::string name() const = 0;
  
 protected:
  int seed_ = static_cast<int>(
      std::time(nullptr)); 
};
}  // namespace numeric_utils

#endif  // _NUMERIC_UTILS_H_