たいていのMPIライブラリには、伝統的にmpi{cc,f90,cxx}といったコンパイラのラッパーが付いています。-Iとか-Lとか-lを適当に設定したうえで、コンパイラを呼び出してくれる便利なものですが、 これが最終的に呼び出すコンパイラを切り替える方法がライブラリによってバラバラなので、ざっくりとまとめておきます。 IntelMPIの場合 コマンド名で呼び出すコンパイラを切り替えています。 mpicc: gcc or cc mpicxx: g++ mpifc: gfortran mpigcc: gcc mpigxx: g++ mpif77: g77 mpif90: gfortran mpiicc: icc mpiicpc: icpc mpiifort: ifort gccのかわりにgcc44とかを呼びだければPATHの設定で先にくるようにしておけば良いはず*1 mpich2, mpich

CやC++で書かれたプログラムをMakeを使ってビルドする、というのはUnix/Linuxではよく行われていることだが、ちゃんとしたMakefileを書くのは意外と難しい。 例えば以下の3つのファイルからなるプログラムを考える。 foo.h: 関数fooの宣言がある。 foo.c: 関数fooの実装がある。 main.c: 関数fooを呼び出す。 /* foo.h */ void foo(int a); /* foo.c */ #include "foo.h" #include <stdio.h> void foo(int a){ printf("%d\n", a); } /* main.c */ #include "foo.h" int main(int argc, char **argv){ foo(10); return 0; } Makefileは例えば以下のように書ける。 PRO

==23525== Invalid free() / delete / delete[] / realloc() ==23525== at 0x1000089DF: free (in /usr/local/Cellar/valgrind/HEAD/lib/valgrind/vgpreload_memcheck-amd64-darwin.so) ==23525== by 0x100000F4F: main (sample.c:8) ==23525== Address 0x10081f510 is 0 bytes inside a block of size 4 free'd ==23525== at 0x1000089DF: free (in /usr/local/Cellar/valgrind/HEAD/lib/valgrind/vgpreload_memcheck-amd64-darwi

はじめに MPIを使って並列処理をする際に必要な準備やコマンド、ホストファイルの書き方を勉強したので、ここにざっとメモっときます。 以下の説明ではCentOS7で実行する場合を想定しています。 (追記: 2/22) ホストファイルの書き方の説明を一部修正しました。 (プロセス割り当ての設定の書き方が、MPICHとOpenMPIとで違っていました) MPIとは MPI (Message Passing Interface) は、並列コンピューティングを利用するための標準化された規格のことです。 並列処理には共有メモリ型と分散メモリ型の2通りの処理形態がありますが、MPIは分散メモリ型の並列処理を行います。 (共有メモリ型の並列処理を行うものにはOpenMPなどがあります) MPIの主な実装としてはMPICH (エムピッチ) とOpenMPI (オープンエムピーアイ) などがありますが、基本的

土曜日に後輩に頼まれてCのプログラムを見たところ、次の様なコードが。 main() { /* 処理1 */ printf("!"); /* 処理2 */ printf("!"); /* 処理3 */ }これで、セグメンテーションフォルトが起こるので、!の数を数えてみると、結果として1つだけ出力される。 そのため、処理2にエラーがあると思ってバグを探したけれど、バグが見つからなかったそうだ。 このように、printfなどで出力結果を見ながらデバッグを行う方法は、それなりに使える。しかし、そのためにはprintfの仕様を理解していなくてはならない。 printfの出力結果は、一般的にバッファへ保持される。 つまり、関数を呼べばすぐに結果がコンソールに出力されるわけではない。 バッファの内容をきちんとコンソールへ出力することを保証しなければならないのである。 コンソールへ出力するための方法として

CTF（セキュリティ関係のクイズ）でExploit（プログラムの脆弱性を突く）問題がさっぱり解けないので、練習した。 問題 こんな問題を自作した。 [kusano@madoka exploittest]$ ll 合計 16 -rwsrwxr-x 1 secret secret 5078 2011-10-05 23:50 crypt -rw-rw-r-- 1 secret secret 141 2011-10-05 23:49 crypt.cpp -rw------- 1 secret secret 88 2011-10-05 23:56 secret.txtcrypt.cpp #include <stdio.h> #include <string.h> int main() { char buf[256]; gets( buf ); memfrob( buf, strlen(buf) );

