JP3123920B2 - マルチプロセッサシステムにおける乱数発生装置および乱数発生方法 - Google Patents
マルチプロセッサシステムにおける乱数発生装置および乱数発生方法Info
- Publication number
- JP3123920B2 JP3123920B2 JP08083416A JP8341696A JP3123920B2 JP 3123920 B2 JP3123920 B2 JP 3123920B2 JP 08083416 A JP08083416 A JP 08083416A JP 8341696 A JP8341696 A JP 8341696A JP 3123920 B2 JP3123920 B2 JP 3123920B2
- Authority
- JP
- Japan
- Prior art keywords
- generated
- random
- processor
- random numbers
- processors
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Multi Processors (AREA)
- Complex Calculations (AREA)
Description
【発明の属する技術分野】本発明は、マルチプロセッサ
システムを構成する各プロセッサに乱数を発生させるマ
ルチプロセッサシステムにおける乱数発生装置および乱
数発生方法に関する。
システムを構成する各プロセッサに乱数を発生させるマ
ルチプロセッサシステムにおける乱数発生装置および乱
数発生方法に関する。
【０００１】
【従来の技術】「Lehmer D.H., Proceedings of a Secon
d Symposium on Large-Scale DigitalCalculating Mach
inery, Ann. Comput. Lab. Harvard U. 26 (1951), pp.
141-146」には、乱数発生法の一つである合同法が記載さ
れている。
d Symposium on Large-Scale DigitalCalculating Mach
inery, Ann. Comput. Lab. Harvard U. 26 (1951), pp.
141-146」には、乱数発生法の一つである合同法が記載さ
れている。
【０００２】合同法は、合同式を用いた漸化式 Ｘi+1 ＝ ａＸi＋ｃ （ｍｏｄ ｍ） すなわち、ある数Ｘiにある定数ａを乗算してある数ｃ
を加算したものを、ある数ｍで除算した余りをＸi+1と
するという整数演算を行う漸化式によって、周期的に反
復する乱数の数列｛Ｘi｝を求める方法である。
を加算したものを、ある数ｍで除算した余りをＸi+1と
するという整数演算を行う漸化式によって、周期的に反
復する乱数の数列｛Ｘi｝を求める方法である。
【０００３】この漸化式においてｃ＝０の場合、漸化式 Ｘi+1 ＝ ａＸi （ｍｏｄ ｍ） の右辺が乗算だけになるので、特に乗算合同法と呼び、
ｃ≠０の場合には混合合同法と呼ぶ。
ｃ≠０の場合には混合合同法と呼ぶ。
【０００４】特開昭６０−１３６８７６号公報には、合
同法のアルゴリズムに基づいて、並列に乱数を発生させ
る方法が記載されている。
同法のアルゴリズムに基づいて、並列に乱数を発生させ
る方法が記載されている。
【０００５】
【発明が解決しようとする課題】ところで汎用的に利用
されるマルチプロセッサシステムにおいて乱数を発生す
る場合には、該乱数を発生するプログラム、発生された
乱数を使用する汎用プログラム、さらにこれらとは全く
関係のない他のプログラム群がマルチプロセッサシステ
ムを構成するプロセッサ群を共有して利用する。このよ
うな場合、乱数を発生するプログラムが利用することが
できるプロセッサの数は動的に変化する。
されるマルチプロセッサシステムにおいて乱数を発生す
る場合には、該乱数を発生するプログラム、発生された
乱数を使用する汎用プログラム、さらにこれらとは全く
関係のない他のプログラム群がマルチプロセッサシステ
ムを構成するプロセッサ群を共有して利用する。このよ
うな場合、乱数を発生するプログラムが利用することが
できるプロセッサの数は動的に変化する。
【０００６】特開昭６０−１３６８７６号公報において
は、乱数発生開始後に乱数を発生させるプロセッサの数
を変更することについては、示唆されておらず、このよ
うな汎用的に利用されるマルチプロセッサシステムにお
いては、効率良く乱数を発生することができないという
問題点がある。
は、乱数発生開始後に乱数を発生させるプロセッサの数
を変更することについては、示唆されておらず、このよ
うな汎用的に利用されるマルチプロセッサシステムにお
いては、効率良く乱数を発生することができないという
問題点がある。
【０００７】本発明の目的は、マルチプロセッサシステ
ムにおいて、容易に乱数を発生することができるように
することにある。
ムにおいて、容易に乱数を発生することができるように
することにある。
【０００８】本発明の他の目的は、マルチプロセッサシ
ステムにおいて、乱数発生中に乱数を発生させるプロセ
ッサの数を変更することができるようにすることにあ
る。
ステムにおいて、乱数発生中に乱数を発生させるプロセ
ッサの数を変更することができるようにすることにあ
る。
【０００９】本発明の他の目的は、マルチプロセッサシ
ステムにおいて実行されるさまざまなプログラム群と共
存して効率良く乱数を発生することができるようにする
ことにある。
ステムにおいて実行されるさまざまなプログラム群と共
存して効率良く乱数を発生することができるようにする
ことにある。
【００１０】
【課題を解決するための手段】本発明の第１のマルチプ
ロセッサシステムにおける乱数発生装置は、複数のプロ
セッサのそれぞれにより発生された乱数列を組み合わせ
た乱数列が一つのプロセッサにより発生された乱数列と
同一になるように、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭただ
し、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎはプロ
セッサの番号、**はべき算からなる式により算出した初
期値と乱数発生数＝（発生総数／プロセッサ数）−１た
だし、乱数発生数はそれぞれのプロセッサに発生させる
乱数の個数、発生総数は全部のプロセッサに発生させる
乱数の総数、プロセッサは乱数を発生させるプロセッサ
数からなる式により算出した乱数発生数とを設定して、
該初期値と該乱数発生数とに基づいてＸn+rt＝（α**
ｒ）＊Ｘn+r(t-1)ただし、Ｘは整数、αは初期値、ｒは
プロセッサ数、ｎはプロセッサの番号、ｔは自然数で、
１≦ｔ≦乱数発生数、**はべき算からなる式により乱数
列を前記複数のプロセッサのそれぞれに発生させ、Ｘn
を発生乱数列格納部内のｎ番目の領域に格納し、Ｘn+rt
を該発生乱数列格納部内のn+rt番目の領域に格納するこ
とによって、一つの乱数列にするプロセッサ管理部から
構成されている。
ロセッサシステムにおける乱数発生装置は、複数のプロ
セッサのそれぞれにより発生された乱数列を組み合わせ
た乱数列が一つのプロセッサにより発生された乱数列と
同一になるように、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭただ
し、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎはプロ
セッサの番号、**はべき算からなる式により算出した初
期値と乱数発生数＝（発生総数／プロセッサ数）−１た
だし、乱数発生数はそれぞれのプロセッサに発生させる
乱数の個数、発生総数は全部のプロセッサに発生させる
乱数の総数、プロセッサは乱数を発生させるプロセッサ
数からなる式により算出した乱数発生数とを設定して、
該初期値と該乱数発生数とに基づいてＸn+rt＝（α**
ｒ）＊Ｘn+r(t-1)ただし、Ｘは整数、αは初期値、ｒは
プロセッサ数、ｎはプロセッサの番号、ｔは自然数で、
１≦ｔ≦乱数発生数、**はべき算からなる式により乱数
列を前記複数のプロセッサのそれぞれに発生させ、Ｘn
を発生乱数列格納部内のｎ番目の領域に格納し、Ｘn+rt
を該発生乱数列格納部内のn+rt番目の領域に格納するこ
とによって、一つの乱数列にするプロセッサ管理部から
構成されている。
【００１１】本発明の第２のマルチプロセッサシステム
における乱数発生装置は、複数のプロセッサのそれぞれ
により発生された乱数列を組み合わせた乱数列が１つの
プロセッサにより発生された乱数列と同一になるよう
に、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭただし、α、Ｘ、お
よびＭは整数で、０≦Ｘ＜Ｍ、ｎはプロセッサの番号、
**はべき算からなる式により算出した初期値と乱数発生
数＝（発生総数／プロセッサ数）−１ただし、乱数発生
数はそれぞれのプロセッサに発生させる乱数の数、発生
総数は全部のプロセッサに発生させる乱数の総数、プロ
セッサは乱数を発生させるプロセッサ数からなる式によ
り算出した乱数発生数とを設定して、該初期値と該乱数
発生数とに基づいて、Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1)
ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式により乱数列を前記複数のプ
ロセッサのそれぞれに発生させている最中に乱数を発生
させるプロセッサを追加し、追加したプロセッサを含め
る前記複数のプロセッサのそれぞれにＸq+n＝Ｘq＊（α
**ｎ）ｍｏｄＭただし、α、Ｘ、およびＭは整数で、０
≦Ｘ＜Ｍ、ｎはプロセッサの番号、ｑは既設のすべての
プロセッサが発生した乱数の合計、**はべき算からなる
式により算出した初期値と乱数発生数＝｛（発生総数−
発生済数）／プロセッサ数｝−１ただし、乱数発生数は
それぞれのプロセッサに発生させる乱数の数、発生総数
は全部のプロセッサに発生させる乱数の総数、発生済数
は発生済みの乱数の個数、プロセッサは乱数を発生させ
るプロセッサ数からなる式により算出した乱数発生数と
を再設定して、該初期値と該乱数発生数とに基づいてＸ
n+q+rt＝（α**ｒ）＊Ｘn+q+r(t-1)ただし、Ｘは整数、
αは初期値、ｒはプロセッサ数、ｎはプロセッサの番
号、ｑは既設のすべてのプロセッサが発生した乱数の合
計、ｔは自然数で、１≦ｔ≦乱数発生数、**はべき算か
らなる式による乱数列を発生させ、Ｘnを発生乱数列格
納部内のｎ番目の領域に格納し、Ｘn+rtを該発生乱数列
格納部内のｎ＋ｒｔ番目の領域に格納して、さらに、Ｘ
q+nを該発生乱数列格納部内のｑ＋ｎ番目の領域に格納
し、Ｘn+q+rtを該発生乱数列格納部内のｎ＋ｑ＋ｒｔ番
目の領域に格納することによって、一つの乱数列にする
プロセッサ管理部から構成されている。
における乱数発生装置は、複数のプロセッサのそれぞれ
により発生された乱数列を組み合わせた乱数列が１つの
プロセッサにより発生された乱数列と同一になるよう
に、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭただし、α、Ｘ、お
よびＭは整数で、０≦Ｘ＜Ｍ、ｎはプロセッサの番号、
**はべき算からなる式により算出した初期値と乱数発生
数＝（発生総数／プロセッサ数）−１ただし、乱数発生
数はそれぞれのプロセッサに発生させる乱数の数、発生
総数は全部のプロセッサに発生させる乱数の総数、プロ
セッサは乱数を発生させるプロセッサ数からなる式によ
り算出した乱数発生数とを設定して、該初期値と該乱数
発生数とに基づいて、Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1)
ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式により乱数列を前記複数のプ
ロセッサのそれぞれに発生させている最中に乱数を発生
させるプロセッサを追加し、追加したプロセッサを含め
る前記複数のプロセッサのそれぞれにＸq+n＝Ｘq＊（α
**ｎ）ｍｏｄＭただし、α、Ｘ、およびＭは整数で、０
≦Ｘ＜Ｍ、ｎはプロセッサの番号、ｑは既設のすべての
プロセッサが発生した乱数の合計、**はべき算からなる
式により算出した初期値と乱数発生数＝｛（発生総数−
発生済数）／プロセッサ数｝−１ただし、乱数発生数は
それぞれのプロセッサに発生させる乱数の数、発生総数
は全部のプロセッサに発生させる乱数の総数、発生済数
は発生済みの乱数の個数、プロセッサは乱数を発生させ
るプロセッサ数からなる式により算出した乱数発生数と
を再設定して、該初期値と該乱数発生数とに基づいてＸ
n+q+rt＝（α**ｒ）＊Ｘn+q+r(t-1)ただし、Ｘは整数、
αは初期値、ｒはプロセッサ数、ｎはプロセッサの番
号、ｑは既設のすべてのプロセッサが発生した乱数の合
計、ｔは自然数で、１≦ｔ≦乱数発生数、**はべき算か
らなる式による乱数列を発生させ、Ｘnを発生乱数列格
納部内のｎ番目の領域に格納し、Ｘn+rtを該発生乱数列
格納部内のｎ＋ｒｔ番目の領域に格納して、さらに、Ｘ
q+nを該発生乱数列格納部内のｑ＋ｎ番目の領域に格納
し、Ｘn+q+rtを該発生乱数列格納部内のｎ＋ｑ＋ｒｔ番
目の領域に格納することによって、一つの乱数列にする
プロセッサ管理部から構成されている。
【００１２】本発明の第１のマルチプロセッサシステム
における乱数発生方法は、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄ
Ｍただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎ
はプロセッサの番号、**はべき算からなる式により算出
した初期値と乱数発生数＝（発生総数／プロセッサ数）
−１ただし、乱数発生数はそれぞれのプロセッサに発生
させる乱数の個数、発生総数は全部のプロセッサに発生
させる乱数の総数、プロセッサは乱数を発生させるプロ
セッサ数からなる式により算出した乱数発生数とを設定
する設定ステップと、該初期値と該乱数発生数とに基づ
いてＸn+rt＝（α**ｒ）＊Ｘn+r(t-1)ただし、Ｘは整
数、αは初期値、ｒはプロセッサ数、ｎはプロセッサの
番号、ｔは自然数で、１≦ｔ≦乱数発生数、**はべき算
からなる式により乱数列を前記複数のプロセッサのそれ
ぞれに発生させる発生ステップと、Ｘnを発生乱数列格
納部内のｎ番目の領域に格納し、Ｘn+rtを該発生乱数列
格納部内のｎ＋ｒｔ番目の領域に格納することによっ
て、一つの乱数列にする格納ステップとを含むことを特
徴とする。
における乱数発生方法は、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄ
Ｍただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎ
はプロセッサの番号、**はべき算からなる式により算出
した初期値と乱数発生数＝（発生総数／プロセッサ数）
−１ただし、乱数発生数はそれぞれのプロセッサに発生
させる乱数の個数、発生総数は全部のプロセッサに発生
させる乱数の総数、プロセッサは乱数を発生させるプロ
セッサ数からなる式により算出した乱数発生数とを設定
する設定ステップと、該初期値と該乱数発生数とに基づ
いてＸn+rt＝（α**ｒ）＊Ｘn+r(t-1)ただし、Ｘは整
数、αは初期値、ｒはプロセッサ数、ｎはプロセッサの
番号、ｔは自然数で、１≦ｔ≦乱数発生数、**はべき算
からなる式により乱数列を前記複数のプロセッサのそれ
ぞれに発生させる発生ステップと、Ｘnを発生乱数列格
納部内のｎ番目の領域に格納し、Ｘn+rtを該発生乱数列
格納部内のｎ＋ｒｔ番目の領域に格納することによっ
て、一つの乱数列にする格納ステップとを含むことを特
徴とする。
【００１３】本発明の第２のマルチプロセッサシステム
における乱数発生方法は、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄ
Ｍただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎ
はプロセッサの番号、**はべき算からなる式により算出
した初期値と乱数発生数＝（発生総数／プロセッサ数）
−１ただし、乱数発生数はそれぞれのプロセッサに発生
させる乱数の数、発生総数は全部のプロセッサに発生さ
せる乱数の総数、プロセッサは乱数を発生させるプロセ
ッサ数からなる式により算出した乱数発生数とを設定す
る設定ステップと、該初期値と該乱数発生数とに基づい
て、Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1)ただし、Ｘは整
数、αは初期値、ｒはプロセッサ数、ｎはプロセッサの
番号、ｔは自然数で、１≦ｔ≦乱数発生数、**はべき算
からなる式により乱数列を前記複数のプロセッサのそれ
ぞれに発生させる発生ステップと、前記発生ステップに
より乱数を発生させている最中に乱数を発生させるプロ
セッサを追加する追加ステップと、前記追加したプロセ
ッサを含める前記複数のプロセッサのそれぞれにＸq+n
＝Ｘq＊（α**ｎ）ｍｏｄＭただし、α、Ｘ、およびＭ
は整数で、０≦Ｘ＜Ｍ、ｎはプロセッサの番号、ｑは既
設のすべてのプロセッサが発生した乱数の合計、**はべ
き算からなる式により算出した初期値と乱数発生数＝
｛（発生総数−発生済数）／プロセッサ数｝−１ただ
し、乱数発生数はそれぞれのプロセッサに発生させる乱
数の数、発生総数は全部のプロセッサに発生させる乱数
の総数、発生済数は発生済みの乱数の個数、プロセッサ
は乱数を発生させるプロセッサ数からなる式により算出
した乱数発生数とを再設定する再設定ステップと、該初
期値と該乱数発生数とに基づいてＸn+q+rt＝（α**ｒ）
＊Ｘn+q+r(t-1)ただし、Ｘは整数、αは初期値、ｒはプ
ロセッサ数、ｎはプロセッサの番号、ｑは既設のすべて
のプロセッサが発生した乱数の合計、ｔは自然数で、１
≦ｔ≦乱数発生数、**はべき算からなる式による乱数列
を再度発生させる再発生ステップと、Ｘnを発生乱数列
格納部内のｎ番目の領域に格納し、Ｘn+rtを該発生乱数
列格納部内のｎ＋ｒｔ番目の領域に格納して、さらに、
Ｘq+nを該発生乱数列格納部内のｑ＋ｎ番目の領域に格
納し、Ｘn+q+rtを該発生乱数列格納部内のｎ＋ｑ＋ｒｔ
番目の領域に格納することによって、一つの乱数列にす
る格納ステップとを含むことを特徴とする。
における乱数発生方法は、Ｘn＝Ｘ＊（α**ｎ）ｍｏｄ
Ｍただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎ
はプロセッサの番号、**はべき算からなる式により算出
した初期値と乱数発生数＝（発生総数／プロセッサ数）
−１ただし、乱数発生数はそれぞれのプロセッサに発生
させる乱数の数、発生総数は全部のプロセッサに発生さ
せる乱数の総数、プロセッサは乱数を発生させるプロセ
ッサ数からなる式により算出した乱数発生数とを設定す
る設定ステップと、該初期値と該乱数発生数とに基づい
て、Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1)ただし、Ｘは整
数、αは初期値、ｒはプロセッサ数、ｎはプロセッサの
番号、ｔは自然数で、１≦ｔ≦乱数発生数、**はべき算
からなる式により乱数列を前記複数のプロセッサのそれ
ぞれに発生させる発生ステップと、前記発生ステップに
より乱数を発生させている最中に乱数を発生させるプロ
セッサを追加する追加ステップと、前記追加したプロセ
ッサを含める前記複数のプロセッサのそれぞれにＸq+n
＝Ｘq＊（α**ｎ）ｍｏｄＭただし、α、Ｘ、およびＭ
は整数で、０≦Ｘ＜Ｍ、ｎはプロセッサの番号、ｑは既
設のすべてのプロセッサが発生した乱数の合計、**はべ
き算からなる式により算出した初期値と乱数発生数＝
｛（発生総数−発生済数）／プロセッサ数｝−１ただ
し、乱数発生数はそれぞれのプロセッサに発生させる乱
数の数、発生総数は全部のプロセッサに発生させる乱数
の総数、発生済数は発生済みの乱数の個数、プロセッサ
は乱数を発生させるプロセッサ数からなる式により算出
した乱数発生数とを再設定する再設定ステップと、該初
期値と該乱数発生数とに基づいてＸn+q+rt＝（α**ｒ）
＊Ｘn+q+r(t-1)ただし、Ｘは整数、αは初期値、ｒはプ
ロセッサ数、ｎはプロセッサの番号、ｑは既設のすべて
のプロセッサが発生した乱数の合計、ｔは自然数で、１
≦ｔ≦乱数発生数、**はべき算からなる式による乱数列
を再度発生させる再発生ステップと、Ｘnを発生乱数列
格納部内のｎ番目の領域に格納し、Ｘn+rtを該発生乱数
列格納部内のｎ＋ｒｔ番目の領域に格納して、さらに、
Ｘq+nを該発生乱数列格納部内のｑ＋ｎ番目の領域に格
納し、Ｘn+q+rtを該発生乱数列格納部内のｎ＋ｑ＋ｒｔ
番目の領域に格納することによって、一つの乱数列にす
る格納ステップとを含むことを特徴とする。
【００１４】
【００１５】
【００１６】
【００１７】
【００１８】
【００１９】
【００２０】
【００２１】
【００２２】
【００２３】
【００２４】
【００２５】
【００２６】
【００２７】
【００２８】
【発明の実施の形態】以下、本発明の実施例について、
図を参照しながら詳細に説明する。
図を参照しながら詳細に説明する。
【００２９】図１を参照すると、本発明の第１の実施例
であるマルチプロセッサシステムにおける乱数発生装置
は、乱数を発生するプロセッサ３１、３２、３３、およ
び３４と、各プロセッサにおいて乱数を発生するための
初期値およびパラメータ群を設定するプロセッサ管理部
１と、乱数を発生するためのパラメータ群を各プロセッ
サ毎に格納するプロセッサ管理テーブル２と、各プロセ
ッサによって発生された乱数を格納する発生乱数列格納
部４とから構成されている。
であるマルチプロセッサシステムにおける乱数発生装置
は、乱数を発生するプロセッサ３１、３２、３３、およ
び３４と、各プロセッサにおいて乱数を発生するための
初期値およびパラメータ群を設定するプロセッサ管理部
１と、乱数を発生するためのパラメータ群を各プロセッ
サ毎に格納するプロセッサ管理テーブル２と、各プロセ
ッサによって発生された乱数を格納する発生乱数列格納
部４とから構成されている。
【００３０】以下に、本発明の第１の実施例であるマル
チプロセッサシステムにおける乱数発生装置について、
図１〜図３、および図７〜図１１を参照して詳細に説明
する。
チプロセッサシステムにおける乱数発生装置について、
図１〜図３、および図７〜図１１を参照して詳細に説明
する。
【００３１】本実施例においては、各プロセッサが、初
期値を基にして４００個ずつの乱数を発生させることと
する。
期値を基にして４００個ずつの乱数を発生させることと
する。
【００３２】まず、プロセッサ管理部１は、各プロセッ
サ３１〜３４の初期値を次式に基づいて設定する。
サ３１〜３４の初期値を次式に基づいて設定する。
【００３３】Ｘn ＝ Ｘ＊αnｍｏｄＭ （ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎ
はプロセッサの番号とする。） したがって、 プロセッサ３１は、Ｘ1 ＝ Ｘ＊αｍｏｄＭ プロセッサ３２は、Ｘ2 ＝ Ｘ＊α2ｍｏｄＭ プロセッサ３３は、Ｘ3 ＝ Ｘ＊α3ｍｏｄＭ プロセッサ３４は、Ｘ4 ＝ Ｘ＊α4ｍｏｄＭ が初期値となる。
はプロセッサの番号とする。） したがって、 プロセッサ３１は、Ｘ1 ＝ Ｘ＊αｍｏｄＭ プロセッサ３２は、Ｘ2 ＝ Ｘ＊α2ｍｏｄＭ プロセッサ３３は、Ｘ3 ＝ Ｘ＊α3ｍｏｄＭ プロセッサ３４は、Ｘ4 ＝ Ｘ＊α4ｍｏｄＭ が初期値となる。
【００３４】プロセッサ管理部１は、図２に示すよう
に、これらの初期値を発生乱数列格納部４に格納する。
に、これらの初期値を発生乱数列格納部４に格納する。
【００３５】プロセッサ管理部１は、図７に示すよう
に、これらの初期値をプロセッサ管理テーブル２にも格
納する。
に、これらの初期値をプロセッサ管理テーブル２にも格
納する。
【００３６】次に、プロセッサ管理部１は、各プロセッ
サの乱数発生数を次式に基づいて設定する。
サの乱数発生数を次式に基づいて設定する。
【００３７】 乱数発生数 ＝ （発生総数／プロセッサ数）−１ したがって、 乱数発生数 ＝ （４００／４）−１ ＝ ９９ プロセッサ管理部１は、図８に示すように、この乱数発
生数をプロセッサ管理テーブル２の乱数発生数の欄に格
納する。該計算において、除算に剰余が出た場合には、
プロセッサ番号の小さい順に、乱数発生数に該剰余から
１ずつ加える。
生数をプロセッサ管理テーブル２の乱数発生数の欄に格
納する。該計算において、除算に剰余が出た場合には、
プロセッサ番号の小さい順に、乱数発生数に該剰余から
１ずつ加える。
【００３８】次に、プロセッサ管理部１は、各プロセッ
サが乱数生成に使用する乗算合同法の乗数を次式によっ
て決定し、図９に示すように、プロセッサ管理テーブル
２の乗数の欄に格納する。
サが乱数生成に使用する乗算合同法の乗数を次式によっ
て決定し、図９に示すように、プロセッサ管理テーブル
２の乗数の欄に格納する。
【００３９】乗数 ＝ α＊＊ｒ (＝ αr) （αは初期値の値、＊＊はべき算、ｒはプロセッサ数で
ある。） したがって、本実施例では、乗数 ＝ α4 である。
ある。） したがって、本実施例では、乗数 ＝ α4 である。
【００４０】プロセッサ管理部１は、各プロセッサが発
生する乱数を直列に発生乱数列格納部４へ格納する際の
結果格納位置間隔を、図１０に示すように、プロセッサ
管理テーブル２の結果格納位置間隔の欄に格納する。本
実施例では、プロセッサ数が４であることにより、結果
格納位置間隔も４となる。
生する乱数を直列に発生乱数列格納部４へ格納する際の
結果格納位置間隔を、図１０に示すように、プロセッサ
管理テーブル２の結果格納位置間隔の欄に格納する。本
実施例では、プロセッサ数が４であることにより、結果
格納位置間隔も４となる。
【００４１】プロセッサ管理部１は、プロセッサ管理テ
ーブル２内の各プロセッサの発生数カウンタに０を格納
する。
ーブル２内の各プロセッサの発生数カウンタに０を格納
する。
【００４２】次に、プロセッサ管理部１は、各プロセッ
サに対して乱数発生開始を指示する。
サに対して乱数発生開始を指示する。
【００４３】各プロセッサは、プロセッサ管理部１から
乱数の発生開始を指示されると、図１１の状態のプロセ
ッサ管理テーブル２から、自プロセッサのプロセッサ番
号に対応する乗数、初期値、および乱数発生数の数値を
得て、次式の計算を行って乱数を発生する。さらに、各
プロセッサは、発生した乱数を図３に示すように発生乱
数列格納部４に結果格納間隔にしたがって格納するとと
もに、プロセッサ管理テーブル２内の自プロセッサのプ
ロセッサ番号に対応する発生数カウンタの値に１加算す
る。
乱数の発生開始を指示されると、図１１の状態のプロセ
ッサ管理テーブル２から、自プロセッサのプロセッサ番
号に対応する乗数、初期値、および乱数発生数の数値を
得て、次式の計算を行って乱数を発生する。さらに、各
プロセッサは、発生した乱数を図３に示すように発生乱
数列格納部４に結果格納間隔にしたがって格納するとと
もに、プロセッサ管理テーブル２内の自プロセッサのプ
ロセッサ番号に対応する発生数カウンタの値に１加算す
る。
【００４４】Ｘn+rt ＝ αr＊Ｘn+r(t-1) （ただし、α、Ｘは整数で、ｎはプロセッサの番号とす
る。） この式で、ｒ＝４、ｔは１から始まる自然数とすると、
プロセッサ３１は、乱数Ｘ1+4＝α4＊Ｘ1、Ｘ1+4*2＝α
4＊Ｘ1+4、…、プロセッサ３２は、乱数Ｘ2+4＝α4＊Ｘ
2、Ｘ2+4*2＝α4＊Ｘ2+4、…、プロセッサ３３は、乱数
Ｘ3+4＝α4＊Ｘ3、Ｘ3+4*2＝α4＊Ｘ3+4、…、プロセッ
サ３４は、乱数Ｘ4+4＝α4＊Ｘ4、Ｘ4+4*2＝α4＊
Ｘ4+4、…、を発生する。
る。） この式で、ｒ＝４、ｔは１から始まる自然数とすると、
プロセッサ３１は、乱数Ｘ1+4＝α4＊Ｘ1、Ｘ1+4*2＝α
4＊Ｘ1+4、…、プロセッサ３２は、乱数Ｘ2+4＝α4＊Ｘ
2、Ｘ2+4*2＝α4＊Ｘ2+4、…、プロセッサ３３は、乱数
Ｘ3+4＝α4＊Ｘ3、Ｘ3+4*2＝α4＊Ｘ3+4、…、プロセッ
サ３４は、乱数Ｘ4+4＝α4＊Ｘ4、Ｘ4+4*2＝α4＊
Ｘ4+4、…、を発生する。
【００４５】各プロセッサは、この乱数の計算をプロセ
ッサ管理テーブル２内の乱数発生数が示す回数だけ実行
する。
ッサ管理テーブル２内の乱数発生数が示す回数だけ実行
する。
【００４６】以上により、本発明の第１の実施例である
マルチプロセッサシステムにおける乱数発生装置の処理
が完了する。
マルチプロセッサシステムにおける乱数発生装置の処理
が完了する。
【００４７】本発明の第１の実施例であるマルチプロセ
ッサシステムにおける乱数発生装置の特徴は、プロセッ
サ管理部１が各プロセッサにおける乱数発生のための初
期値およびパラメータ群を設定してプロセッサ管理テー
ブル２に格納し、各プロセッサが該プロセッサ管理テー
ブル２に格納された乱数発生のための初期値およびパラ
メータ群に基づいて、乗算合同法により乱数を発生する
ようにしたことにある。
ッサシステムにおける乱数発生装置の特徴は、プロセッ
サ管理部１が各プロセッサにおける乱数発生のための初
期値およびパラメータ群を設定してプロセッサ管理テー
ブル２に格納し、各プロセッサが該プロセッサ管理テー
ブル２に格納された乱数発生のための初期値およびパラ
メータ群に基づいて、乗算合同法により乱数を発生する
ようにしたことにある。
【００４８】本発明の第１の実施例であるマルチプロセ
ッサシステムにおける乱数発生装置は、このような特徴
を備えたことにより、マルチプロセッサシステムにおい
て、効率良く容易に乱数を発生することができる効果を
有している。
ッサシステムにおける乱数発生装置は、このような特徴
を備えたことにより、マルチプロセッサシステムにおい
て、効率良く容易に乱数を発生することができる効果を
有している。
【００４９】図４を参照すると、本発明の第２の実施例
であるマルチプロセッサシステムにおける乱数発生装置
は、第１の実施例であるマルチプロセッサシステムにお
ける乱数発生装置の構成に加えて、乱数を発生するプロ
セッサ３５が追加されている。
であるマルチプロセッサシステムにおける乱数発生装置
は、第１の実施例であるマルチプロセッサシステムにお
ける乱数発生装置の構成に加えて、乱数を発生するプロ
セッサ３５が追加されている。
【００５０】以下に、本発明の第２の実施例であるマル
チプロセッサシステムにおける乱数発生装置について、
図４〜図６、および図１２〜図１８を参照して詳細に説
明する。
チプロセッサシステムにおける乱数発生装置について、
図４〜図６、および図１２〜図１８を参照して詳細に説
明する。
【００５１】本実施例においては、プロセッサ管理テー
ブル２の内容が図１２に示す状態である環境において、
プロセッサ３５を追加する場合について説明する。
ブル２の内容が図１２に示す状態である環境において、
プロセッサ３５を追加する場合について説明する。
【００５２】プロセッサ管理部１は、プロセッサ管理テ
ーブル２を参照し、発生数カウンタの値が最も大きい、
すなわちその時点で乱数を最も多く発生しているプロセ
ッサの番号と該プロセッサの発生数カウンタとを認識す
る。本実施例においては、プロセッサ番号“３”と発生
数カウンタの値“３１”とが認識される。
ーブル２を参照し、発生数カウンタの値が最も大きい、
すなわちその時点で乱数を最も多く発生しているプロセ
ッサの番号と該プロセッサの発生数カウンタとを認識す
る。本実施例においては、プロセッサ番号“３”と発生
数カウンタの値“３１”とが認識される。
【００５３】プロセッサ管理部１は、認識した発生数カ
ウンタの値に予め設定した遅延数を加えた値（以下、発
生一時停止数とする）を求め、既存のプロセッサ３１〜
３４に対して、発生一時停止数に発生数カウンタの値が
達した時点で乱数の発生処理を停止させる。この遅延数
は、プロセッサ管理部１が乱数発生の中断指示を出して
から実際にプロセッサが乱数発生を中断するまでに発生
する乱数の個数より大きい最小の整数値である。本実施
例においては、該遅延数は“１”であり、発生一時停止
数は“３２”となる。
ウンタの値に予め設定した遅延数を加えた値（以下、発
生一時停止数とする）を求め、既存のプロセッサ３１〜
３４に対して、発生一時停止数に発生数カウンタの値が
達した時点で乱数の発生処理を停止させる。この遅延数
は、プロセッサ管理部１が乱数発生の中断指示を出して
から実際にプロセッサが乱数発生を中断するまでに発生
する乱数の個数より大きい最小の整数値である。本実施
例においては、該遅延数は“１”であり、発生一時停止
数は“３２”となる。
【００５４】次に、プロセッサ管理部１は、プロセッサ
管理テーブル２のすべてのプロセッサの発生カウンタの
値が“３２”になるまで、該プロセッサ管理テーブル２
を監視する。プロセッサ管理部１は、各プロセッサの発
生数カウンタが“３２”になると、該プロセッサに乱数
発生処理を停止させる。この時のプロセッサ管理テーブ
ル２の状態を図１３に示す。
管理テーブル２のすべてのプロセッサの発生カウンタの
値が“３２”になるまで、該プロセッサ管理テーブル２
を監視する。プロセッサ管理部１は、各プロセッサの発
生数カウンタが“３２”になると、該プロセッサに乱数
発生処理を停止させる。この時のプロセッサ管理テーブ
ル２の状態を図１３に示す。
【００５５】プロセッサ管理部１は、この時点でプロセ
ッサ管理テーブル２の内容を変更する。すなわち、この
時点の乱数発生合計数は初期値を加えて、 ４＋（３２×４）＝１３２ であり、プロセッサ数が“５”になってからの各プロセ
ッサの初期値を次式に基づいて設定する。
ッサ管理テーブル２の内容を変更する。すなわち、この
時点の乱数発生合計数は初期値を加えて、 ４＋（３２×４）＝１３２ であり、プロセッサ数が“５”になってからの各プロセ
ッサの初期値を次式に基づいて設定する。
【００５６】Ｘ132+n ＝ Ｘ132＊αnｍｏｄＭ （ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎ
はプロセッサの番号とする。） プロセッサ管理部１は、設定した初期値を図５に示すよ
うに発生乱数列格納部４に格納する。
はプロセッサの番号とする。） プロセッサ管理部１は、設定した初期値を図５に示すよ
うに発生乱数列格納部４に格納する。
【００５７】プロセッサ管理部１は、図１４に示すよう
に、これらの初期値をプロセッサ管理テーブル２にも格
納する。
に、これらの初期値をプロセッサ管理テーブル２にも格
納する。
【００５８】次に、プロセッサ管理部１は、乱数の発生
総数４００と中断時の合計発生数１３２とに基づいて、
各プロセッサのこれからの乱数発生数を次式に基づいて
決定する。
総数４００と中断時の合計発生数１３２とに基づいて、
各プロセッサのこれからの乱数発生数を次式に基づいて
決定する。
【００５９】乱数発生数 ＝ ｛（発生総数−発生済数）
／プロセッサ数｝−１ したがって、 乱数発生数 ＝ ｛（４００−１３２）／５｝−１ ＝ ５
２ 該計算において、除算に剰余が出た場合には、プロセッ
サ番号の小さい順に、乱数発生数に該剰余から１ずつ加
える。
／プロセッサ数｝−１ したがって、 乱数発生数 ＝ ｛（４００−１３２）／５｝−１ ＝ ５
２ 該計算において、除算に剰余が出た場合には、プロセッ
サ番号の小さい順に、乱数発生数に該剰余から１ずつ加
える。
【００６０】プロセッサ管理部１は、図１５に示すよう
に、この乱数発生数をプロセッサ管理テーブル２の乱数
発生数の欄に格納する。
に、この乱数発生数をプロセッサ管理テーブル２の乱数
発生数の欄に格納する。
【００６１】次に、プロセッサ管理部１は、各プロセッ
サが乱数生成に使用する乗算合同法の乗数を次式によっ
て決定し、図１６示すように、プロセッサ管理テーブル
２の乗数の欄に格納する。
サが乱数生成に使用する乗算合同法の乗数を次式によっ
て決定し、図１６示すように、プロセッサ管理テーブル
２の乗数の欄に格納する。
【００６２】乗数 ＝ α＊＊ｒ （αは初期値の値、＊＊はべき算、ｒはプロセッサ数で
ある。） したがって、本実施例では、乗数 ＝ α5 である。
ある。） したがって、本実施例では、乗数 ＝ α5 である。
【００６３】プロセッサ管理部１は、各プロセッサが発
生する乱数を直列にして発生乱数列格納部４へ格納する
際の結果格納位置間隔を、図１７に示すように、プロセ
ッサ管理テーブル２の結果格納位置間隔の欄に格納す
る。本実施例では、プロセッサ数が５であることによ
り、結果格納位置間隔も５となる。
生する乱数を直列にして発生乱数列格納部４へ格納する
際の結果格納位置間隔を、図１７に示すように、プロセ
ッサ管理テーブル２の結果格納位置間隔の欄に格納す
る。本実施例では、プロセッサ数が５であることによ
り、結果格納位置間隔も５となる。
【００６４】プロセッサ管理部１は、図１８に示すプロ
セッサ管理テーブル２内の各プロセッサの発生数カウン
タに０を格納する。
セッサ管理テーブル２内の各プロセッサの発生数カウン
タに０を格納する。
【００６５】次に、プロセッサ管理部１は、各プロセッ
サに対して乱数発生開始を指示する。
サに対して乱数発生開始を指示する。
【００６６】各プロセッサは、プロセッサ管理部１から
乱数の発生開始を指示されると、図１８の状態のプロセ
ッサ管理テーブル２から、自プロセッサのプロセッサ番
号に対応する乗数、初期値、および乱数発生数の数値に
したがって、次式の計算を行って乱数を発生する。さら
に、各プロセッサは、該乱数を図６に示すように発生乱
数列格納部４に結果格納間隔にしたがって格納するとと
もに、プロセッサ管理テーブル２内の自プロセッサのプ
ロセッサ番号に対応する発生数カウンタの値に１加算す
る。
乱数の発生開始を指示されると、図１８の状態のプロセ
ッサ管理テーブル２から、自プロセッサのプロセッサ番
号に対応する乗数、初期値、および乱数発生数の数値に
したがって、次式の計算を行って乱数を発生する。さら
に、各プロセッサは、該乱数を図６に示すように発生乱
数列格納部４に結果格納間隔にしたがって格納するとと
もに、プロセッサ管理テーブル２内の自プロセッサのプ
ロセッサ番号に対応する発生数カウンタの値に１加算す
る。
【００６７】Ｘn+q+rt ＝ αr＊Ｘn+q+r(t-1) （ただし、α、Ｘは整数で、ｎはプロセッサの番号とす
る。） この式で、ｒ＝５、ｔは１から始まる自然数、ｑは１３
２であるので、各プロセッサは、乱数Ｘn+132+5＝α5＊
Ｘn+132、Ｘn+132+5*2＝α5＊Ｘn+132+5、…、を発生す
る。
る。） この式で、ｒ＝５、ｔは１から始まる自然数、ｑは１３
２であるので、各プロセッサは、乱数Ｘn+132+5＝α5＊
Ｘn+132、Ｘn+132+5*2＝α5＊Ｘn+132+5、…、を発生す
る。
【００６８】各プロセッサは、この乱数の計算をプロセ
ッサ管理テーブル２内の乱数発生数が示す回数だけ実行
する。
ッサ管理テーブル２内の乱数発生数が示す回数だけ実行
する。
【００６９】以上により、本発明の第２の実施例である
マルチプロセッサシステムにおける乱数発生装置の処理
が完了する。
マルチプロセッサシステムにおける乱数発生装置の処理
が完了する。
【００７０】本発明の第２の実施例であるマルチプロセ
ッサシステムにおける乱数発生装置の特徴は、乱数の発
生途中に新たなプロセッサを追加するために、プロセッ
サ管理部１がプロセッサ管理テーブル２に格納されてい
るパラメータにおいて、既設のプロセッサの乱数の発生
数を同一にし、その時点までに全ての既設のプロセッサ
で発生した乱数の合計数を基とした初期値と以降の乱数
の発生数をプロセッサ追加後のすべてのプロセッサに設
定し、各プロセッサが該プロセッサ管理テーブル２に格
納された乱数発生のための初期値およびパラメータ群に
基づいて、乗算合同法により乱数を発生するようにした
ことにある。
ッサシステムにおける乱数発生装置の特徴は、乱数の発
生途中に新たなプロセッサを追加するために、プロセッ
サ管理部１がプロセッサ管理テーブル２に格納されてい
るパラメータにおいて、既設のプロセッサの乱数の発生
数を同一にし、その時点までに全ての既設のプロセッサ
で発生した乱数の合計数を基とした初期値と以降の乱数
の発生数をプロセッサ追加後のすべてのプロセッサに設
定し、各プロセッサが該プロセッサ管理テーブル２に格
納された乱数発生のための初期値およびパラメータ群に
基づいて、乗算合同法により乱数を発生するようにした
ことにある。
【００７１】本発明の第２の実施例であるマルチプロセ
ッサシステムにおける乱数発生装置は、このような特徴
を備えたことにより、マルチプロセッサシステムにおい
て、乱数発生中に乱数を発生させるプロセッサの数を変
更することができる効果を有している。
ッサシステムにおける乱数発生装置は、このような特徴
を備えたことにより、マルチプロセッサシステムにおい
て、乱数発生中に乱数を発生させるプロセッサの数を変
更することができる効果を有している。
【００７２】
【発明の効果】以上説明したように、本発明は、マルチ
プロセッサシステムにおいて、乱数発生中に乱数を発生
させるプロセッサの数を変更することができるので、マ
ルチプロセッサシステムにおいて実行されるさまざまな
プログラム群と共存して効率良く乱数を発生することが
できる効果を有している。
プロセッサシステムにおいて、乱数発生中に乱数を発生
させるプロセッサの数を変更することができるので、マ
ルチプロセッサシステムにおいて実行されるさまざまな
プログラム群と共存して効率良く乱数を発生することが
できる効果を有している。
【図１】本発明の第１の実施例におけるブロック図であ
る。
る。
【図２】本発明の第１の実施例において、発生乱数列格
納部４に初期値が格納された状態を示す図である。
納部４に初期値が格納された状態を示す図である。
【図３】本発明の第１の実施例において、各プロセッサ
が発生した乱数が発生乱数列格納部４に格納される状態
を示す図である。
が発生した乱数が発生乱数列格納部４に格納される状態
を示す図である。
【図４】本発明の第２の実施例におけるブロック図であ
る。
る。
【図５】本発明の第２の実施例において、発生乱数列格
納部４に初期値が格納された状態を示す図である。
納部４に初期値が格納された状態を示す図である。
【図６】本発明の第２の実施例において、各プロセッサ
が発生した乱数が発生乱数列格納部４に格納される状態
を示す図である。
が発生した乱数が発生乱数列格納部４に格納される状態
を示す図である。
【図７】本発明の第１の実施例におけるプロセッサ管理
テーブル２の内容を示す図である。
テーブル２の内容を示す図である。
【図８】本発明の第１の実施例におけるプロセッサ管理
テーブル２の内容を示す図である。
テーブル２の内容を示す図である。
【図９】本発明の第１の実施例におけるプロセッサ管理
テーブル２の内容を示す図である。
テーブル２の内容を示す図である。
【図１０】本発明の第１の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１１】本発明の第１の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１２】本発明の第２の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１３】本発明の第２の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１４】本発明の第２の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１５】本発明の第２の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１６】本発明の第２の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１７】本発明の第２の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
【図１８】本発明の第２の実施例におけるプロセッサ管
理テーブル２の内容を示す図である。
理テーブル２の内容を示す図である。
１ プロセッサ管理部 ２ プロセッサ管理テーブル ４ 発生乱数列格納部 ３１ プロセッサ ３２ プロセッサ ３３ プロセッサ ３４ プロセッサ ３５ プロセッサ
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平４−242827（ＪＰ，Ａ) 特開 平６−202856（ＪＰ，Ａ) 特開 平４−370827（ＪＰ，Ａ) 特開 昭60−136876（ＪＰ，Ａ) Ｋｎｕｔｈ著，渋谷政昭訳「準数値算 法／乱数」、第９−11頁、株式会社サイ エンス社発行
Claims (4)
- 【請求項１】 複数のプロセッサのそれぞれにより発生
された乱数列を組み合わせた乱数列が一つのプロセッサ
により発生された乱数列と同一になるように、 Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭ ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎは
プロセッサの番号、**はべき算からなる式により算出し
た初期値と乱数発生数＝（発生総数／プロセッサ数）−
１ただし、乱数発生数はそれぞれのプロセッサに発生さ
せる乱数の個数、発生総数は全部のプロセッサに発生さ
せる乱数の総数、プロセッサは乱数を発生させるプロセ
ッサ数からなる式により算出した乱数発生数とを設定し
て、該初期値と該乱数発生数とに基づいて Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1) ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式により乱数列を前記複数のプ
ロセッサのそれぞれに発生させ、Ｘnを発生乱数列格納
部内のｎ番目の領域に格納し、Ｘn+rtを該発生乱数列格
納部内のｎ＋ｒｔ番目の領域に格納することによって、
一つの乱数列にするプロセッサ管理部を備えたことを特
徴とするマルチプロセッサシステムにおける乱数発生装
置。 - 【請求項２】 複数のプロセッサのそれぞれにより発生
された乱数列を組み合わせた乱数列が１つのプロセッサ
により発生された乱数列と同一になるように、 Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭ ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎは
プロセッサの番号、**はべき算からなる式により算出し
た初期値と 乱数発生数＝（発生総数／プロセッサ数）−１ ただし、乱数発生数はそれぞれのプロセッサに発生させ
る乱数の数、発生総数は全部のプロセッサに発生させる
乱数の総数、プロセッサは乱数を発生させるプロセッサ
数からなる式により算出した乱数発生数とを設定して、
該初期値と該乱数発生数とに基づいて、 Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1) ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式により乱数列を前記複数のプ
ロセッサのそれぞれに発生させている最中に乱数を発生
させるプロセッサを追加し、追加したプロセッサを含め
る前記複数のプロセッサのそれぞれに Ｘq+n＝Ｘq＊（α**ｎ）ｍｏｄＭ ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎは
プロセッサの番号、ｑは既設のすべてのプロセッサが発
生した乱数の合計、**はべき算からなる式により算出し
た初期値と 乱数発生数＝｛（発生総数−発生済数）／プロセッサ
数｝−１ ただし、乱数発生数はそれぞれのプロセッサに発生させ
る乱数の数、発生総数は全部のプロセッサに発生させる
乱数の総数、発生済数は発生済みの乱数の個数、プロセ
ッサは乱数を発生させるプロセッサ数からなる式により
算出した乱数発生数とを再設定して、該初期値と該乱数
発生数とに基づいて Ｘn+q+rt＝（α**ｒ）＊Ｘn+q+r(t-1) ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｑは既設のすべてのプロセッサが
発生した乱数の合計、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式による乱数列を発生させ、Ｘ
nを発生乱数列格納部内のｎ番目の領域に格納し、Ｘn+r
tを該発生乱数列格納部内のｎ＋ｒｔ番目の領域に格納
して、さらに、Ｘq+nを該発生乱数列格納部内のｑ＋ｎ
番目の領域に格納し、Ｘn+q+rtを該発生乱数列格納部内
のｎ＋ｑ＋ｒｔ番目の領域に格納することによって、一
つの乱数列にするプロセッサ管理部を備えたことを特徴
とするマルチプロセッサシステムにおける乱数発生装
置。 - 【請求項３】 Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭ ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎは
プロセッサの番号、**はべき算からなる式により算出し
た初期値と 乱数発生数＝（発生総数／プロセッサ数）−１ ただし、乱数発生数はそれぞれのプロセッサに発生させ
る乱数の個数、発生総数は全部のプロセッサに発生させ
る乱数の総数、プロセッサは乱数を発生させるプロセッ
サ数からなる式により算出した乱数発生数とを設定する
設定ステップと、 該初期値と該乱数発生数とに基づいて Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1) ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式により乱数列を前記複数のプ
ロセッサのそれぞれに発生させる発生ステップと、 Ｘnを発生乱数列格納部内のｎ番目の領域に格納し、Ｘn
+rtを該発生乱数列格納部内のｎ＋ｒｔ番目の領域に格
納することによって、一つの乱数列にする格納ステップ
とを含むことを特徴とするマルチプロセッサシステムに
おける乱数発生方法。 - 【請求項４】 Ｘn＝Ｘ＊（α**ｎ）ｍｏｄＭ ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎは
プロセッサの番号、**はべき算からなる式により算出し
た初期値と 乱数発生数＝（発生総数／プロセッサ数）−１ ただし、乱数発生数はそれぞれのプロセッサに発生させ
る乱数の数、発生総数は全部のプロセッサに発生させる
乱数の総数、プロセッサは乱数を発生させるプロセッサ
数からなる式により算出した乱数発生数とを設定する設
定ステップと、 該初期値と該乱数発生数とに基づいて、 Ｘn+rt＝（α**ｒ）＊Ｘn+r(t-1) ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式により乱数列を前記複数のプ
ロセッサのそれぞれに発生させる発生ステップと、 前記発生ステップにより乱数を発生させている最中に乱
数を発生させるプロセッサを追加する追加ステップと、 前記追加したプロセッサを含める前記複数のプロセッサ
のそれぞれに Ｘq+n＝Ｘq＊（α**ｎ）ｍｏｄＭ ただし、α、Ｘ、およびＭは整数で、０≦Ｘ＜Ｍ、ｎは
プロセッサの番号、ｑは既設のすべてのプロセッサが発
生した乱数の合計、**はべき算からなる式により算出し
た初期値と 乱数発生数＝｛（発生総数−発生済数）／プロセッサ
数｝−１ ただし、乱数発生数はそれぞれのプロセッサに発生させ
る乱数の数、発生総数は全部のプロセッサに発生させる
乱数の総数、発生済数は発生済みの乱数の個数、プロセ
ッサは乱数を発生させるプロセッサ数からなる式により
算出した乱数発生数とを再設定する再設定ステップと、 該初期値と該乱数発生数とに基づいて Ｘn+q+rt＝（α**ｒ）＊Ｘn+q+r(t-1) ただし、Ｘは整数、αは初期値、ｒはプロセッサ数、ｎ
はプロセッサの番号、ｑは既設のすべてのプロセッサが
発生した乱数の合計、ｔは自然数で、１≦ｔ≦乱数発生
数、**はべき算からなる式による乱数列を再度発生させ
る再発生ステップと、 Ｘnを発生乱数列格納部内のｎ番目の領域に格納し、Ｘn
+rtを該発生乱数列格納部内のｎ＋ｒｔ番目の領域に格
納して、さらに、Ｘq+nを該発生乱数列格納部内のｑ＋
ｎ番目の領域に格納し、Ｘn+q+rtを該発生乱数列格納部
内のｎ＋ｑ＋ｒｔ番目の領域に格納することによって、
一つの乱数列にする格納ステップとを含むことを特徴と
するマルチプロセッサシステムにおける乱数発生方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP08083416A JP3123920B2 (ja) | 1995-04-11 | 1996-04-05 | マルチプロセッサシステムにおける乱数発生装置および乱数発生方法 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP8552995 | 1995-04-11 | ||
| JP7-85529 | 1995-04-11 | ||
| JP08083416A JP3123920B2 (ja) | 1995-04-11 | 1996-04-05 | マルチプロセッサシステムにおける乱数発生装置および乱数発生方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08339294A JPH08339294A (ja) | 1996-12-24 |
| JP3123920B2 true JP3123920B2 (ja) | 2001-01-15 |
Family
ID=26424442
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP08083416A Expired - Fee Related JP3123920B2 (ja) | 1995-04-11 | 1996-04-05 | マルチプロセッサシステムにおける乱数発生装置および乱数発生方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3123920B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3539870B2 (ja) * | 1998-07-14 | 2004-07-07 | 株式会社日立製作所 | 乱数生成システム及び乱数生成方法 |
| JP4950325B2 (ja) * | 2010-07-02 | 2012-06-13 | みずほ第一フィナンシャルテクノロジー株式会社 | モンテカルロ法の効率的な並列処理手法 |
-
1996
- 1996-04-05 JP JP08083416A patent/JP3123920B2/ja not_active Expired - Fee Related
Non-Patent Citations (1)
| Title |
|---|
| Ｋｎｕｔｈ著，渋谷政昭訳「準数値算法／乱数」、第９−１１頁、株式会社サイエンス社発行 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08339294A (ja) | 1996-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Brent | Factorization of the tenth Fermat number | |
| US5793657A (en) | Random number generating apparatus and random number generating method in a multiprocessor system | |
| JP2003330734A (ja) | 重みを使用してアプリケーションにシステム資源を割当てる方法およびシステム | |
| JPS59160267A (ja) | ベクトル処理装置 | |
| CN108762897A (zh) | 多任务管理方法及智能手表 | |
| JP2004038291A (ja) | クロック周波数の制御方法および電子機器 | |
| Perarnau et al. | Victim selection and distributed work stealing performance: A case study | |
| CN112000484A (zh) | 栈帧地址随机化方法及相关设备 | |
| JP3123920B2 (ja) | マルチプロセッサシステムにおける乱数発生装置および乱数発生方法 | |
| CN110232289A (zh) | 椭圆曲线密码的高速倍点运算方法 | |
| JP6316593B2 (ja) | 半導体装置 | |
| CN109375895B (zh) | 多项式乘法的加速方法及装置、ntru加解密加速方法及装置 | |
| Cucu et al. | Non-preemptive multiprocessor scheduling for strict periodic systems with precedence constraints | |
| JP2007193744A (ja) | 情報処理装置、プログラム、およびスケジューリング方法 | |
| JP2003122251A (ja) | 暗号情報生成方法と暗号情報生成装置、暗号情報生成プログラム及び記録媒体 | |
| JPH11184717A (ja) | スケジューリング解析装置及び方法並びにスケジューリング解析用ソフトウェアを記録した記録媒体 | |
| JP3686261B2 (ja) | 計算機、プログラム変換装置及びプログラム記録媒体 | |
| JP3223530B2 (ja) | データ転送命令生成処理方法 | |
| JPH09167114A (ja) | マルチスレッドタスクのスタック割り当て方法 | |
| JPH11345129A (ja) | インターバルタイマ回路 | |
| Kaviani et al. | On scheduling in multiprocessor systems using fuzzy logic | |
| JP2021149793A (ja) | 情報処理装置、情報処理プログラムおよび情報処理方法 | |
| JPH03163630A (ja) | プロセッサのスタック制御方法 | |
| JPS61182135A (ja) | 処理選択方法 | |
| JPS633336A (ja) | システム生成方式 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071027 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081027 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |