/** Utilities for random values necessary in MOMULPER.
 *
 */
#ifndef RANDOM_HPP_
#define RANDOM_HPP_

// TODO: all random value generation through these functions.

#include<vector>
#include<random>

using namespace std;

/** Randomize vector with values that sum up to 1. */
void randomize1(vector<double> &vec, mt19937 *mt);

#endif