家庭用ゲーム機プレイステーション2による大規模並列流体解析システムの構築

使用家用游戏机 PlayStation 2 构建大型并行流体分析系统

• 批准号: 12555021
• 负责人: 青木 尊之
• 金额: $ 7.81万
• 依托单位: Tokyo Institute of Technology
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Scientific Research (B)
• 财政年份: 2000
• 资助国家: 日本
• 起止时间: 2000 至 2002
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/en/grant/KAKENHI-PROJECT-12555021/
• 关键词: プレイステーション2; Linux; 並列計算; 流体解析; Power PC; G4 Cube; Emotion Engine; Velocity Engine

昨年度導入した20台のプレイステーション2 Linux kitで、並列の数値流体シミュレーションを行った。非圧縮性流体の場合、大量のネットワーク通信を伴うポアソン方程式の解法が最も大きな課題となる。並列化が容易なSOR法で2次中心差分法による離散化式と、局所補間微分オペレータ法による4次精度の離散化式を解いた。Linux kitに付属するGNUのコンパイラは、プレイステーション2のCPUであるEmotion Engine(EE)の高速化アーキテクチャを利用する最適化が行えず、インライン・アセンブラにより直接EEのベクトル・レジスタVU0を利用することにした。ベクトル・レジスタVU0を4並列で効率よく実行する必要があるが、メモリーからデータをロードするのに時間がかかり、その間はVU0が休んでしまう結果となった。これを回避するために、VU0の演算スケジュールを工夫し、データのロード・ストアを含めて450MFlopsの計算速度を達成した。プレイステーション2が高速のRAMBUSメモリーを用いているために、CPUキャッシュに入りきらないデータに対するアクセスに対しても余り速度が低下しないことが明らかになった。EEのもう一つの問題点として指摘されている「VU0が単精度演算しか行えないこと」を解決する必要がある。倍精度計算に対しては浮動小数点演算プロセッサが働かなくなり、計算速度は100分の一以下に低下してしまう。単精度計算で直接従属変数からポアソン方程式を計算すると桁落ちが生じてしまうので、残差から修正方程式を導出し修正量を求めて従属変数に加算・減算して修正する方法を開発した。修正方程式を解く際に桁落ちなどの誤差が入るが、反復ごとに修正値が小さな値に移行し、従属変数の精度範囲をまかなえることが明らかになった。非圧縮性流体の移流部分や圧縮性流体に対しては局所補間微分オペレータ法が少ないメモリーに対して高精度な計算を行うことができ、プレイステーション2に適していることが分かった。また、ルンゲクッタ時間積分を行うことにより安定な計算ができることが明らかになった。

