FortranとCの境界線：NULL終端文字列の「翻訳」を極める
こんにちは。宇宙航空の現場で何十年も数値計算コードの最適化を追いかけてきた身として、これからFortranの世界に足を踏み入れるあなたに、一つだけ大事なことをお伝えします。
Fortranは、いわば「生粋の数学者」です。対してC言語は「ハードウェアを直接操作する職人」です。この二人が会話をするとき、最もトラブルになりやすいのが「文字列の受け渡し」なんです。
C言語の文字列は`’\0’`（NULL文字）で終わるのがお約束ですが、Fortranの文字列は「何文字の箱か」という長さ（LEN）を厳密に管理しています。この文化の差を埋めるための作法を、現場の知見を交えて解説しますね。
—
1. なぜ「NULL終端」が必要なのか？
C言語の関数に文字列を渡すとき、相手は「どこで文字列が終わるか」をNULL文字を探して判断します。一方、Fortranは「メモリのこの番地から10バイト分が文字列だ」と決めてかかります。
ここで、FortranからCへ文字列を渡す際にNULL終端を忘れると、C側はメモリの果てまで文字を探し続け、運が悪ければセグメンテーション違反（Segmentation Fault）でプログラムが落ちます。この「翻訳ミス」を防ぐのが、`iso_c_binding` モジュールの役割です。
—
2. 実践：FortranからCへ文字列を渡す
まずは、C側が読める形に文字列を変換する定石を見てみましょう。`c_char` 型の配列を作り、最後に `c_null_char` を置くのが鉄則です。
subroutine send_to_c(fortran_str)
use iso_c_binding
implicit none
! 入力は通常のFortran文字列
character(len=), intent(in) :: fortran_str
! C言語用のバッファを用意
! +1はNULL終端文字の分です
character(kind=c_char), dimension(len(fortran_str) + 1) :: c_str
integer :: i
! 1文字ずつコピーして、最後にNULLを足す
do i = 1, len(fortran_str)
c_str(i) = fortran_str(i:i)
end do
c_str(len(fortran_str) + 1) = c_null_char
! ここでC言語の関数を呼び出す
! call c_print_function(c_str)
end subroutine
ここがポイント：
なぜ `len(fortran_str) + 1` としているのか。C言語は「メモリの連続性」を信じてNULLが出るまで読み進めます。Fortranの配列はメモリ上で連続していることが保証されているので、この書き方をすれば、C側は安心してそのアドレスをポインタとして扱えるんです。
—
3. パフォーマンスの泥臭い現実：スタックとヒープ
初心者がやりがちなミスは、関数のたびに巨大な文字列を毎回 `allocate` することです。数値計算のメインループでこれをやると、メモリ管理コストで計算速度がガタ落ちします。
- 静的領域を使い倒す: 文字列の長さが最大値として決まっているなら、`character(len=MAX_LEN, kind=c_char)` をあらかじめ定義しておくのが、メモリ断片化を防ぐ最強の最適化です。
- 配列の連続アクセス: Fortranは「列優先」です。文字列もメモリ上で連続していますが、もし複雑な構造体の中の文字列を扱う場合は、メモリアライメントを意識してください。不自然なパディングが入ると、SIMD命令（ベクトル演算）の効率が落ちることがあります。
—
4. 現場からのアドバイス：まずはこれから始めてみよう
最初は、「Fortranの文字列をC言語に渡す関数」を一つ書いてみることから始めてください。
! インターフェースブロックでC側の関数を宣言
interface
subroutine print_c_string(str) bind(c, name=”print_c_string”)
use iso_c_binding
character(kind=c_char), dimension() :: str
end subroutine
end interface
この `dimension()` という書き方は、「長さは教えないから、NULL文字を探して判断してね」というC言語へのメッセージです。
最後に
Fortranは古い言語だと思われがちですが、メモリレイアウトをこれほど直感的に制御できる言語は他にありません。文字列の受け渡し一つとっても、裏側で何が起きているか（ポインタの先頭アドレスを渡しているのか、配列全体をコピーしているのか）を意識できるようになると、あなたはもう中級者の入り口に立っています。
まずはコンパイルを通すこと。そして、`-g` オプションと `valgrind` などのツールを使い、メモリが正しく終端されているか確認する癖をつけてください。その泥臭いデバッグの積み重ねが、将来的にあなたを「計算科学の達人」へと押し上げてくれるはずです。
困ったことがあれば、またいつでも聞きに来てくださいね。一緒にいいコードを書きましょう！
コメント