㸳ḟඖࣈࣛࢯࣇࢥ࣮ࢻ࡟ࡼࡿ↓⾪✺ࣉࣛࢬ࣐ࢩ࣑࣮ࣗࣞࢩࣙࣥࡢᛶ⬟ホ౯ᱵ⏣㝯⾜ྡྂᒇ኱Ꮫኴ㝧ᆅ⌫⎔ቃ◊✲ᡤ῝ἑᆂ୍㑻஑ᕞ኱Ꮫ᝟ሗᇶ┙◊✲㛤Ⓨࢭࣥࢱ࣮㸯㸬ࡣࡌࡵ࡟ᡃࠎࡀఫࡴᏱᐂࡢ 99.99%௨ୖࡢయ✚ࡣࣉࣛࢬ࣐࡜࿧ࡤࢀࡿ㟁㞳Ẽయ࡛༨ࡵࡽࢀ࡚࠸ࡿࠋᏱᐂ✵㛫࡟Ꮡᅾࡍࡿࣉࣛࢬ࣐ࡢ኱㒊ศࡣᐦᗘࡀ㠀ᖖ࡟ᑠࡉࡃ↓⾪✺≧ែ࡟࠶ࡾ㸪Ᏹᐂࣉࣛࢬ࣐㸦↓⾪✺ࣉࣛࢬ࣐㸧ࢆ⌮ゎࡍࡿࡇ࡜ࡣ㸪Ᏹᐂࡢᮏ㉁ⓗ࡞⌮ゎ࡟ࡘ࡞ࡀࡿࠋᆅ⌫࿘㎶ࡢᏱᐂ⎔ቃࡣ㸪ኴ㝧࠿ࡽᨺฟࡉࢀࡓ㧗㏿ࡢࣉࣛࢬ࣐ὶ࡛࠶ࡿ͆ኴ㝧㢼͇ཬࡧ㸪ኴ㝧㢼ࡀ㐠ࡪᝨᫍ㛫✵㛫☢ሙ㸦ኴ㝧ࡢᅛ᭷☢ሙ㸧࡜㸪ᆅ⌫ࡢᅛ᭷☢ሙ࡜ࡢ┦஫స⏝࡟ࡼࡗ࡚」㞧࡞☢Ẽᅪᵓ㐀ࢆᙧᡂࡋ࡚࠸ࡿࠋࣉࣛࢬ࣐ᨺฟ⌧㇟ࢆࡣࡌࡵ࡜ࡍࡿኴ㝧ࡢᵝࠎ࡞ኚື࡟ࡼࡾ㸪ேᕤ⾨ᫍࡢᨾ㞀ࡸᏱᐂ㣕⾜ኈࡢ⿕᭚࡟⧅ࡀࡿᆅ⌫☢Ẽᅪࡢ⎔ቃኚືࡀᘬࡁ㉳ࡇࡉࢀ㸪ࡇࢀࢆᏱᐂኳẼ࡜࿧ࡪࠋ㏆ᖺࡢᅜ㝿Ᏹᐂࢫࢸ࣮ࢩ࡛ࣙࣥࡢάືࡸேᕤ⾨ᫍᡴࡕୖࡆ࡞࡝㸪᪥ᮏ࡟࠾࠸࡚ࡶᏱᐂ฼⏝ࡀ⌧ᐇⓗ࡟࡞ࡗ࡚ࡁ࡚࠾ࡾ㸪ᏱᐂኳẼࡢணሗ࣭ண ࡟⧅ࡀࡿᏱᐂࣉࣛࢬ࣐◊✲ࡣ㠀ᖖ࡟㔜せ࡛࠶ࡿࠋᆅ⌫☢Ẽᅪෆ࡟ࡣ㸪ࣉࣛࢬ࣐ࡢᐦᗘࡸ ᗘ࡞࡝ࡢ≀⌮ࣃ࣓࣮ࣛࢱࡀ␗࡞ࡿᵝࠎ࡞㡿ᇦࡀ⏕ࡌࡿࠋࡑࡢ㡿ᇦ㛫ࡢቃ⏺ᒙ࡛⌧ࢀࡿ୙Ᏻᐃᛶ㸦ᖹ⾮≧ែࡢ◚ࢀ㸧ࡣ㸪☢Ẽᅪࡢኚື࡟኱ࡁ࡞ᙳ㡪ࢆ୚࠼࡚࠸ࡿ࡜⪃࠼ࡽࢀ࡚࠸ࡿࠋࢢ࣮ࣟࣂࣝ☢Ẽᅪᵓ㐀࡟ᑐࡋ࡚㸪ቃ⏺ᒙ୙Ᏻᐃᛶࡣ୰㛫㸦࣓ࢰ㸧ࢫࢣ࣮ࣝ⌧㇟࡜࿧ࡤࢀࡿࠋࡇࢀࡽࡢࢢ࣮ࣟࣂࣝཬࡧ୰㛫ࢫࢣ࣮ࣝࡢ⌧㇟ࡣ㸪ࣈࣛࢯࣇ㸦↓⾪✺࣎ࣝࢶ࣐ࣥ㸧᪉⛬ᘧࡢ 0 ḟ࣭1 ḟ࣭2 ḟࡢ࣮࣓ࣔࣥࢺࢆྲྀࡿࡇ࡜࡟ࡼࡗ࡚ồࡵࡽࢀࡿ☢ẼὶయຊᏛ㸦MHD㸧᪉⛬ᘧ࡟ࡼࡗ࡚グ㏙ࡉࢀࡿࠋࡋ࠿ࡋ㸪㏆ᖺࡢ⛉Ꮫ⾨ᫍ࡟ࡼࡿ㧗⢭ᗘ࡞ࠕࡑࡢሙࠖほ࡛ࡣ㸪୰㛫ࢫࢣ࣮ࣝࡢ୙Ᏻᐃᛶࡣ㸪MHD ᪉⛬ᘧ࡛グ㏙࡛ࡁࡿ࣐ࢡࣟ≀⌮㐣⛬࡜⢏Ꮚࡢ㐠ືㄽ᪉⛬ᘧ࡟ࡼࡗ࡚グ㏙࡛ࡁࡿ࣑ࢡࣟ≀⌮㐣⛬ࡀ⤖ྜࡋ࡚࠸ࡿࡇ࡜ࢆ♧၀ࡋ࡚࠸ࡿࠋࡇࢀࡽࡢ࣐ࣝࢳࢫࢣ࣮ࣝࡢ☢Ẽᅪኚື࡛࠶ࡿᏱᐂኳẼࢆ┿࡟⌮ゎࡍࡿࡓࡵ࡟ࡣ㸪඲࡚ࡢࢫࢣ࣮ࣝࢆࢩ࣮࣒ࣞࢫ࡟ᢅ࠼ࡿ㐠ືㄽ᪉⛬ᘧ㸦➨୍ཎ⌮㸧࡟ࡼࡿࢩ࣑࣮ࣗࣞࢩࣙࣥࡀᮏ㉁ⓗ࡛࠶ࡿࠋࣉࣛࢬ࣐ࡢ㐠ືㄽࢩ࣑࣮ࣗࣞࢩࣙࣥ࡟ࡣ 2 ࡘࡢᡭἲࡀ࠶ࡿࠋ1 ࡘࡣ㸪ࣉࣛࢬ࣐⢏Ꮚ࡛࠶ࡿ࢖࢜ࣥࡸ㟁Ꮚ࡞࡝ࡢಶࠎࡢⲴ㟁⢏Ꮚࡢ㐠ືࢆ㸪ࢽ࣮ࣗࢺࣥ㸫࣮ࣟࣞࣥࢶ᪉⛬ᘧ࡟ࡼࡾゎࡁ㐍ࡵࡿ
PIC㸦Particle-In-Cell㸧ἲ࡛࠶ࡿࠋ᱁Ꮚ㸦Cell㸧ୖ࡟ᐃ⩏ࡉࢀࡓ㟁☢ሙ୰ࢆ⢏Ꮚࡀືࡁࡲࢃࡿࡇ࡜࠿ࡽ㸪ࡇࡢࡼ࠺࡟࿧ࡤࢀ࡚࠸ࡿࠋᏱᐂ✵㛫࡟Ꮡᅾࡍࡿ⭾኱࡞ᩘࡢⲴ㟁⢏Ꮚࢆ᭷㝈ࡢィ⟬ᶵ㈨※࡛ᢅ࠺ࡇ࡜ࡣ୙ྍ⬟࡛࠶ࡿࡓࡵ㸪࠶ࡿ⛬ᗘࡲ࡜ࡲࡗࡓᩘࡢⲴ㟁⢏Ꮚࡢ㞟ᅋࢆ 1 ࡘࡢ͆㉸͇⢏Ꮚ࡜ࡋ࡚ᢅ࠺ࠋPIC ἲࡣࡑࡢᩘ್ゎἲࡢ᏶ᡂᗘࡀ㧗ࡃ㸪ࣉࣛࢬ࣐⛉Ꮫศ㔝࡛ࡣᗈࡃ⏝࠸ࡽࢀ࡚࠸
ࡿࠋࡋ࠿ࡋ㸪㉸⢏Ꮚ࡜ࡋ࡚ᢅ࠺ࡇ࡜࡟ࡼࡾ⇕㞧㡢ࡀ኱ࡁࡃ࡞ࡿࡇ࡜㸪㟁Ⲵᐦᗘࡸ㟁ὶᐦᗘ࡞࡝
ࡢⲴ㟁⢏Ꮚࡢ㐠ື࡟㉳ᅉࡍࡿሙࡢ㔞ࢆ᱁Ꮚୖ࡟๭ࡾ᣺ࡿ㝿࡟⏕ࡌࡿ㧗Ἴᩘ࣮ࣔࢻࡀᩘ್ㄗᕪ࡜
ࡋ࡚⵳✚ࡍࡿࡇ࡜㸪ࡉࡽ࡟୪ิ໬ࡢ㝿࡟㈇Ⲵࡢࣂࣛࣥࢫ㸦ྛࣉࣟࢭࢫෆࡢ⢏Ꮚᩘࡢᆒ୍ᛶ㸧ࢆ

