skip_ahead#

スキップアヘッド方式を使用してエンジンの状態を進めます。skip_ahead 関数は、スキップアヘッド方式を適用するため次のインターフェイスをサポートしています。

共通インターフェイス
スキップされた要素のパーティション化された数とのインターフェイス

説明 #

skip_ahead 関数は、エンジンによって提供されるランダムシーケンス内の指定された数の要素をスキップします。この機能は、元のエンジンからの乱数を異なる計算ノードに分散させる場合に特に役立ちます。計算ノードで使用される乱数の最大数が num_to_skip である場合、元の乱数シーケンスは skip_ahead によって num_to_skip サイズの重複しないブロックに分割され、各ブロックがそれぞれの計算ノードに対応するようになります。計算ノードの数は無制限です。この方法は、ブロック分割方式またはスキップアヘッド方式として知られています (次の図を参照)。

ブロック分割法

スキップアヘッド法は、スキップアヘッドによる要素のスキップを許可する基本生成器でのみサポートされます。これは、生成器で要素を単純に生成してその後手動でスキップするよりも効率的です。詳細については VS ノートを参照してください。

準ランダム基本生成器の場合、スキップアヘッド法は、準ランダムベクトル全体ではなく、準ランダムベクトルのコンポーネントで機能することに注意してください。したがって、NS 準ランダムベクトルをスキップするには、num_to_skip パラメーターを num_to_skip *dim に設定します。ここで、dim は準ランダムベクトルの次元です。

スキップされた要素の数が $2^{63}$ より大きい場合、スキップされた要素の分割された数を持つインターフェイスが使用されます。

関数の事前呼び出しは、以下に示すように、サイズ n のリストでスキップされた要素の数を表します。

num_to_skip[0]+ num_to_skip[1]*2⁶⁴+ num_to_skip[2]* 2¹²⁸+ … +num_to_skip[n-1]*2^64*(n-1) ;

スキップされる要素の数が $2^{63}$ 未満の場合、両方のインターフェイスを使用できます。

次のコードは、skip_ahead 関数を使用して 3 つの独立したストリームを初期化する方法を示しています。

ブロック分割法のコード#

     ...
  // 3 つの同一のエンジンを作成 
  oneapi::mkl::rng::mcg31m1 engine_1(queue, seed); 
  oneapi::mkl::rng::mcg31m1 engine_2(engine_1); 
  oneapi::mkl::rng::mcg31m1 engine_3(engine_2); 
 
 
  // 7 要素先へスキップする 2 番目のエンジン 
  oneapi::mkl::rng::skip_ahead(engine_2, 7); 
 
 
 // 14 要素先へスキップする 3 番目のエンジン 
 oneapi::mkl::rng::skip_ahead(engine_3, 14); 
 ...

パーティション化された要素数を持つブロック分割法のコード#

  // 最初のエンジンの作成 
  oneapi::mkl::rng::mrg32k3a engine_1(queue, seed); 
 
 
  // ランダムストリーム内の 2^64 個の要素をスキップするには、スキップする要素の数は 
  //  num_to_skip = 2^64 = 0 + 1 * 2^64 と表す必要があります。
  std::initializer_list<std::uint64_t> num_to_skip = {0, 1}; 
 
 
 
 
 // 1 番目のエンジンに基づいて 2 番目のエンジンを作成。 2^64 スキップ 
 oneapi::mkl::rng::mrg32k3a engine_2(engine_1); 
 oneapi::mkl::rng::skip_ahead(engine_2, num_to_skip); 
 
 
 // 2 番目のエンジンをベースに 3 番目のエンジンを作成。 2^64 スキップ 
 oneapi::mkl::rng::mrg32k3a engine_3(engine_2); 
 oneapi::mkl::rng::skip_ahead(engine_3, num_to_skip); 
 ...

API #

構文#

共通インターフェイス

namespace oneapi::mkl::rng { 
  template<typename Engine> 
  void skip_ahead (Engine& engine, 
                   std::uint64_t num_to_skip) 
}

スキップされた要素のパーティション数を持つインターフェイス

namespace oneapi::mkl::rng { 
  template<typename Engine> 
  void skip_ahead (Engine& engine, 
                   std::initializer_list<std::uint64_t> num_to_skip) 
}

インクルード・ファイル#

oneapi/mkl/rng.hpp

入力パラメーター#