■ [PG] char 型を 30bit ビットシフトする C言語の話です。プログラミングのお話。 ここで言う char 型は 8bit の型定義としますが、これを 30 bit シフトするというのは問題ないのかどうか、という話。まあ仕事の最中にそんな記述を見て「これなんでコンパイルエラーにならんの？」となったのが発端なのですが。 C言語においてのビットシフトとは、まあ下みたいに書く訳です。 char x = 1; x = x << 2; char 型の変数 x （初期値1）を2bit分左にシフトした値にする、の意味です。結果的に x = 4 になるわけですが、まあその後の値はさておき。 ここでは、以下の式は許されるのかという話。 char x = 1; x = x << 30; ここで x はいくつになるかというと、 x = 0 です。 char 型は 8bit なので、当然ながらビット

■ [Game] ジルオール∞+(2) ウェイトレスさんイベントを終えて、でもそこまでとそこからがまだまだ長く、ようやく終盤かなあといったところ。 INT 上げまくったナッジが強くなりすぎて、デュアルスペルのホーリー×2で3000くらい出るので、ほぼ敵なし状態に。ただ AGI を上げてないのでターンが回ってくるのが遅いのはご愛敬。まあその間に全滅したりするんですけどね！ 取り敢えず主力は主人公とウェイトレスさんの二大高 AGI + 高 STR コンビ。ダブルでソニックアタック連発してるだけで中ボスクラスなら余裕で沈められる。それで沈められなくても先のホーリー×2が待ってるので、大抵は1ターンで戦闘終わり。まあアークナイトのソウルが強すぎですよっていう話ですかね。 VIT と MIN が低いから喰らうと涙目ではあるのですが。 一周目なので、ストーリーは後手後手の展開に。ゼネテスを助けられなか

■ [PG] volatile（C言語） 話題に困ってかつCDのレビューを書く気力もないときは適当にプログラミング話を綴っていけば良いんじゃないかと閃いた次第。まあ基本はやっぱりC言語になっちゃうんですけど。しかも初歩的な話だったりしますが。 そんなわけで、volatileです。 volatileって、あんま使わないですよね。たぶんVCとかで開発やってる限り、ほっとんどお世話にならないのではないかと。そもそもマルチスレッドプログラミングなんていきなりやらされることもなく、そんな理由からも新人研修とかでも飛ばされたりして。でも知らないと困るときが絶対的にあります。マルチスレッドプログラミングが必要なときですね。 volatileってのが何をするかっていうと、コンパイル時の最適化を抑止します。最適化ってのは、例えば変数aに0を代入して、次のステートメントで即座に変数aに100を代入するような処

プログラミングの入門書に書いてある通りに勉強していても、 自分のプログラムの結果がコンソールに出ているだけではなんだか味気なく感じてしまうことはありませんか？ ある程度プログラミングのやり方が分かってくると、画像を扱いたくなってくる人も多いでしょう。 それに、結果が画像の形で現れてくれると、プログラミングがなんだか一層楽しくなります。 でも、画像ファイルの扱い方は入門書にはほとんど書かれていませんね、 私も画像ファイルを扱いたいと思ったときどうすればいいのか調べるのに苦労しました。 ここでは、その画像ファイルのC言語での扱い方を解説したいと思います。 ここでやる内容は、画像ファイルのデータを書き出して、それをビューワなどで開いて確認するという形になります。 画像ファイルをC言語のプログラムで表示したりするというのとはちょっと違いますのでご注意を。 (こいつをするにはOSに依存したライブラリ

■ [PG] 可変引数を持つ関数マクロ（C言語） GNU Cの言語拡張とはいえ、そんなの出来たのか！ いや出来そうな気はしてたんだけど、じゃあそれってどう使うんだ？ って考えて挫折しました。まあ当然ながら、GCCのリファレンスに書いてあるわけで。（たまたま覗いてたら見つけた） http://www.sra.co.jp/wingnut/gcc/gcc-j.html#Macro%20Varargs GNU C では、関数と同じように、マクロも可変数の引数を取ることができる。その場合のマクロの定義構文は関数の場合に使われるものに良く似ている。次に例を示す。 #define eprintf(format, args...) \ fprintf (stderr, format , ## args) ##というプリプロセッサディレクティブがあるのは、マイナーながら知っていた。前後の文字列を合体させると