スーパーコンピューティングニュース- 47 -

Vol. 15, No. 12013ಖࡘࡓࡵ࡟≉Ṧ࡞ࢹ࣮ࢱศ๭ࡀᚲせ࡟࡞ࡿࡇ࡜࡞࡝ࡢḞⅬࡀ࠶ࡿࠋ୍᪉ࡶ࠺ 1 ࡘࡢ㐠ືㄽⓗᡭἲ࡛࠶ࡿࣈࣛࢯࣇἲࡣ㸪఩⨨㸫㏿ᗘ఩┦✵㛫࡟ᐃ⩏ࡉࢀࡓࣉࣛࢬ࣐⢏Ꮚࡢศᕸ㛵ᩘࡢ᫬㛫Ⓨᒎࢆ㸪ࣈࣛࢯࣇ᪉⛬ᘧ࡟ࡼࡾ┤᥋ゎࡁ㐍ࡵࡿ᪉ἲ࡛࠶ࡿࠋ᱁Ꮚୖ࡟ᐃ⩏ࡉࢀࡓศᕸ㛵ᩘࡣ⇕㞧㡢ࢆᣢࡓࡎ㸪ࡲࡓὶయࢩ࣑࣮ࣗࣞࢩࣙࣥ࡜ྠᵝ࡟㡿ᇦศ๭࡟ࡼࡿ୪ิィ⟬ࡶᐜ࡛᫆࠶ࡿࠋࡋ࠿ࡋ㸪ࣈࣛࢯࣇ᪉⛬ᘧࡣᐇ✵㛫 3 ḟඖཬࡧ㏿ᗘ✵㛫 3 ḟඖࡢィ 6 ḟඖࢆᢅ࠺᪉⛬ᘧ࡛࠶ࡾ㸪ࢥࣥࣆ࣮ࣗࢱ࡛ゎࡃ࡟ࡣ⭾኱࡞ࣜࢯ࣮ࢫࢆᚲせ࡜ࡍࡿࠋࡇࡢࡓࡵ㸪ࡑࡢᡭἲࡢ㛤Ⓨࡣ࠶ࡲࡾ㐍ࢇ࡛࠸࡞࠸ࠋᐇ㝿㸪ࡇࡇᩘᖺࡢ HPC ࣉࣟࢪ࢙ࢡࢺ࡟ࡼࡿィ⟬ᶵ⎔ቃࡢ㣕㌍ⓗ࡞ྥୖ࡟ࡼࡗ࡚ᡭἲࡢ㛤Ⓨࡀ㐍ࡳ㸪ᐇ✵㛫 2 ḟඖཬࡧ㏿ᗘ✵㛫 3 ḟඖࡢ 5 ḟඖࢸࢫࢺࢩ࣑
࣮ࣗࣞࢩࣙࣥࡀࡼ࠺ࡸࡃᐇ⾜࡛ࡁࡿࡼ࠺࡟࡞ࡗࡓẁ㝵࡛࠶ࡿࠋᮏ◊✲ࡢ᭱⤊ⓗ࡞┠ⓗࡣ㸪ࣉࣛࢬ࣐ࢩ࣑࣮ࣗࣞࢩࣙࣥ࡜ࡋ࡚ࡣࠕḟࠎࠎࠖୡ௦ࡢᢏ⾡࡟࠶ࡓ
ࡿ 6 ḟඖ➨୍ཎ⌮ࣈࣛࢯࣇࢩ࣑࣮ࣗࣞࢩࣙࣥᡭἲࢆୡ⏺࡟ඛ㥑ࡅ࡚☜❧ࡋ㸪ࣉࣛࢬ࣐⛉Ꮫ࡟ᇶ࡙࠸ࡓᏱᐂኳẼࡢᐇ⌧࡟㈉⊩ࡍࡿࡇ࡜࡟࠶ࡿࠋࡑࡢࡓࡵ࡟ᮏࣉࣟࢪ࢙ࢡࢺ࡛ࡣ㸪5 ḟඖၥ㢟ࢆᢅ࠺ࣈࣛࢯࣇࢥ࣮ࢻࡢᛶ⬟ホ౯ࢆ㸪ᅜෆ࡛⌧Ꮡࡍࡿ㉸୪ิィ⟬ᶵࡢ୰࡛ࡶ≉࡟ࢥ࢔ᩘࡢከ࠸ᮾ
ி኱Ꮫ᝟ሗᇶ┙ࢭࣥࢱ࣮ࡢ FX10 ୖ࡛⾜࠺ࠋᮏ◊✲ࢢ࣮ࣝࣉ࡛ࡣࡇࢀࡲ࡛㸪ᮾி኱Ꮫ㸪ྡྂᒇ኱Ꮫ㸪஑ᕞ኱Ꮫ➼ࡢࢫ࣮ࣃ࣮ࢥࣥࣆ࣮ࣗࢱࢩࢫࢸ࣒࡟࠾࠸࡚඲ࣀ࣮ࢻィ ࢆ⾜ࡗ࡚ࡁ࡚࠾ࡾ㸪ࡲࡓ☢Ẽࣜࢥࢿࢡࢩࣙࣥ㸪ࢣࣝࣅࣥ㸫࣒࣊ࣝ࣍ࣝࢶ୙Ᏻᐃᛶ㸪ኴ㝧㢼㸫㠀☢໬ኳయ┦஫స⏝࡞࡝ࡢ≀⌮ィ⟬࡬ࡢ㐺⏝࡟㸪ୡ⏺࡛၏୍ᡂຌࡋ࡚࠸ࡿ[Umeda et al., 2010a,b; Umeda et al., 2011;Umeda, 2012; Umeda et al., 2012a]ࠋ2012 ᖺ 5 ᭶࡟ࡣ㸪ᮾி኱Ꮫ᝟ሗᇶ┙ࢭࣥࢱ࣮ࡢ኱つᶍ HPC
ࢳࣕࣞࣥࢪ࡟࠾࠸࡚ FX10 ࢆ⏝࠸ࡓィ ࢆ⾜ࡗࡓ㸦㟁☢ὶయࢥ࣮ࢻ࡟ࡼࡿᝨᫍ☢Ẽᅪࢩ࣑࣮ࣗࣞࢩࣙࣥᛶ⬟ᐃ㸪௦⾲㸸῝ἑ㸧ࠋ⾲ 1 ࡟ᐩኈ㏻♫〇ࢩࢫࢸ࣒ࢆ⏝࠸࡚ࡇࢀࡲ࡛࡟ィࢆ⾜ࡗ࡚ࡁࡓࣈࣛࢯࣇࢥ࣮ࢻࡢ༢ࢥ࢔ᛶ⬟ࢆ♧ࡍࠋᮏࣉࣟࢪ࢙ࢡࢺ࡟ඛ❧ࡕ㸪ࡇࢀࡽࡢ඲ࣀ࣮ࢻィ ࡟࠾࠸࡚᪤▱࡜࡞ࡗࡓၥ㢟Ⅼࢆ௨ୗ࡟♧ࡍࠋ
࣭FX1 ࡜ FX10 ࡛㸪ᐇຠຠ⋡࡟኱ࡁ࡞ᕪࡀ࠶ࡿࠋFX1 ࡢ 15%࡟ᑐࡋ࡚㸪FX10 ࡛ࡣ 11.6%ࡋ࠿ฟ
࡚࠸࡞࠸㸦5/22 ࡢィ㸧ࠋ࣭CX400 ୖ࡟࠾࠸࡚㸪ᐩኈ㏻♫〇ࢥࣥࣃ࢖ࣛ࡜࢖ࣥࢸࣝ♫〇ࢥࣥࣃ࢖࡛ࣛᛶ⬟ᕪࡀ࠶ࡿࠋ࢖ࣥ
ࢸࣝࢥࣥࣃ࢖ࣛࢆ⏝࠸ࡓሙྜࡢ 22.5%࡟ᑐࡋ࡚㸪ᐩኈ㏻ࢥࣥࣃ࢖ࣛࢆ⏝࠸ࡓሙྜ࡛ࡣ 14.4%
ࡋ࠿ฟ࡚࠸࡞࠸ࠋ
࣭FX10 ࡛ࡣ㸪1440 ࣀ࣮ࢻࡲ࡛ࡢィ࡟ࡣᡂຌࡋࡓࡀ㸪2400 ࣀ࣮ࢻ௨ୖࡣᐇ⾜࡛ࡁ࡞࠿ࡗࡓࠋ

