C言語との邂逅：Fortranで「言語の壁」を壊すための第一歩
こんにちは。宇宙航空の現場で数値計算の泥沼と格闘し続けてきた身として、これからFortranに触れる皆さんに一つだけ伝えたいことがあります。
「Fortranは古い」なんて言葉は、ただの迷信です。今もスパコンの心臓部を動かし、物理シミュレーションの最前線で君臨しているのは、この言語の持つ「計算に対する純粋な効率」があるからです。
今回は、C言語やPythonを触ってきた皆さんが最初にぶつかる壁、「他言語とのデータ連携」についてお話しします。C言語の構造体やポインタを、Fortranの世界にどう持ち込むか。そのための鍵が `ISO_C_BINDING` モジュールです。
—
なぜ「型」を合わせる必要があるのか？
C言語の `int` と Fortranの `integer`。これらは、実はコンパイラやCPUの気まぐれでサイズが変わることがあります。例えば、Cの `int` が4バイトでも、Fortran側のデフォルト設定では8バイトになっていたり……。もしこの整合性が崩れたら、メモリ空間でデータが「ズレ」を起こし、計算結果はゴミになります。
これを防ぐための「共通言語」が `ISO_C_BINDING` です。これは、メモリ上のデータの「並び方」と「サイズ」を、C言語の規格と強制的に握り合わせるための魔法の道具箱だと思ってください。
—
実践：ISO_C_BINDINGで橋を架ける
まずは、C言語とFortranで共通して使える型定義の例を見てみましょう。
module my_data_types
use, intrinsic :: iso_c_binding
implicit none
! C言語の int に対応する型
integer(c_int), parameter :: my_int = 100
! C言語の double に対応する型
real(c_double), parameter :: my_val = 3.1415926535d0
end module my_data_types
ここがポイント！
- `use, intrinsic :: iso_c_binding`: これを宣言することで、C言語と互換性のある定数（`c_int`, `c_double` など）が使えるようになります。
- 型指定: `integer(c_int)` と書くことで、「この変数は環境によらず、C言語側の `int` と同じサイズだよ」とコンパイラに強く指示を出しているのです。
—
現場の知見：配列の「列優先順位」という罠
C言語経験者が最も苦しむのが「配列の並び順」です。
- C言語: 行優先（Row-major order）。`a[row][col]` のとき、メモリ上では `row` が固定で `col` が連続して並びます。
- Fortran: 列優先（Column-major order）。`a(row, col)` のとき、メモリ上では `col` が固定で `row` が連続して並びます。
もし、C言語からFortranへ配列のポインタを渡す際、この順序を考慮せずにループを回すとどうなるか。CPUのキャッシュミスが多発し、計算速度が10倍以上遅くなることも珍しくありません。
最適化の鉄則：
Fortranで配列を操作する際は、必ず「左側の添字」が内側のループに来るように回すこと。これが、スパコンで生き残るための最低限の作法です。
—
コンパイル時の注意：最適化フラグの魂
最後に、現場でよく使うビルド設定のヒントを置いておきます。皆さんが使っている `gfortran` や `ifort`（Intel Fortran）で、C言語との連携を考慮しつつ性能を引き出すフラグの例です。
GFortranでのビルド例
gfortran -O3 -march=native -ffree-form -o simulation_app main.f90
- `-O3`: 最適化の極み。これがないとFortranを使う意味が半減します。
- `-march=native`: あなたが今使っているPCのCPU命令セットをフル活用します。
- なぜこれが必要か: Fortranの配列操作は非常に最適化が効きやすいのですが、明示的にフラグを立てないと、コンパイラは安全側に倒して性能を抑えてしまうことがあるからです。
—
最後に：まずは「動く」ことから始めよう
最初は `ISO_C_BINDING` を使って、C言語で計算した値をFortranのサブルーチンに渡して出力する……そんな小さな一歩から始めてみてください。
Fortranは、あなたの意図をそのままCPUに伝えるための、極めて純粋で強力な道具です。言語の壁に恐れず、メモリの並びを意識し、計算の美しさを追求してください。
何かわからないことがあれば、またいつでも聞いてくださいね。皆さんの研究が、素晴らしい成果を生むことを応援しています！
コメント