とりあえず本家
http://www-unix.mcs.anl.gov/mpi/
http://www.h7.dion.ne.jp/~matsu/pc_cluster/mpi/intro/intro.html
http://p-grp.nucleng.kyoto-u.ac.jp/~uchida/mpi/
高校生でもこれほどの環境を手に入れられる時代になったのですね。
素晴らしい。
http://www.lam-mpi.org/
SMPP > PCクラスタ超入門
http://mikilab.doshisha.ac.jp/dia/smpp/
｢スカラー･チューニング講習会｣ ｢並列プログラミング(MPI)講習会｣で配布したテキスト
http://accc.riken.jp/HPC/training/text.html
OOMPI
http://www.osl.iu.edu/research/oompi/
とりあえず9までにいくつかある
MPI自体はほんとに難しいことはないから、リンク先でも読めばすぐにわかるよ。
mpiccを実行した後、mpirunすると
0 - MPI_SEND : Invalid rank 2
p0_23956: p4_error: : 8262
[0] Aborting program !
[0] Aborting program!
パイプ破壊
とエラーがでてしまいます。
何がいけないのでしょうか？；；
他のMPIプログラムはちゃんと動いたの？
もし動いたのなら、今回動かなかったプログラムのソースは載せられる？
> 0 - MPI_SEND : Invalid rank 2
Invalid rankなんだから、"2" に対応するpeがMPI_Sendで使用したコミュニケー
タに無いんじゃないの？
>>14
他のプログラム(πの計算）などはできました。
今回のは↓にＵＰしました。
http://www.geocities.jp/cirque55510/cellmpi.c
>>15
ちょっと
>>15
知識が足りなくて理解できない…
本でも買って出直してきます；；
for(source=1;source<mb;source++){
多分ここでしょ。でもmbをpに直しても間違いのような・・・
突っ込みどころたくさんあるんだけど、もうちょい勉強してみて。
これでも読んでください。
同志社大学の人が授業で使っている講義資料かな？
あそこは設備があるしね。。
http://mikilab.doshisha.ac.jp/dia/smpp/cluster2000/PDF/chapter02.pdf
ありがとうございます。
周りに聞ける人いなくて助かりましたー
勉強してがんばります！
「MPI並列プログラミング」（P.パチェコ著）も買ってみたんですが
結構難しくてちゃんと理解できるか不安に・・・。
日時：12月28日 13:30-
場所：同志社大学 京田辺校地 香知館 3F会議場
ttp://www.is.doshisha.ac.jp/SMPP/meetings/2004/041228.html
かなり行き詰まってきた；；
<NO ERROR MESSAGE> : Internal MPI error! Attribute in commun
icator is not a valid attribute
Special bit pattern 0 in attribute is incorrect.
[0] Aborting program !
[0] Aborting program!
原因がわかりません…ご教授のほどを
他のプログラムはプロセス数に限らず動いたんだけどなぁ
多分MPIの通信ルーチンの所だと思うけどね。
どういう状況でとまるのか書かないとわかりにくいッす。
愚考ながら、１のときと３のときでは止まる原因が違う、ということも。
ところで、MPICHのWIN版設定ってやってるページ少ないですね。
GUIのmpirunは接続異常は表示しないし…
MPICHって、GlobalArrays実装してますか？
英文必死こいて読んで見たけど、らしき記述が見当たらないのです。
・・・いや、職場で異様に自信満々で「共有メモリ使える」とってる香具師がいるもんで。
そもそも分散型プロセス通信の為の規格っつーか仕様ですよね？＞MPI
まさか「共有メモリ」なんて無いですよ。MPIって "Message Passing Interface" ですもん。
プロセッサ間のやりとりはすべてメッセージの送受信ですよ。
> ・・・いや、職場で異様に自信満々で「共有メモリ使える」とってる香具師がいるもんで。
ずいぶんと不幸な職場環境ですね……
・・・やっぱり？
おかしいとは思ったんですわ。
まあ、もうちょっと読み進めては見ますがせっかくだし。
それはそれとして、RMAのある実装って何かありますか？
MPICHにはないようで。
うまく使えばかなりおいしそうなんで、いじってみたい・・・
ハード的にはネットワークの使えるパソコン数台と、ルーター
ソフト系は上のほうの紹介サイトぐるっと回ってライブラリをとってくる。
あとはメインのプログラムの開発環境、これぐらい?
何がいけないのでしょうか？解決方法を教えてください。
/tmp/ccGU5ztX.o(.eh_frame+0x11): undefined reference to `__gxx_personality_v0'
collect2: ld returned 1 exit status
g++でコンパイルやリンクしてる？
このクラスのデストラクタがランク0以外働かないのですが、何故なのでしょうか？
mpiCCでいける
まずはMPIの仕様書を読もうよ
コンパイルする詳しい方法分かりますか？
include 'mpif.h'で定義ファイルを読み込むんですよね？
Unixは使えないのですか？
使えないです。unixでないとできませんか？
threadを使ったプログラムは経験があるのですが、MPIでの作成方法がよくわかりません。
なんでwhile(!done){}あたりのギミックが必要なの?
( ∀`) 彡 >>56 教えてやってもいいが
((≡O ⊃
| | おまい！左腕を上下 右腕を前後
し⌒J 同時にやってみろ！話はそれからだ！
他のｎ個のCPUが処理が完全同期で終わる保証はどこにもないだろ。
だったら全部の処理が終わるまで待つのは当たり前だろ。
それともそんなことすら解らないのか？プログラム組むの辞めろ。お前は
もう、日本から出て行け
> なんでwhile(!done){}あたりのギミックが必要なの?
単にコードの書き方がまずいだけと思われ。
この内容ならもっと見やすくシンプルにできるよね。
なんかわかったような気がする。
コメントアウトの部分のソースを無視していました。
頭悪くてすみませんでした。
>>58さん
MPI_Reduce()で全プロセスが終わるの待ってるのかと思ってたんですけど
違うんですか？
p5_30033: p4_error: interrupt SIGSEGV: 11
Broken pipe
Broken pipe
というエラーが出ます。これは何がいけないのでしょうか？
個人的には、ファイルアクセスで問題が起きていると思っているのですが、エラーの原因が未だつかめません。
>>58 そいつ、単なるバカだから相手にすんな
配列の範囲外にアクセスしているとか
MPI以外のエラーを疑え
問題ないと思っていたのですが、gcc4にアップデートして以来、原因不明で0番プロセス以外が実行中に落ちます。
gcc3 を使ってた方が安全だよ。
mpd.confにsecretwordを書いた後どうすればいいかちょっと分かりません
教えてもらえないでしょうか？
○ ちょっと教えてもらえないでしょうか？
C++でサーバアプリで計算しながら、クライアントアプリにソケット通信するようなシステムに
適用できますか？
サーバー側は１０００件/分ぐらいのトランザクションで、クライアント側はサマリー情報を
半リアルで表示させるだけといったかんじです。
http://www-lab.imr.tohoku.ac.jp/~t-nissie/computer/wolfchen/doc/HomeWolfchenProject.html
mpirun -c 16 myprog.exe
で実行すると、16個のプロセスが動き始めるのですが、問題があります。
それは、プログラムを実行すると、データファイルを読み込んで、処理が始まるのですが、データファイルはいくつかの大きさのものを用意しているのですが、
あるデータファイルでは、正しく計算が完了するのですが、あるデータファイルでは
計算開始直後にエラーで終了してしまいます。調べたところ、MPI_Bcast()で
エラーしていました。
また、正しく計算できるデータファイルを使って、
mpirun -c 8 myprog.exe
のように変えるとエラーする場合もあります。
そこで教えていただきたいのですが、そもそも1CPUマシンで、このような並列計算を
実行することは、何か問題がありますでしょうか？
よろしく。
メモリ不足とかじゃなきゃ、プログラムのバグだろうね
>基本的に動かないとまずい。
>メモリ不足とかじゃなきゃ、プログラムのバグだろうね
やはりそうですか。
実際、プログラムを単純化して、単に１６プロセスでMPI_Bcast()するだけというプログラムでは、問題なく動きます。
１６以外のどんな数字でも問題ありません。なので、プログラムのバグの可能性が高いのですが、一つ気になる点は、
エラーするケースでも、エラーしないケースでも、いずれの場合でもプロセスをps -ef などでモニタしていると、
実行中に１６プロセスそれぞれの親プロセスIDがころころ変わっているのが確認できるのですが、何かそのあたり
に関係しているのではと思うのですが、どなたか詳しい方、ヒントなどお願いします。
プログラムのバグ。
MPI_Barrier()で、きちんと同期とって
MPI_Bcast()で起こってるなら
MPI実装を疑ってもいい。
まあ、なんだ。
漏れには「FedoraのMPIでエラー終了」じゃ、
どのMPI環境なのかも
どんな異常終了かもサパーリですよ。
(1) 単純化したプログラムと、エラーが起こる実プログラムで、
送信データの大きさは同じですか？
(2)どういうエラーメッセージが出ますか?
(3)mpi_bcastの前後にprintfでランクの値を書き出した場合、
どうなりますか？
（全プロセスがmpi_bcast前しか表示しないのか、
それともプロセスによってはmpi_bcast後も表示するのか）
(4)エラーが起こるデータファイルで、プロセス数が16より少ないときは
どうなりますか？
>(1) 単純化したプログラムと、エラーが起こる実プログラムで、
>送信データの大きさは同じですか？
はい。MPI_INTを一個です。
>(2)どういうエラーメッセージが出ますか?
今自宅なので、正確には覚えていませんが、英語のメッセージで、
『どれかのプロセスでエラーしました』とか言う感じだったと思います。
>(3)mpi_bcastの前後にprintfでランクの値を書き出した場合、
>どうなりますか？
正しく、自分のランク値が表示されます。
>（全プロセスがmpi_bcast前しか表示しないのか、
>それともプロセスによってはmpi_bcast後も表示するのか）
全プロセスで、MPI_Bcast前に正しいランク値が表示されますが、
MPI_Bcast後には、どのプロセスも何も表示されません。
>どうなりますか？
ええっと、その後、いろいろ試してみたら以下のようになりました。
分割数 大きいデータ 小さいデータ
2 OK エラー
4 OK OK
5 エラー エラー
6 エラー エラー
7 エラー エラー
8 エラー OK
9 エラー OK
10 エラー OK
11 エラー OK
12 OK OK
13 エラー OK
14 OK OK
15 エラー OK
16 エラー OK
20 エラー OK
というわけで、よく分かりません。
ちなみに、MPI_Barrier同期は、まったく問題なく動きます。
今のところ、MPI_Bcastを使うのをやめて、MPI_Send/MPI_Recvで同等の処理を
することで、エラーを回避できたようなのですが、MPI_Bcastがエラーする
原因が分からないので、ちょっと気がかりです。
>はい。MPI_INTを一個です。
確認ですが、一個というのは、下記の2つ目の引数が1という意味ですか？
MPI_Bcast(a,1,MPI_INT,〜);
(6) 上記MPI_BCASTの4つ目の引数(送信元プロセスのランク値)は、
全プロセスで同じ値になっていますか？
また、2,3つ目の引数(送受信バッファーの個数と型)も
全プロセスで同じ値になっていますか？
(7)下記の表の意味ですが、1CPUで2プロセス(mpirun -np 2)で実行し、
小さいデータを入力させた場合でもエラーになるという意味ですか？
またこれは何度やっても同じ結果ですか？
（つまり同じケースでもOKになったりエラーになったりすることはないか）
分割数 大きいデータ 小さいデータ
2 OK エラー
(9) 実プログラムが下記のような構成になっているとして
(a)の部分をすべてコメントにすれば、
（単純化したプログラムではうまくいくとのことなので）
うまく動くと思いますが、
その後で、(a)の部分を前半から少しずつ
コメントをはずしてはテストしを繰り返せば、
どこかで再びエラーとなるはずです。
それによって、エラーの原因となっている個所が特定できるかもしれません。
あるいは逆に、(b)と同等の通信を(送信バッファーの中身は空でも可)
まず(a)の直前で行い、うまくいったら、
(a)内を少しずつ後ろへずらしていっても同じテストができます。
MPI_Init();
(a) 入力ファイルの読み込み等
(b) MPI_Bcast(a,1,MPI_INT,〜);
(c) 計算
MPI_Finalize();
>はい。MPI_INTを一個です。
確認ですが、一個というのは、下記の2つ目の引数が1という意味ですか？
MPI_Bcast(a,1,MPI_INT,〜);
また、下記の大きいデータの場合も小さいデータの場合も1ですか？
分割数 大きいデータ 小さいデータ
2 OK エラー
>(5)の訂正
>確認ですが、一個というのは、下記の2つ目の引数が1という意味ですか？
>MPI_Bcast(a,1,MPI_INT,〜);
はい、そうです。
>また、下記の大きいデータの場合も小さいデータの場合も1ですか？
>分割数 大きいデータ 小さいデータ
>2 OK エラー
はい。
>(6) 上記MPI_BCASTの4つ目の引数(送信元プロセスのランク値)は、
>全プロセスで同じ値になっていますか？
はい、そうなっています。printfで確認しました。
>また、2,3つ目の引数(送受信バッファーの個数と型)も
>全プロセスで同じ値になっていますか？
はい、同じです。
>(7)下記の表の意味ですが、1CPUで2プロセス(mpirun -np 2)で実行し、
>小さいデータを入力させた場合でもエラーになるという意味ですか？
はいそうです。
>またこれは何度やっても同じ結果ですか？
>（つまり同じケースでもOKになったりエラーになったりすることはないか）
再現されます。同じデータ、同じ分割数で、数回実験しましたが、その時は再現しました。
でも、何十回も同じ条件でテストしたわけではありません。
>(8) 問題のMPI_Bcastに到達する前に、MPIの他の通信は行っていますか？
他に行っているのは、MPI_Barrierのみです。同期をたくさん入れれば問題が解決するんじゃあ
ないかと思って、たくさん入れました。
>(9) 実プログラムが下記のような構成になっているとして
>(a)の部分をすべてコメントにすれば、
これらのテスト結果の報告は、後日させていただきます。
オプションをちょっと教えていただけないでしょうか？
osはfedora core4 コンパイラーにintel fortran compilar 9.0
mpich-1.2.7p1 を使用しています。
よろしくおねがいします。
if ( myid .eq. 0 ) then
call MPI_IRecv(retunum,90,MPI_REAL8,1,itag(myid+2),MPI_COMM_WORLD,ireq2,ierr)
elseif ( myid .eq. 1 ) then
call MPI_ISend(retunum,60,MPI_REAL8,0,itag(myid+1),MPI_COMM_WORLD,ireq2,ierr)
end if
if ( myid .eq. 0 ) then
call MPI_IRecv(iretunum,90,MPI_INTEGER,1,itag(myid+2),MPI_COMM_WORLD,istaus,ierr)
elseif ( myid .eq. 1 ) then
call MPI_ISend(iretunum,60,MPI_INTEGER,0,itag(myid+1),MPI_COMM_WORLD,ierr)
end if
以下のように、何故ですか?わかる人教えて下さい。
test-calcul2.f:99: warning:
call MPI_IRecv(retunum,90,MPI_REAL8,
1
test-calcul2.f:106: (continued):
call MPI_IRecv(iretunum,90,MPI_INTEGER,
2
Argument #1 of `mpi_irecv' is one type at (2) but is some other type at (1) [info -f g77 M GLOBALS]
> if ( myid .eq. 0 ) then
> call MPI_IRecv(iretunum,90,MPI_INTEGER,1,itag(myid+2),MPI_COMM_WORLD,istaus,ierr)
> elseif ( myid .eq. 1 ) then
> call MPI_ISend(iretunum,60,MPI_INTEGER,0,itag(myid+1),MPI_COMM_WORLD,ierr)
> end if
if ( myid .eq. 0 ) then
call MPI_IRecv(iretunum,90,MPI_INTEGER,1,itag(myid+2),MPI_COMM_WORLD,ireq2,ierr)
elseif ( myid .eq. 1 ) then
call MPI_ISend(iretunum,60,MPI_INTEGER,0,itag(myid+1),MPI_COMM_WORLD,ireq2,ierr)
end if
です。
もし同じエラーが出たとしたら、
同じサブルーチンの同じ引数に異なる型を指定しているという
警告で、MPIとは関係ないと思います。
program test
call sub(a)
call sub(i)
end
subroutine sub(a)
a=1
end
-2つのif文を１つにまとめた方がすっきりすると思います。
-例えばランク0プロセスが同じireq2を(mpi_waitする前に)2回使うのはまずいです。
片方を例えばireq1にしてください。
また2回のcall mpi_xxxの後でmpi_waitが2回必要となります。
あるいはmpi_sendrecvを使う手もあります。
まさにそれでした。
>>94
なるほど、了解です。
ところで、送信バッファーに３次元の配列を持たせたいのですが
送受信がうまくいきません。何故かわかる方いますか？
if ( myid .eq. 0 ) then
call MPI_Recv(test(i,j,1),3,MPI_REAL8,1,itag,
& MPI_COMM_WORLD,istaus,ierr)
else
call MPI_Send(test(i,j,1),3,MPI_REAL8,0,itag,
& MPI_COMM_WORLD,ierr)
endif
test(1〜50,1〜100,1〜3)です。
(2) 送る要素数は3個ですか？
(3) istatusは大きさmpi_status_sizeで配列宣言してますか？
(4) itagには何が指定されてますか？
(5) 実行するとどういうエラーになりますか？
(1) test(i,j,1)のiとjには何が指定されてますか？
i=1,100とj=1,300です。DOループで指定されます。
(2) 送る要素数は3個ですか？
はい、南北方向・東西方向・上下方向の意味なので、３個です。
(3) istatusは大きさmpi_status_sizeで配列宣言してますか？
integer istatus(MPI_STATUS_SIZE)で宣言しています。
(4) itagには何が指定されてますか？
itag=1です。
(5) 実行するとどういうエラーになりますか？
テストプログラムで行った結果、エラー文は出ませんでしたが
入っているはずの値が入っていませんでした。
たとえば、10.5が入るはずなのに、0.0のままであったりします。
（Fortranの場合は(A)が送られます。）
(A)
test(1,1,1)
test(2,1,1)
test(3,1,1)
(B)
test(1,1,1)
test(1,1,2)
test(1,1,3)
余談ですが、通信回数が多いと通信時間がかかるので、
普通は送りたいデータをまとめて1回で送るようにします。
従ってループ内で通信を行うということはあまりしません。
TOP カテ一覧 スレ一覧 削除依頼 ▲
・ 次のスレ
【会津】パソコン甲子園2004【若松】
0から始める2chブラウザfor超漢字 "2ch de BTRON"
Ｃ言語は関数ができなくても、理解可能か？
Cで書くかアセンブラで書くか・・・