ᮏࣉࣟࢪ࢙ࢡࢺ࡛ࡣ㸪௨ୖࡢၥ㢟Ⅼࢆඞ᭹ࡋ㸪4800 ࣀ࣮ࢻࡢィ ࢆ⾜࠺ࡇ࡜ࢆ┠ⓗ࡜ࡍࡿࠋ⾲㸯㸸ࣈࣛࢯࣇࢥ࣮ࢻࡢ༢ࢥ࢔ᛶ⬟ࠋ
᱁ᏊⅬᩘ 40x20x30x30x30/core ࡢᙅ࠸ࢫࢣ࣮ࣜࣥࢢ࡟ࡼࡿィ ⤖ᯝࠋ᫬㛫ࢫࢸࢵࣉᩘࡣ 5ࠋࢩࢫࢸ࣒⤒㐣᫬㛫ᐇຠᛶ⬟ᐇຠຠ⋡

FX10㸦ᮾ኱㸧5/22

227.7sec

1.71GFlops

11.6%

FX10㸦ᮾ኱㸧10/25㸦ᮏࣉࣟࢪ࢙ࢡࢺ㸧200.0sec

2.54GFlops

17.2%

FX1㸦ྡ኱㸧

290.0sec

1.50GFlops

15.0%

HX600㸦ྡ኱㸧

250.1sec

1.74GFlops

17.4%

RX200㸦஑኱㸧

167.1sec

2.60GFlops

22.2%

CX400㸦஑኱㸧ᐩኈ㏻ࢥࣥࣃ࢖ࣛ

153.2sec

3.11GFlops

14.4%

CX400㸦஑኱㸧࢖ࣥࢸࣝࢥࣥࣃ࢖ࣛ

97.7sec

4.87GFlops

22.5%

スーパーコンピューティングニュース- 48 -

Vol. 15, No. 1

2013㸰㸬ィ⟬ᡭἲࡢᴫせ ↓⾪✺ࣉࣛࢬ࣐ࡢ᣺ࡿ⯙࠸ࡣ㸪௨ୗࡢࣈࣛࢯࣇ㸦↓⾪✺࣎ࣝࢶ࣐ࣥ㸧᪉⛬ᘧ࡟ࡼࡗ࡚グ㏙ࡉࢀࡿࠋ

�f s � �f s q s � � � �f s
�v � �( E � v � B) � � 0     (1)�t�r ms�v
� ��� ��ࡇࡇ࡛ E 㸪 B 㸪 r ཬࡧ v ࡣࡑࢀࡒࢀ㟁ሙ㸪☢ሙ㸪఩⨨㸪㏿ᗘࢆ⾲ࡍࠋࡲࡓ㸪 f s (r , v , t ) ࡣ఩⨨㸫㏿ᗘ఩┦✵㛫࡟࠾ࡅࡿࣉࣛࢬ࣐⢏Ꮚࡢศᕸ㛵ᩘ࡛࠶ࡾ㸪 s ࡣ࢖࢜ࣥࡸ㟁Ꮚ࡞࡝✀㢮ࢆ♧ࡍࠋq s ࡜ ms ࡣࡑࢀࡒࢀ㟁Ⲵ࡜㉁㔞ࢆ⾲ࡍࠋ ࣉࣛࢬ࣐⢏Ꮚࡢศᕸ㛵ᩘࡣ㸪㟁☢ሙ࡟ࡼࡗ࡚ኚᙧࡍࡿࠋ㟁☢ሙࡢ᫬✵㛫Ⓨᒎࡣ௨ୗࡢ࣐ࢵࢡࢫ࢙࢘ࣝ᪉⛬ᘧ࡟ࡼࡗ࡚グ㏙ࡉࢀࡿࠋ��� 1 �E� � B � �0 J � 2
         (2.1)
� c �t�
�B�� E � �         (2.2)
�t� �        (2.3)
��E ��0�      (2.4)
�� B � 0�
ࡇࡇ࡛ J ࡣ㟁ὶᐦᗘ㸪 � ࡣ㟁Ⲵᐦᗘ㸪 � 0 ࡣ┿✵୰ࡢ㏱☢⋡㸪 � 0 ࡣ┿✵୰ࡢㄏ㟁⋡㸪 c ࡣග㏿ࢆ♧ࡍࠋࣈࣛࢯࣇ᪉⛬ᘧ(1)ࢆ㏿ᗘ✵㛫࡛✚ศࡍࡿ࡜㸪௨ୗࡢ㟁ⲴಖᏑ๎ࡀᚓࡽࢀࡿࠋ���
��� J � 0�t                  (3)�㟁☢ሙࢆኚ໬ࡉࡏࡿࡓࡵࡢ㟁ὶᐦᗘ J ࡣࣉࣛࢬ࣐ࡢ㐠ື࡟ࡼࡗ࡚⏕ࡌࡿࡀ㸪ࡇࢀࡣࣈࣛࢯࣇ᪉�
⛬ᘧ(1)ࡢ➨஧㡯࡟࠶ࡿᐇ✵㛫ࡢὶ᮰ v f s ࢆ㏿ᗘ✵㛫࡛✚ศࡍࡿࡇ࡜࡟ࡼࡗ࡚ồࡲࡾ㸪㟁ὶᐦᗘ�J ࡀ㟁ⲴಖᏑ๎(3)ࢆ‶㊊ࡍࡿ㝈ࡾ㸪࣏࢔ࢯࣥ᪉⛬ᘧ(2.3)ࡣ⮬ືⓗ࡟‶ࡓࡉࢀࡿࠋ௨ୖࡢ᪉⛬ᘧࡣ㸪ࣈࣛࢯࣇࢥ࣮ࢻ࡟࠾࠸࡚ゎ࠸࡚࠸ࡿࣉࣛࢬ࣐⢏Ꮚࡢ㐠ືㄽ᪉⛬ᘧ࡛࠶ࡾ㸪↓⾪✺ࣉࣛࢬ࣐ࡢ➨୍ཎ⌮࡜࿧ࡪࠋ
ࣈࣛࢯࣇ᪉⛬ᘧࡣ 4 ḟඖ௨ୖࡢࠕ㉸ከḟඖࠖࢆᢅ࠺᪉⛬ᘧ࡛࠶ࡾ㸪ࡑࡢࡲࡲࡢᙧ࡛ከḟඖᩘ್✚ศࢆ⾜࠺ࡢࡣ㠀ᖖ࡟ᅔ㞴࡛࠶ࡿࡓࡵ㸪₇⟬Ꮚศ㞳㸦operator splitting㸧ἲࡀྂࡃ࠿ࡽ⏝࠸ࡽࢀ࡚ࡁࡓ[e.g., Cheng & Knorr, 1976]ࠋ㐣ཤࡢ◊✲࡛ࡣ㸪ྛḟඖࡑࢀࡒࢀࢆ 1 ḟඖ⛣ὶ᪉⛬ᘧ࡟ศゎࡍࡿ᪉ἲࡀ᥇⏝ࡉࢀ࡚࠸ࡓࡀ㸪ᮏ◊✲࡛ࡣ㸪௨ୗࡢࡼ࠺࡟ᐇ✵㛫⛣ὶ㸪㏿ᗘ✵㛫⛣ὶ㸪㏿ᗘ✵㛫ᅇ㌿ࡢ 3 ࡘࡢ≀⌮ⓗ࡞₇⟬Ꮚ࡟ศ㞳ࡍࡿᡭἲࢆ㛤Ⓨࡋࡓ[Umeda et al., 2009]ࠋ