Yesterday's annual import し た 20 の プ レ イ ス テ ー シ ョ ン 2 kit で, tied for Linux の the numerical fluid シ ミ ュ レ ー シ ョ ン を line っ た. In non-compressible fluid <s:1> situations, a large amount of <s:1> ネットワ and <s:1> communication を accompanied by うポアソ, the solution method of <s:1> equations が is the most <s:1> major <s:1> な topic となる. Parallel で が easy な SOR method two center difference method に よ る discretization type differential between と, bureau オ ペ レ ー タ method に よ る four precision の discretization type を solution い た. Linux kit に pay genus す る GNU の コ ン パ イ ラ は, プ レ イ ス テ ー シ ョ ン 2 の CPU で あ る Emotion High speed Engine (EE) の ア ー キ テ ク チ ャ を using す る line optimization が え ず, イ ン ラ イ ン · ア セ ン ブ ラ に よ り EE directly の ベ ク ト ル · レ ジ ス タ VU0 を using す る こ と に し た. ベ ク ト ル · レ ジ ス タ VU0 を 4 parallel で sharper rate よ く line be す る necessary が あ る が, メ モ リ ー か ら デ ー タ を ロ ー ド す る の に time が か か り, そ の between は VU0 が Hugh ん で し ま う results と な っ た. こ れ を avoid す る た め に, VU0 の calculus ス ケ ジ ュ ー ル を time し, デ ー タ の ロ ー ド · ス ト ア を containing め て 450 mflops を の calculation speed reached し た. プ レ イ ス テ ー シ ョ ン 2 が high-speed の RAMBUS メ モ リ ー を with い て い る た め に, CPU キ ャ ッ シ ュ に into り き ら な い デ ー タ に す seaborne る ア ク セ ス に し seaborne て も yu り low speed が し な い こ と が Ming ら か に な っ た. EE の も う a つ の problem point と し て blame さ れ て い る "VU0 が 単 precision calculus し か line え な い こ と" を solve す る necessary が あ る. Double precision calculation に し seaborne て は floating point calculations プ ロ セ ッ サ が 働 か な く な は り, calculation speed below 100 points の a low に し て し ま う. 単 precision calculation of で 従 directly belong to - か ら ポ ア ソ を ン equations calculating す る と girder falling ち が raw じ て し ま う の で, residual か ら correction equation is derived を し correction を o め て 従 genera - several に addition, reduced calculate し て correction す る method を open 発 し た. Correction equation を solution く interstate に girder falling ち な ど が の error into the る が, repeated ご と に fixed numerical が small さ な numerical に transitional し, 従 - several の precision van 囲 を ま か な え る こ と が Ming ら か に な っ た. Non 圧 shrinkage fluid の advection part や 圧 shrinkage fluid に し seaborne て は bureau to fill between differential オ ペ レ ー タ method が less な い メ モ リ ー に し seaborne て high-precision な count を う こ と が で き, プ レ イ ス テ ー シ ョ ン 2 に optimum し て い る こ と が points か っ た. ま た, ル ン ゲ ク ッ タ time line integral を う こ と に よ り settle な computing が で き る こ と が Ming ら か に な っ た.

项目成果

T.Aoki: "3D Simulation for Falling Papers"Comput. Phys. Comm.. Vol.142, No.1-3. 326-329 (2002)

T.Aoki：“纸张掉落的 3D 模拟”计算。

保坂光紀,塙雅典,高原幹夫: "非線形伝送の原理を応用した高密度WDMシステムにおける長距離伝送の検討"電子情報通信学会論文誌B. J83-B,No.12. 1779-1781 (2000)

Mitsunori Hosaka、Masanori Hanawa、Mikio Takahara：“应用非线性传输原理的高密度 WDM 系统中的长距离传输研究” IEICE Transactions B. J83-B，No. 12. 1779-1781 (2000)

今井陽介, 青木尊之, 小原徹也: "IDO法によるLMR(Local Mesh Refinement)を用いた血流シミュレーション"日本機械学会 第15回 バイオエンジニアリング講演集. 15. 167-168 (2003)

Yosuke Imai、Takayuki Aoki、Tetsuya Ohara：“使用 IDO 方法使用 LMR（局部网格细化）进行血流模拟”第 15 届日本机械工程师学会生物工程会议记录（2003 年）。

Takayuki Aoki, Tetsuya Kobara, Mamoru Tanahashi: "Interpolated Differential Operator (IDO) Scheme and Application Test to Turbulent Flow"5th China-Japan Workshop on Turbulent Flows, Tianjin, China. 5. 1-2 (2002)

Takayuki Aoki、Tetsuya Kobara、Mamoru Tanahashi：“插值微分算子（IDO）方案及其在湍流中的应用测试”第五届中日湍流研讨会，中国天津。

K.Kato,T.Aoki,T.Sekine: "A Numerical Simulation of State of Compressed Diamond"The 8^<th>NIRIM International Symposium on Advanced Materials(3月4日〜8日つくば). (発表予定). (2001)

K.Kato、T.Aoki、T.Sekine：“压缩金刚石状态的数值模拟”第 8 届 NIRIM 国际先进材料研讨会（3 月 4 日至 8 日筑波）（预定演讲）。