random_dirichlet.h

// Description: Dirichlet-distributed random numbers

#ifndef PASTELSYS_RANDOM_DIRICHLET_H
#define PASTELSYS_RANDOM_DIRICHLET_H

#include "pastel/sys/mytypes.h"
#include "pastel/sys/vector.h"

namespace Pastel
{

    //! Generates a vector from the Dirichlet distribution.
    /*!
   Preconditions:
   allGreater(shape, 0)

   Dirichlet distribution is a multi-dimensional 
   generalization of the beta-distribution.
   */
    template <typename Real, integer N>
    Vector<Real, ModifyN<N, N - 1>::Result> randomDirichlet(
        const Vector<Real, N>& shape);

}

#include "pastel/sys/random/random_dirichlet.hpp"

#endif