�f s � �f s
�v � � 0     �t�r
�f s qs � �f s
�E � �0 
�t ms �v
�f s qs � � �f s
�(v � B ) � � 0   
�t ms
�v

(4.1)
  (4.2)
 (4.3)ࡇࡢ₇⟬Ꮚศ㞳ࡣ㸪PIC ἲ࡟࠾࠸࡚ࢽ࣮ࣗࢺࣥ㸫࣮ࣟࣞࣥࢶᘧ㸦Ⲵ㟁⢏Ꮚࡢ㐠ື᪉⛬ᘧ㸧ࢆ᫬
㛫 2 ḟ⢭ᗘ࡛ゎࡃᡭἲ࡜ࡋ࡚ᗈࡃ⏝࠸ࡽࢀ࡚࠸ࡿ Boris ࢔ࣝࢦࣜࢬ࣒[Boris, 1970]࡟ᇶ࡙࠸࡚࠸
ࡿࠋ��

�

�

ᘧ(4.1)ཬࡧ(4.2)ࡣከḟඖࡢ⛣ὶ᪉⛬ᘧ࡛࠶ࡾ㸪 v ཬࡧ E ࡀࡑࢀࡒࢀ r ཬࡧ v ࡟౫Ꮡࡋ࡞࠸⥺ᙧ᪉⛬ᘧ࡛࠶ࡿࠋࣈࣛࢯࣇ᪉⛬ᘧࡢᩘ್✚ศ࡛ࡣ㸪ྂࡃ࠿ࡽࢭ࣑ࣛࢢࣛࣥࢪࣗἲࡀ᥇⏝ࡉࢀ࡚
スーパーコンピューティングニュース- 49 -

Vol. 15, No. 12013

ࡁࡓࠋ༶ࡕ㸪᱁ᏊⅬୖ࡟୚࠼ࡽࢀ࡚࠸ࡿศᕸ㛵ᩘࡢ್ࢆ᭦᪂ࡍࡿ㝿࡟㸪ḟࡢ᫬้࡟ࡑࡢ᱁ᏊⅬୖ࡟᮶ࡿ࡭ࡁศᕸ㛵ᩘࡢ್ࡢᗙᶆࢆࣂࢵࢡࢺ࣮ࣞࢫࡋ㸪㞄᥋ࡍࡿ᱁ᏊⅬࡼࡾ್ࢆෆᤄࡍࡿࠋ2000ᖺ௦๓༙ࡲ࡛ࡣ㸪ࢫࣉࣛ࢖ࣥ⿵㛫ࡸ CIP㸦Constraint Interpolation Profile㸧ἲ࡞࡝ࡢ㧗ḟ⿵㛫ἲࢆ⏝࠸ࡓࢭ࣑ࣛࢢࣛࣥࢪࣗἲࡀ⏝࠸ࡽࢀࡓ㸪ᐇ✵㛫 1 ḟඖ㸪㏿ᗘ✵㛫 1 ḟඖࡢ㸪2 ḟඖ఩┦✵㛫ࢩ࣑࣮ࣗࣞࢩࣙࣥࡀ୺ὶ࡛࠶ࡗࡓࠋᮏ◊✲ࢆࡣࡌࡵ࡜ࡍࡿ㏆ᖺࡢ◊✲࡛ࡣಖᏑᆺࢭ࣑ࣛࢢࣛࣥࢪࣗἲࢆ᥇⏝ࡋ㸪㉸ከḟඖࢩ࣑࣮ࣗࣞࢩࣙࣥ࡟ᡂຌࡋ࡚࠸ࡿࠋࡇࢀࡣ㸪ಖᏑ๎࡛࠶ࡿࣈࣛࢯࣇ᪉⛬ᘧ࡟࠾࠸࡚㸪ࣉࣛࢬ࣐ࡢᩘᐦᗘࡢಖᏑࢆᩘ್ⓗ࡟ಖドࡋ㸪ᩘ್ㄗᕪࢆ㍍ῶࡋࡼ࠺࡜ࡍࡿࡶ
ࡢ࡛࠶ࡿࠋPIC ࢩ࣑࣮ࣗࣞࢩ࡛ࣙࣥࡣ㸪㟁ⲴಖᏑ๎(3)ࢆཝᐦ࡟‶ࡓࡍࡇ࡜࡟ࡼࡗ࡚ᩘ್ㄗᕪࡢ
㍍ῶࢆ⾜ࡗ࡚࠾ࡾ㸪ࣈࣛࢯࣇࢩ࣑࣮ࣗࣞࢩ࡛ࣙࣥࡣಖᏑᆺゎἲࢆ⏝࠸ࡿࡇ࡜࡟ࡼࡾ㸪ᘧ(3)ࡀ⮬ືⓗ࡟‶ࡓࡉࢀࡿࡓࡵ㸪ಖᏑᆺゎἲࢆ᥇⏝ࡍࡿ࣓ࣜࢵࢺࡣ኱ࡁ࠸ࠋࡲࡓᘧ(4.1)ཬࡧ(4.2)ࡢከḟඖࡢ⥺ᙧ⛣ὶ᪉⛬ᘧࢆゎࡃୖ࡛㸪ィ⟬⢭ᗘࢆୖࡆࡿࡇ࡜ࡣ㠀ᖖ࡟㔜せ࡛࠶ࡿࡀ㸪2 ḟ⢭ᗘ௨ୖࡢᡭἲࢆ⏝࠸ࡓ㝿࡟⌧ࢀࡿᩘ್᣺ືࢆᢚไࡍࡿࡓࡵࡢ࣑ࣜࢵࢱࢆಖᏑᆺゎἲ࡟ᑟධࡍࡿࡢࡣẚ㍑ⓗᐜ࡛᫆࠶ࡿࠋᮏ◊✲࡛ࡣ㸪₇⟬Ꮚศ㞳࡟ࡼࡿᩘ್ᣑᩓࢆᢚไࡍࡿࡓࡵ࡟㸪ከḟඖࡢ⥺ᙧ⛣ὶ᪉⛬ᘧ࡟ᑐࡍࡿ₇⟬Ꮚ㠀ศ㞳㸦unspliting㸧ἲࢆ᪂ࡓ࡟㛤Ⓨࡋ࡚࠸ࡿ[Umeda et al., 2009]ࠋࡲࡓᮏ◊✲࡛ࡣ㸪↓᣺ືᛶཬࡧṇ್ᛶࢆಖドࡍࡿ࣑ࣜࢵࢱࢆ᪂ࡓ࡟㛤Ⓨࡋ㸪ᩘ್᣺ືࡢᢚไࢆ⾜ࡗ࡚࠸ࡿ[Umeda,2008; Umeda et al., 2012b]ࠋࡇࡇ࡛↓᣺ືࢫ࣮࣒࢟࡜ࡣ㸪࠶ࡿ༊㛫࡟࠾࠸࡚᪂ࡓ࡞ᴟ್㸦ᴟ኱㸪ᴟᑠ㸧ࢆ⏕ࡌࡎ㸪᪤࡟Ꮡᅾࡍࡿᴟᆅࢆ㸦࡛ࡁࡿࡔࡅ㸧ῶ⾶ࡉࡏ࡞࠸ࢫ࣮࣒࡛࢟࠶ࡿࠋᮏ◊✲࡛⾜ࡗࡓ࣋ࣥࢳ࣐࣮ࢡࢸࢫࢺ࡛ࡣ㸪㛤Ⓨ㏵ୖࡢ 5 ḟከ㡯ᘧࢆ⏝࠸ࡓ⿵㛫ࢆࢸࢫࢺⓗ࡟᥇⏝ࡋ࡚࠸ࡿࠋ
ᘧ(4.3)ࡣⲴ㟁⢏Ꮚࡢ㏿ᗘࡀ☢ຊ⥺࡟ࡼࡾ㐠ື࢚ࢿࣝࢠ࣮ࢆಖࡗࡓࡲࡲኚ໬ࡍࡿᅇ㌿᪉⛬ᘧࢆ⾲ࡍࠋ┤஺ᗙᶆ⣔࡟࠾ࡅࡿᅇ㌿᪉⛬ᘧࡣ๛యᅇ㌿ၥ㢟࡜➼౯࡛࠶ࡾ㸪⥺ᙧ⛣ὶၥ㢟࡜ྠᵝ࡟㸪ᩘ್ィ⟬࡟࠾࠸࡚᭱ࡶᇶᮏⓗ࡛࠶ࡿࡀ㸪ィ⟬⢭ᗘࡀ㔜せ࡜࡞ࡿၥ㢟࡛࠶ࡿࠋᮏ◊✲࡛᥇⏝ࡋ࡚࠸ࡿ back-substitution ἲ[Schmitz & Grauer, 2006]࡛ࡣ㸪Boris ࢔ࣝࢦࣜࢬ࣒[Boris, 1970]࡟ᇶ࡙࠸࡚㏿ᗘ✵㛫࡛ࡢ⢏Ꮚࡢ㌶㐨ࢆࣂࢵࢡࢺ࣮ࣞࢫࡋ㸪vx㸪 vy㸪 vz ᪉ྥࡑࢀࡒࢀࡢ₇⟬Ꮚࢆศ㞳ࡋ࡚ᅇ㌿㐠ືࢆゎ࠸࡚࠸ࡿࠋ๛యᅇ㌿ၥ㢟࡛ࡣ㸪⣔ࡢእഃ㸦㏿ᗘ✵㛫࡛ࡣ㏿ᗘࡀ㏿࠸ᡂศ㸧࡟࡞ࢀࡤ࡞ࡿ࡯࡝⛣ື㔞ࡀ኱ࡁࡃ࡞ࡾ㸪ࢡ࣮ࣛࣥ᮲௳ࡢᙳ㡪ࢆཷࡅࡸࡍࡃ࡞ࡿⅬ࡟ὀពࡀᚲせ࡛࠶ࡾ㸪௒ᚋ㸪㝜ゎἲࡸ₇⟬Ꮚ㠀ศ㞳ἲࡢ㛤Ⓨࡀᚲせ࡛࠶ࡿࠋ௨ୖࡢࡼ࠺࡟㸪ࣈࣛࢯࣇ᪉⛬ᘧࡢᩘ್ゎἲࡣᮍࡔⓎᒎ㏵ୖ࡛࠶ࡿࠋࡇࡢ኱ࡁ࡞ཎᅉࡣ㸪ࣈࣛࢯࣇࢥ࣮ࢻ࡛ᢅ࠺ḟඖࡀከ࠸ࡓࡵ࡛࠶ࡾ㸪㛤Ⓨࡸࢹࣂࢵࢢ࡟ࡓࡵ࡟኱ᐜ㔞ࡢඹ᭷࣓ࣔࣜ⎔ቃࡀᚲせ࡜࡞ࡿ࠿ࡽ࡛࠶ࡿࠋ୍᪉㸪࣐ࢵࢡࢫ࢙࢘ࣝ᪉⛬ᘧ(2.1)ཬࡧ(2.2)ࡣ㸪FDTD㸦Finite Difference Time Domain㸧ἲ࡜࿧ࡤࢀࡿ㟁☢ሙゎᯒἲࢆ⏝࠸࡚ゎࡃࠋFDTD ἲ࡛ࡣ㸪Yee ᱁Ꮚ[Yee, 1966]࡜࿧ࡤࢀࡿ staggered ᱁Ꮚࢆ⏝࠸࡚࠾ࡾ㸪ᘧ(2.4)ࡀ⮬ືⓗ࡟‶ࡓࡉࢀࡿࡼ࠺࡟≀⌮㔞ࡀ㓄⨨ࡉࢀ࡚࠸ࡿࠋࡲࡓ leap-frog࢔ࣝࢦࣜࢬ࣒࡟ᇶ࡙࠸࡚㟁ሙ࡜☢ሙࡢ᫬้ࡶ༙ࢱ࢖࣒ࢫࢸࢵࣉࡎࡽࡋ࡚࠾ࡾ㸪✵㛫⢭ᗘ࡜᫬㛫⢭ᗘࡣඹ࡟ 2 ḟ࡛࠶ࡿࠋࣈࣛࢯࣇࢩ࣑࣮ࣗࣞࢩ࡛ࣙࣥࡣ㠀ᖖ࡟ከࡃࡢ࣓ࣔࣜࢆᚲせ࡜ࡍࡿࡓࡵ㸪୪ิィ⟬ࡀᚲ㡲࡜࡞ࡿࠋࣈࣛࢯࣇࢥ࣮ࢻ࡛౑⏝ࡍࡿ≀⌮㔞ࡣ඲࡚᱁ᏊⅬୖ࡛୚࠼ࡽࢀ࡚࠾ࡾ㸪୪ิ໬࡟࠾࠸࡚ࡣ㡿ᇦศ๭ἲࡀ᭷ຠ࡛࠶ࡿࠋᅗ 2 ࡣᐇ✵㛫 2 ḟඖཬࡧ㏿ᗘ✵㛫 3 ḟඖࢆ౑⏝ࡍࡿࣈࣛࢯࣇࢥ࣮ࢻ࡟࠾ࡅࡿ୪ิ໬ࡢᴫᛕࢆ♧ࡍࠋᡃࠎࡢ┠ࡣ 4 ḟඖ௨ୖࡢ✵㛫ࢆㄆ㆑࡛ࡁ࡞࠸ࡀ㸪2 ḟඖᐇ✵㛫ࡢ

