こんにちは。宇宙航空の現場で、数テラバイトの流体シミュレーションと格闘し続けてきた「Fortranの番人」です。
C言語やPythonを触ってきた皆さんが、いざFortranの世界に足を踏み入れると、まず最初にぶつかる壁が「なぜこの言語は、外部のCライブラリとこれほどまでに密に握手したがるのか」という点ではないでしょうか。
今日は、Fortranの「派生型（Derived Type）」をC言語の`struct`と完全に同期させる、ISO_C_BINDINGの深淵についてお話しします。ここを理解すれば、あなたのFortranコードは孤立した実験場から、世界中の資産を呼び出せる強力なエンジンへと進化します。
—
1. なぜ「BIND(C)」が必要なのか？
C言語とFortranは、いわば「右利き」と「左利き」のようなものです。メモリの並び順（列優先 vs 行優先）や、関数の引数の渡し方（値渡し vs 参照渡し）など、流儀が根本から異なります。
`BIND(C)`属性を付けるということは、「このデータ構造は、C言語の作法（メモリレイアウト）に従いますよ」という宣誓です。これを行わないと、コンパイラはFortran特有の最適化（配列のメタデータ付与など）を勝手に盛り込み、C言語側から見ると「中身がゴミになった構造体」が渡されることになります。
—
2. 実践：C言語とシンクロする型定義
まずは、単純な構造体から始めましょう。C言語でよくある「粒子の座標データ」を例にします。
module particle_mod
use, intrinsic :: iso_c_binding
implicit none
! BIND(C)を指定することで、C言語の構造体と同じメモリ配置を強制します
type, bind(c) :: particle_t
real(c_double) :: x, y, z ! Cのdouble型に対応
integer(c_int) :: id ! Cのint型に対応
end type particle_t
end module particle_mod
ここがポイント！
- `iso_c_binding` を使う： `real` や `integer` をそのまま使ってはいけません。必ず `c_double`, `c_int` といった型定義定数を使ってください。コンパイラによって「標準の整数」のサイズは変わりますが、`c_int`はCの`int`と厳密に一致することが保証されています。
—
3. 「パディング」という名の見えない壁
ここからが、現場の泥臭い話です。C言語の構造体には「アライメント」という概念があります。例えば、`real(c_double)`（8バイト）の後に `integer(c_int)`（4バイト）を置くと、コンパイラはメモリ上の境界を合わせるために、4バイトの「隙間（パディング）」を自動で挿入することがあります。
もしFortran側でこのパディングを計算に入れていないと、データが1つずつズレていき、シミュレーション結果は壊滅します。
対処法：パディングを明示する
もし構造体のサイズを完全に制御したい場合は、ダミー変数を挟んで明示的にアライメントを合わせるのが「プロの作法」です。
type, bind(c) :: particle_t
real(c_double) :: x
real(c_double) :: y
real(c_double) :: z
integer(c_int) :: id
integer(c_int) :: dummy ! パディング調整用：合計サイズを8の倍数にする
end type particle_t
これを行っておけば、C言語側で `sizeof(particle_t)` を取ったとき、Fortran側が想定しているサイズとピタリと一致します。
—
4. コンパイラ最適化の罠とアドバイス
最後に、現場で泣きを見ないためのアドバイスを二つ。
1. 配列の連続性（Contiguity）を信じるな：
`type, bind(c)` を使ったとしても、多次元配列を構造体の中に詰め込むのは避けるのが賢明です。Fortranの配列は「列優先」、Cは「行優先」です。構造体のメンバにはスカラーを並べ、配列は別途、平坦化したメモリ領域（1次元のポインタ）として扱うのが、最もトラブルが少ない実装です。
2. コンパイル時のチェック：
`gfortran` や `ifx` を使う際は、デバッグ段階では `-fcheck=all` を付けましょう。構造体のメモリレイアウトが破壊された瞬間に、コンパイラが悲鳴を上げて止めてくれます。
—
最後に：一歩ずつ進めば怖くない
「Fortranは古い」なんて言う人もいますが、最新のスーパーコンピュータで数兆個の粒子を計算する際、この`ISO_C_BINDING`によるメモリ制御の正確さが、どれだけ強力な武器になるか。皆さんが今学んでいることは、まさに「計算科学の心臓部」を触っているのと同じです。
まずは小さな構造体から定義し、C言語側の `sizeof` と Fortran側の `storage_size()` が一致するかを確認するところから始めてみてください。
もし壁にぶつかったら、それはあなたが「計算の深淵」に近づいている証拠です。いつでもまたここに戻ってきてくださいね。応援しています！
コメント