Scythe-1.0.3
template<class RNGTYPE>
template<matrix_order O, matrix_style S>
Matrix<double, O, Concrete> scythe::rng< RNGTYPE >::rdirich ( const Matrix< double, O, S > &  alpha) [inline]

Generate a Dirichlet distributed random variate Matrix.

This function returns a pseudo-random matrix-valued variate drawn from the Dirichlet disribution described by the vector alpha.

Parameters:
alphaA vector of non-negative reals.
Exceptions:
scythe_invalid_arg(Level 1)
scythe_dimension_error(Level 1)