スーパーコンピューティングニュース- 50 -

Vol. 15, No. 12013

ྛ᱁Ꮚୖ࡟ 3 ḟඖ㏿ᗘ✵㛫㸦㏿ᗘศᕸ㛵ᩘ㸧ࡀᐃ⩏ࡉࢀ࡚࠸ࡿ࡜⪃࠼ࡿ࡜ศ࠿ࡾࡸࡍ࠸ࠋᮏ◊✲࡛ࡣᅗ 1 ࡢࡼ࠺࡟ᐇ✵㛫㸦x㸫y ᖹ㠃㸧࡟࠾࠸࡚ࡢࡳ㡿ᇦศ๭ࢆ⾜࠸㸪㏿ᗘ✵㛫ࡢ㡿ᇦศ๭ࡣ⾜ࢃ࡞࠸[Umeda et al., 2012a]ࠋࡇࢀࡣ㸪㟁Ⲵᐦᗘࡸ㟁ὶᐦᗘ࡞࡝ࡢ࣮࣓ࣔࣥࢺ㔞ࢆィ⟬ࡍࡿ㝿࡟ᚲせ࡞㏿ᗘ✵㛫ࡢ✚ศ࡟࠾࠸࡚㸪ྛᐇ✵㛫࡛ࡢࣜࢲࢡࢩࣙࣥฎ⌮ࢆ⾜ࢃ࡞࠸ࡼ࠺࡟ࡍࡿࡓࡵ࡛࠶ࡿࠋࡲࡓ࢜ࣉࢩࣙࣥ࡜ࡋ࡚㸪ࣉࣟࢭࢫ࢚࣓ࣞࣥࢺ㸦PE㸧ẖ࡟ OpenMP ࡟ࡼࡿࢫࣞࢵࢻ୪ิ໬ࡶチྍࡋ࡚࠸ࡿࠋ

ᅗ 1㸸5 ḟඖࣈࣛࢯࣇࢥ࣮ࢻ࡟࠾ࡅࡿ✵㛫㡿ᇦศ๭[Umeda et al., 2012a]ࠋ

