gamma

ガンマ分布の乱数を生成します。

説明

gamma クラス・オブジェクトは、generate 関数で使用され、形状 \(\alpha\)、変位 \(a\)、スケール・パラメーター \(\beta\) を持つガンマ分布の乱数を生成します。ここで、\(a, \alpha, \beta \in R; \alpha > 0; \beta > 0\) です。

確率分布は次のように表されます。

\[\begin{split}f_{a, \alpha, \beta}(x) = \left\{ \begin{array}{rcl} \frac{1}{\Gamma(\alpha)\beta^{\alpha}}(x - a)^{\alpha - 1}e^{-(x - a) / \beta}, x \ge a \\ 0, x < a \end{array}\right.\end{split}\]

累積分布関数は次のようになります。

\[\begin{split}F_{a, \alpha, \beta}(x) = \left\{ \begin{array}{rcl} \int^x_a\frac{1}{\Gamma(\alpha)\beta^{\alpha}}(y - a)^{\alpha - 1}e^{-(y - a) / \beta}dy, x \ge a \\ 0, x < a \end{array}\right.\end{split}\]

API

構文

namespace oneapi::mkl::rng::device { 
  template<typename RealType, typename Method> 
  class gamma { 
  public: 
    using method_type = Method; 
    using result_type = RealType; 

    gamma() : gamma((RealType)1.0, (RealType)0.0, (RealType)1.0){} 
    explicit gamma(RealType alpha, RealType a, RealType beta); 
    explicit gamma(const param_type& pt); 

    RealType alpha() const; 
    RealType a() const; 
    RealType beta() const; 
    param_type param() const; 
    std::size_t count_rejected_numbers() const; 
    void param(const param_type& pt); 
  }; 
}

インクルード・ファイル

• oneapi/mkl/rng/device.hpp

テンプレート・パラメーター

typename Type	生成された値のタイプ。固有の値は次のとおりです: float double
typename Method	生成メソッド。固有の値は次のとおりです: oneapi::mkl::rng::device::gamma_method::by_default oneapi::mkl::rng::device::gamma_method::marsaglia oneapi::mkl::rng::device::gamma_method::marsaglia_accurate 配布テンプレート・パラメーター・メソッドの説明を参照してください。

注

ベクトル数学関数に基づく追加の最適化を有効にするには、MKL_RNG_USE_BINARY_CODE マクロを定義し、インテル® oneAPI マス・カーネル・ライブラリー・リンクラインのアドバイザー (英語) で生成された出力に従って oneMKL 静的ライブラリーをリンクします。一部のハードウェアではパフォーマンスが向上する可能性があります。

入力パラメーター

名前	タイプ	説明
alpha	RealType (float, double)	形状 \(\alpha\)。
a	RealType (float, double)	変位 a。
beta	RealType (float, double)	スケール係数 \(\beta\)。

特別なメソッド

std::size_t count_rejected_numbers() const

最後の generate 関数呼び出し中に拒否された乱数の数を返します。generate 呼び出しがない場合、0 を返します。