㸱㸬ᛶ⬟ホ౯
3.1 ࢩࢫࢸ࣒㛫ཬࡧࢥࣥࣃ࢖ࣛ㛫ࡢᛶ⬟ᕪ࡟ࡘ࠸࡚ᮏࣉࣟࢪ࢙ࢡࢺ࡛⏝࠸ࡿࣈࣛࢯࣇࢥ࣮ࢻ࡛ࡣ㸪ᵝࠎ࡞⤌ࡳ㎸ࡳ㛵ᩘࢆከ⏝ࡋ࡚࠸ࡿࡓࡵ㸪ࡲࡎ㸪FX10, FX1 ཬࡧ CX400 ୖ࡟࠾࠸࡚㸪⾲ 2 ࡟♧ࡍ⤌ࡳ㎸ࡳ㛵ᩘࡢᛶ⬟ホ౯ࢆ⾜ࡗࡓࠋᅗ 2 ࡟♧ࡍࡼ࠺࡞⡆༢࡞ࣉࣟࢢ࣒ࣛࢆ౑⏝ࡋ㸪1 ⛊࠶ࡓࡾࡢ࿧ࡧฟࡋᅇᩘ㸦Mcall/sec㸧ࢆᐃࡋࡓ⤖ᯝࢆ⾲ 2 ࡟♧ࡍࠋࡇࡇ࡛㸪FX10, FX1 ཬࡧ CX400 ୖ࡛ࡢᐩኈ㏻♫〇ࢥࣥࣃ࢖࡛ࣛࡣfrt -Kfast
ࢆ㸪CX400 ୖ࡛ࡢ࢖ࣥࢸࣝ♫〇ࢥࣥࣃ࢖࡛ࣛࡣ
ifort –ip –ipo –O3 –xAVX
ࢆᣦᐃࡋࡓࠋ ࡲࡎ㸪FX1/FX10 ୖ࡛ࡣ CX400 ࡜ࡢࢡࣟࢵࢡ࿘Ἴᩘẚ(FX10㸸1.848GHZ㸪CX400㸸2.7GHz)
࡟ᑐࡋ࡚㸪ᩚᩘ࡟㛵㐃ࡋࡓ㛵ᩘࡢᛶ⬟ࡀᝏࡃ࡞ࡿࡇ࡜ࡀศ࠿ࡗࡓࠋ≉࡟㸪floor 㛵ᩘࡢᛶ⬟ࡀᴟ➃࡟ᝏࡃ㸪FX10 ࡟⮳ࡗ࡚ࡣ FX1 ࡼࡾࡶᛶ⬟ࡀᝏ࠸ࡇ࡜ࡀศ࠿ࡿࠋࡲࡓ CX400 ࡟࠾࠸࡚ࡶ㸪floor㛵ᩘࡢᛶ⬟࡟࠾࠸࡚ࠊᐩኈ㏻ࢥࣥࣃ࢖ࣛ࡜࢖ࣥࢸࣝࢥࣥࣃ࢖࡛ࣛᕪࡀฟࡓࠋࡇࢀ࡟㛵ࡋ࡚ࡣ㸪
floor 㛵ᩘᙉ໬࢜ࣉࢩࣙࣥ㸦-W0,-zmpa=Ailfunc-for-floor㸸㠀බᘧࠋ࡞࠾㸪ḟᅇࡢࢥࣥࣃ࢖ࣛࣂ࣮ࢪࣙࣥ࢔ࢵࣉ࡛㸪࢜ࣉࢩࣙࣥ࡞ࡋ࡛ࡶ floor 㛵ᩘࡀᙉ໬ࡉࢀࡿணᐃ㸧ࢆᣦᐃࡍࡿࡇ࡜࡟ࡼࡗ࡚࠶ࡿ⛬ᗘᨵၿࡍࡿࡇ࡜ࡀศ࠿ࡗࡓࡀ㸪FX10/FX1 ୖ࡛ࡣ㸪࢔ࢻࣞࢫィ⟬࡟ฟ᮶ࡿࡔࡅᩚᩘࡢ⤌ࡳ㎸ࡳ㛵ᩘࢆ⏝࠸࡞࠸ࡼ࠺࡟ࡋࡓ᪉ࡀ㸪ᛶ⬟ࡀྥୖࡍࡿࡇ࡜ࡀศ࠿ࡗࡓࠋᮏࣉࣟࢪ࢙ࢡࢺ࡛⏝࠸ࡿࣈࣛࢯࣇࢥ࣮ࢻ࡛ࡣ㸪ᩚᩘ࡛࢔ࢻࣞࢫィ⟬ࡋ࡚࠸ࡿ㒊ศࢆ඲࡚㸪ᆺኚ᥮ࢆྵࡴಸ⢭ᗘ࡛ィ⟬ࡍࡿࡼ࠺࡟ኚ᭦ࡋ㸪ᛶ⬟ࡀྥୖࡋࡓ㸦⾲ 1 ཧ↷㸧ࠋ

スーパーコンピューティングニュース- 51 -

Vol. 15, No. 1

2013real(kind=8) :: din(nx)integer(kind=4) :: nout(nx)real(kind=8) :: t1, t2
integer(kind=8) :: ii
call cpu_time(t1)do ii = 1,nx
nout(ii) = floor(din(ii))
end docall cpu_time(t2)print*,dble(nx/1000000)/(t2-t1),'M calls/sec'ᅗ 㸸㛵ᩘᛶ⬟ᐃࣉࣟࢢ࣒ࣛࡢ౛ࠋ

ࡇࡇ࡛ࡣ㸪IORRU 㛵ᩘ㸦ᘬᩘࡣಸ⢭ᗘᐇᩘ㸪ᡠࡾ್ࡣ༢⢭ᗘᩚᩘ㸧ࡢሙྜࢆ♧ࡋ࡚࠸ࡿࠋ⾲ 㸸ྛࢩࢫࢸ࣒࡟࠾ࡅࡿ㸪⤌ࡳ㎸ࡳ㛵ᩘ࿧ࡧฟࡋࡢ༢ࢥ࢔ᛶ⬟㸦༢఩ࡣ࣓࢞ᅇ⛊㸧ࠋᩚᩘ㛵ᩘࡣ༢⢭ᗘ㸪ᐇᩘ㛵ᩘࡣಸ⢭ᗘࢆ౑⏝ࡋࡓࠋ

FX10FX1

CX400㸦Fujitsu㸧 CX400㸦Intel㸧RAND

15.8M/s

19.6M/s

61.8M/s

97.7M/s

AINT

384.6M/s

286.8M/s

721.3M/s

724.7M/s

ANINT

714.3M/s

224.6M/s

721.3M/s

714.4M/s

ABS

833.3M/s

427.9M/s

728.3M/s

724.7M/s

SIGN

833.3M/s

435.6M/s

728.3M/s

724.8M/s

FLOOR

42.0M/s э 416.7M/s

92.9M/s

824.3M/s

1041.8M/s

NINT

416.7M/s

139.1M/s

903.8M/s

926.1M/s

INT

416.7M/s

318.9M/s

1000.2M/s

1020.6M/s

3.2 ๓ᅇ࡟඲ࣀ࣮ࢻィࢆኻᩋࡋࡓཎᅉ5 ᭶ 22 ᪥࡟ᐇ᪋ࡋࡓ኱つᶍ HPC ࢳࣕࣞࣥࢪ࡛ࡣ㸪1440 ࣀ࣮ࢻ㸦23040 ࢥ࢔㸧ࡲ࡛ࡢィ ࡟ࡣᡂຌࡋࡓࡀ㸪2400 ࣀ࣮ࢻ㸦38400 ࢥ࢔㸧௨ୖࡣᐇ⾜࡛ࡁࡎ㸪ࡲࡓཎᅉࡢ≉ᐃ࡟ࡣ⮳ࡽ࡞࠿ࡗࡓࠋ௒ᅇ㸦10 ᭶ 25 ᪥㸧ࡣ㸪2400 ࣀ࣮ࢻࢪࣙࣈࡢᐇ⾜ࢹ࣮ࢱࢆศᯒࡋ㸪ึᮇ≧ែ࡛᪤࡟್ࡀ NaN࡟࡞ࡗ࡚࠸ࡿ㒊ศࡀ࠶ࡿࡇ࡜ࢆ✺ࡁṆࡵࡓࠋᮏࣉࣟࢪ࢙ࢡࢺ࡛⏝࠸ࡿࣈࣛࢯࣇࢥ࣮ࢻ࡛ࡣ㸪ึᮇ≧ែ࡛ࣁࣜࢫ㟁ὶᒙ࡜࿧ࡤࢀࡿ☢໬ࣉࣛࢬ࣐ࡢᖹ⾮ゎ[Harris, 1962]ࢆタᐃࡋ࡚࠾ࡾ㸪ࡑࡢ୰࡛⏝࠸࡚࠸ࡿ 1/cosh2(x) ࡟࠾࠸࡚㸪NaN ࡀ᳨ฟࡉࢀ࡚࠸ࡿࡇ࡜ࢆ✺ࡁṆࡵࡓࠋ༶ࡕ㸪1.0E+154 ௨ୖࡢ್ࡢ 2 ஌࡛ Inf ࡀᡠࡾ㸪ࡑࡢ㏫ᩘ࡜
ࡋ࡚ NaN ࡀᡠࡿࠋࡇࢀࡣ FX10 ୖࡢᐩኈ㏻ࢥࣥࣃ࢖ࣛ≉᭷ࡢ⌧㇟ࡔ࡜ᛮࢃࢀࡿࡀ㸦CX400 ࡛ࡣ㉳ࡇࡽ࡞࠸ࡇ࡜ࡣ☜ㄆ῭㸧㸪ィ⟬㡿ᇦࡀ኱ࡁࡃ࡞ࡾᘬᩘࡀ኱ࡁࡃ࡞ࡗࡓࡇ࡜࡛⏕ࡌࡓ㸪኱つᶍィ⟬࡟ࡼࡗ࡚ึࡵ࡚ศ࠿ࡗࡓၥ㢟࡛࠶ࡿࠋ

スーパーコンピューティングニュース- 52 -

Vol. 15, No. 12013

3.3 ඲ࣀ࣮ࢻィ⤖ᯝၥ㢟Ⅼࡢಟṇࢆ⾜࠸㸪ࢥ࣮ࢻࡢᛶ⬟ ᐃࢆ⾜ࡗࡓࠋ ḟඖ௨ୖࡢ㉸ከḟඖၥ㢟ࡢࢩ࣑࣮ࣗࣞࢩࣙࣥ࡟ࡣ㠀ᖖ࡟ከࡃࡢ࣓ࣔࣜᐜ㔞ࢆᚲせ࡜ࡍࡿࡓࡵ㸪ࢥ࢔࠶ࡓࡾ࠶ࡿ࠸ࡣࣀ࣮ࢻ࠶ࡓࡾࡢ࣓ࣔࣜᐜ㔞ࢆ㐺ษ࡟タᐃࡍࡿᚲせࡀ࠶ࡿࠋ⌧Ꮡࡍࡿィ⟬ᶵୖ࡛ᐇ㝿࡟ࢩ࣑࣮ࣗࣞࢩࣙࣥࢆᐇ⾜ࡍࡿ㝿࡟ࡣ㸪303—603 ⛬ᗘࡢ㏿ᗘ✵㛫ࢆ౑⏝ࡋ࡚࠸ࡿࠋࡇࢀࡣ㸪ࢥ࢔࠶ࡓࡾࡢ࣓ࣔࣜᐜ㔞࡜ࡋ࡚
1GB—4GB ࢆ᝿ᐃࡋ࡚࠾ࡾ㸪㏿ᗘ✵㛫ࡢ᱁Ꮚᩘ࡜ᐇ✵㛫ࡢゎീᗘࢆᅛᐃࡋࡓࡲࡲ㸪ᐇ✵㛫ࡢ᱁Ꮚᩘ㸦ィ⟬㡿ᇦ㸧ࢆቑࡸࡍ㸪ᙅ࠸ࢫࢣ࣮ࣜࣥࢢ࡟ࡼࡾᐇ㝿ࡢィ⟬ࢆ⾜ࡗ࡚࠸ࡿࡇ࡜ࢆព࿡ࡍࡿࠋᮏࣉࣟࢪ࢙ࢡࢺ࡛ࡣ㸪ࢥ࢔࠶ࡓࡾࡢ᱁Ꮚᩘࢆᐇ✵㛫࡛ࡣ 40×20㸪㏿ᗘ✵㛫࡛ࡣ 30×30×30 ࡟タᐃࡍࡿࠋࡇࢀࡣ࣓ࣔࣜᐜ㔞࡛⣙ 1GB ࡟┦ᙜࡋ㸪ᐇ㝿ࡢィ⟬࡛౑⏝ࡋ࡚࠸ࡿ᱁Ꮚᩘ࡜࡯ࡰྠࡌ࡛࠶ࡿࠋ
ᮏࣉࣟࢪ࢙ࢡࢺ࡛⏝࠸ࡿࣈࣛࢯࣇࢥ࣮ࢻ࡛ࡣ MPI ࡜ OpenMP/⮬ື୪ิࡢࣁ࢖ࣈࣜࢵࢻ୪ิ
ࢆ᥇⏝ࡋ࡚࠸ࡿࡓࡵ㸪ࡲࡎ㸪ࣉࣟࢭࢫ୪ิ࡜ࢫࣞࢵࢻ୪ิࡢ࡝ࡢ⤌ࡳྜࢃࡏࡀ᭱㐺࠿ࢆ 1 ࣀ࣮ࢻ㸦16 ࢥ࢔㸧࡛ᐃࡋࡓࠋ⾲ 3 ࡟♧ࡍ࡜࠾ࡾ㸪2 ࢫࣞࢵࢻࡲࡓࡣ 4 ࢫࣞࢵࢻ࡛ᐇ⾜ࡋࡓሙྜࡀ

᭱ࡶ㏿ࡃ㸪ࡲࡓࣇࣛࢵࢺ MPI ࡢሙྜࡼࡾࡶ㧗㏿࡛࠶ࡿࡇ࡜ࡀศ࠿ࡿࠋࡉࡽ࡟⾲ 1 ࡜ࡢẚ㍑࡟ࡼ
ࡾ㸪ࣁ࢖ࣈࣜࢵࢻ୪ิࡢᐇຠᛶ⬟ࡀ༢ࢥ࢔ࡢᐇຠຠ⋡ࢆ㉸࠼࡚࠸ࡿࡇ࡜ࡀศ࠿ࡗࡓࠋ௨ୖࡢ⤖ᯝ࠾ࡼࡧࣉࣟࢭࢫᩘࡢ㍍ῶࡼࡾ㸪௨㝆ࡣ 4 ࢫࣞࢵࢻࡢࣁ࢖ࣈࣜࢵࢻ୪ิ࡟ࡼࡾᛶ⬟ィࢆ⾜ࡗࡓࠋ
1 ࢥ࢔࠶ࡓࡾ 1GB ࡟࣓ࣔࣜ౑⏝㔞ࢆᅛᐃࡋࡓ࡜ࡁࡢ㸪5 ḟඖࣈࣛࢯࣇࢥ࣮ࢻࡢᙅ࠸ࢫࢣ࣮ࣜ
ࣥࢢᛶ⬟ࢆᅗ 3 ཬࡧᅗ 4 ࡟♧ࡍࠋࡲࡓཧ⪃್࡜ࡋ࡚㸪HA8000 ཬࡧ FX1 ࡛ࡢィ⤖ᯝࡶ♧ࡍࠋ

ᅗ 3 ࡣࢥ࢔ᩘ࡟ᑐࡍࡿᐇຠ㏿ᗘࢆ⾲ࡋ㸪ᅗ 4 ࡣࢥ࢔ᩘ࡟ᑐࡍࡿᐇຠ୪ิ໬⋡㸦N ࢥ࢔ࡢᐇຠຠ⋡ࢆ 1 ࢥ࢔ࡢᐇຠຠ⋡࡛๭ࡗࡓ㔞㸧ࢆ⾲ࡍࠋࡲࡎ㸪࡝ࡢࢩࢫࢸ࣒࡟࠾࠸࡚ࡶ㸦ࡓࡔࡋ㸪ᮾி኱
Ꮫ HA8000 ࢆ㝖ࡃ㸧㸪ᛶ⬟ࡀ࡯ࡰࢫࢣ࣮ࣝࡋ࡚࠸ࡿࡇ࡜ࡀศ࠿ࡿࠋHA8000 ୖ࡛ࡣ㸪4096 ࢥ࢔࡜8192 ࢥ࢔࡛ᐇຠ୪ิ໬ຠ⋡ࡀ࡯ࡰྠࡌ࡛࠶ࡾ㸪ࡇࢀࡣࣀ࣮ࢻ㛫㏻ಙࡢࢺࣛࣇ࢕ࢵࢡࡀ㣬࿴ࡋ࡚࠸ࡿ࡜⪃࠼ࡽࢀࡿࠋFX1 ୖ࡛ࡣᐇຠ୪ิ໬⋡ࡀ 90%௨ୖ࡛࠶ࡗࡓࡀ㸪඲ࣀ࣮ࢻࢆ⏝࠸࡚ィࡋࡓሙྜ࡟ᴟ➃࡟ᛶ⬟ࡀຎ໬ࡍࡿ⌧㇟ࡀぢࡽࢀࡓ㸦ྡ኱ࡢ 3072 ࢥ࢔㸸ᅗ 3 ࡢ FX1 ࡢᐇ⥺ࢹ࣮ࢱ㸪
JAXA ࡢ 12032 ࢥ࢔㸸ᅗ 3 ࡢ FX1 ࡢⅬ⥺ࢹ࣮ࢱ㸧ࠋ୍᪉ FX10 ୖ࡛ࡣ㸪256 ࢥ࢔ࡲ࡛ࡣ 100%ࢆ㉸࠼࡚࠾ࡾ㸦๓㏙ࡢ㏻ࡾ㸪4 ࢫࣞࢵࢻ໬࡟ࡼࡾ༢ࢥ࢔ࡼࡾࡶᛶ⬟ྥୖࡀぢࡽࢀ࡚࠸ࡿࡓࡵ㸧
㸪1000ࢥ࢔௨ୖ࡛ࡶ 97%௨ୖࡢᐇຠ୪ิ໬⋡࡛ࢫࢣ࣮ࣝࡋ࡚࠾ࡾ㸪1PFlops ௨ୖ࡛ࡶ༑ศ࡟ࢫࢣ࣮ࣝࡍࡿࡇ࡜ࡀண᝿࡛ࡁࡿࠋࡲࡓᅗ 4 ࡟ぢࡽࢀࡿࡼ࠺࡟୪ิᛶ⬟࡟ከᑡࡢࣈࣞࡀ࠶ࡿࡶࡢࡢ㸪඲ࣀ࣮ࢻࢆ⏝࠸ࡓሙྜ࡛ࡶ FX1 ࡛ぢࡽࢀࡓࡼ࠺࡞ᛶ⬟ຎ໬ࡣ඲ࡃ㉳ࡇࡽ࡞࠿ࡗࡓࠋ᭱⤊ⓗ࡞⤖ᯝ࡜ࡋ࡚ࡣ㸪76,800 ࢥ࢔㸦4800 ࣀ࣮ࢻ㸧ࢆ౑⏝ࡋࡓ㝿ࡢᐇ⾜ᛶ⬟ཬࡧᐇຠຠ⋡ࡣࡑࢀࡒࢀ㸪192GFlops
࡜ 16.9%࡛࠶ࡾ㸪┠ᶆ࡜ࡋ࡚࠸ࡓຠ⋡ࢆୖᅇࡗࡓࠋ
⾲ 3㸸ࣈࣛࢯࣇࢥ࣮ࢻࡢ FX10 ༢ࣀ࣮ࢻᛶ⬟ࠋ
᱁ᏊⅬᩘ 80x40x30x30x30 ࡢィ ⤖ᯝࠋ᫬㛫ࢫࢸࢵࣉᩘࡣ 5ࠋࢫࣞࢵࢻᩘ124816

ࣉࣟࢭࢫᩘ

168421⤒㐣᫬㛫

199.8 sec

197.3 sec

197.3 sec

203.6 sec

205.9 sec

スーパーコンピューティングニュース- 53 -

Vol. 15, No. 1

2013㸲㸬ࡲ࡜ࡵ
H24 ᖺᗘ➨ 2 ᅇ኱つᶍ HPC ࢳࣕࣞࣥࢪ࡟࠾࠸࡚㸪Ᏹᐂࣉࣛࢬ࣐ࡢ➨୍ཎ⌮ࢩ࣑࣮ࣗࣞࢩࣙࣥ
ᡭἲ࡛࠶ࡿࣈࣛࢯࣇࢥ࣮ࢻࡢ 5 ḟඖࣂ࣮ࢪࣙࣥࡢᛶ⬟ホ౯ࢆ⾜ࡗࡓࠋ30,000 ௨ୖࡢ୪ิィ⟬ࢆ⾜࠼ࡿࢩࢫࢸ࣒ࡣ᪥ᮏ࡛ࡣᑡ࡞ࡃ㸪ࡲࡓ኱つᶍィ⟬࡟ࡼࡗ࡚ึࡵ࡚㟢ࢃ࡟࡞ࡿࢥ࣮ࢻࡢࣂࢢࡸ㸪ࢩࢫࢸ࣒㸦ࢥࣥࣃ࢖ࣛ㸧ࡢ㞀ᐖࡶ࠶ࡿࡓࡵ㸪ࡇࡢࡼ࠺࡞඲ࣀ࣮ࢻィࡣ HPC ศ㔝࡟࡜ࡗ࡚ᴟࡵ࡚㔜せ࡛࠶ࡿ࡜ゝ࠼ࡿࠋᮏࣉࣟࢪ࢙ࢡࢺࡢᡂᯝ࡜ࡋ࡚㸪ࡲࡎ㸪FX10 ࡣ㏻ಙࡀ㠀ᖖ࡟㏿ࡃ㸪ிࢥࣥࣆ࣮ࣗࢱࢡࣛࢫࡢ
୪ิᩘ࡛ࡶ༑ศ࡟ᛶ⬟ࡀࢫࢣ࣮ࣝࡍࡿྍ⬟ᛶࡀ࠶ࡿࡇ࡜ࡀศ࠿ࡗࡓࠋḟ࡟㸪FX10㸦SPARC64ࢳࢵࣉ㸧ࡢ≉ᚩ࡜ࡋ࡚㸪ᩚᩘ₇⟬ࡀ࢖ࣥࢸࣝ♫〇 CPU ࡟ᑐࡋ࡚㐜࠸ࡇ࡜ࡀᣲࡆࡽࢀ㸪ᩚᩘ₇⟬ࢆฟ᮶ࡿࡔࡅ⾜ࢃ࡞࠸ࡼ࠺࡟ࢳ࣮ࣗࢽࣥࢢࡍࡿࡇ࡜ࡀ㸪㧗࠸ᐇຠᛶ⬟ࢆᚓࡿࢥࢶࡢ 1 ࡘࡔ࡜ゝ࠼ࡿࠋࡲࡓ㸪୪ิᩘࢆୖࡆࡓ࡜ࡁ࡟ᛴ࡟ື࠿࡞ࡃ࡞ࡿሙྜ㸪㛵ᩘࡢᘬᩘࡀ኱ࡁࡃ࡞ࡗࡓࡇ࡜࡟
ࡼࡿ࣮࢜ࣂ࣮ࣇ࣮ࣟࡶཎᅉ࡜࡞ࡾᚓࡿࡢ࡛ࠊࢹࣂࢵࢢ᫬ࡢཧ⪃࡟ࡋ࡚㡬ࡅࢀࡤᖾ࠸࡛࠶ࡿࠋ

ᅗ 3㸸1GB/core ౑⏝᫬ࡢ 5 ḟඖࣈࣛࢯࣇࢥ࣮ࢻࡢ

ᅗ 4㸸1GB/core ౑⏝᫬ࡢ 5 ḟඖࣈࣛࢯࣇࢥ࣮ࢻࡢ

ᙅ࠸ࢫࢣ࣮ࣜࣥࢢᛶ⬟(1)ࠋ

ᙅ࠸ࢫࢣ࣮ࣜࣥࢢᛶ⬟(2)ࠋ

ࢥ࢔ᩘ࡟ᑐࡍࡿᐇຠ㏿ᗘࢆ♧ࡋ࡚࠸ࡿࠋ

ࢥ࢔ᩘ࡟ᑐࡍࡿᐇຠ୪ิ໬⋡ࢆ♧ࡋ࡚࠸ࡿࠋཧ ⪃ ᩥ ⊩Boris, J. P.: Relativistic plasma simulation-optimization of a hybrid code, Proc. Fourth Conf. Num. Sim.
Plasmas, ed. by J. P. Boris and R. A. Shanny, pp.3—67, Naval Research Laboratory, Washington D.
C. (Nov. 1970).
Cheng, C. Z., Knorr, G. : The integration of the Vlasov equation in configuration space, J. Comput. Phys.,Vol.22, No.3, 330—351 (1976).
Harris, E. G., On a plasma sheath separating regions of oppositely directed magnetic fields, Nuovo
Cimento, Vol.23, No.1, 115—121 (1962).
Schmitz, H., R. Grauer, R.: Comparison of time splitting and backsubstitution methods for integrating
Vlasov's equation with magnetic fields, Comput.Phys. Commun., Vol.175, No.2, 86—92 (2006).
Umeda, T.: A conservative and non-oscillatory scheme for Vlasov code simulations, Earth Planets Space,
Vol.60, No.7, 773—779 (2008).

スーパーコンピューティングニュース- 54 -Vol. 15, No. 1

2013

Umeda, T.: Effect of ion cyclotron motion on the structure of wakes: A Vlasov simulation, Earth Planets
Space, Vol.64, No.2, 231—236 (2012).
Umeda, T., Fukazawa, K., Nariyuki, Y., Ogino, T.: A scalable full electromagnetic Vlasov solver forcross-scale coupling in space plasma, IEEE Trans. Plasma Sci., Vol.40, No.5, 1421—1428 (2012a).
Umeda, T., Kimura, T., Togano, K., Fukazawa, K., Matsumoto, Y., Miyoshi, T., Terada, N., Nakamura, T.
K. M., Ogino, T.: Vlasov simulation of the interaction between the solar wind and a dielectric body,
Phys. Plasmas, Vol.18, No.1, 012908(7pp.) (2011).
Umeda, T., Miwa, J., Matsumoto, Y., Nakamura, T. K. M., Togano, K., Fukazawa, K., Shinohara, I.: Fullelectromagnetic Vlasov code simulation of the Kelvin-Helmholtz instability, Phys. Plasmas, Vol.17,No.5, 052311(10pp.) (2010a).
Umeda, T., Nariyuki, Y., Kariya, D.: A non-oscillatory and conservative semi-Lagrangian scheme with
fourth-degree polynomial interpolation for solving the Vlasov equation, Comput. Phys. Commun.,
Vol.183, No.5, 1094—1100 (2012b).
Umeda, T., Togano, K., Ogino, T.: Two-dimensional full-electromagnetic Vlasov code with conservative
scheme and its application to magnetic reconnection, Comput. Phys. Commun., Vol.180, No.3,365—374 (2009).
Umeda, T., Togano, K., Ogino, T.: Structures of diffusion regions in collisionless magnetic reconnection,
Phys. Plasmas, Vol.17, No.5, 052103(6pp.) (2010b).
Yee, K. S., Numerical solution of initial boundary value problems involving Maxwell's equations inisotropic media, IEEE Trans. Antenn. Propagat., Nol.AP-14, No.3, 302—307 (1966).

スーパーコンピューティングニュース- 55 -Vol. 15, No. 12013