Fortranは「隠す」ことで速くなる：MODULEにおけるカプセル化の真実
数値計算の世界へようこそ。CやPythonを使いこなしてきた皆さんなら、「なぜ今さらFortranなのか」と疑問に思うかもしれません。答えはシンプル。「計算機が持つポテンシャルを、もっとも素直に引き出せる言語だから」です。
しかし、Fortranを書き始めると最初にぶつかる壁が「名前空間の汚染」と「ビルド時間の増大」です。今日は、モダンFortranの心臓部である`MODULE`を使い、`PRIVATE/PUBLIC`属性で堅牢かつ爆速なコードベースを作るための「現場の作法」をお話ししましょう。
—
なぜ「隠す」必要があるのか？
Pythonなら`from module import `で全てを取り込んでも、せいぜいメモリを食う程度かもしれません。しかし、Fortranの世界では事情が違います。
コンパイラ（Intel FortranやGNU Fortranなど）は、`PUBLIC`と宣言された変数や手続きに対し、「どこからでもアクセスされる可能性がある」という前提で最適化を行います。対して`PRIVATE`に設定されたものは、「このモジュールの外からは一切触らせない」という契約を結ぶことになります。
この契約により、コンパイラは「この変数はモジュール内でしか使われていないから、レジスタに保持し続けよう」「この関数はここだけでしか呼ばれないから、インライン展開してしまおう」といった大胆な最適化を安心して行えるようになります。
ステップ1：まずは「全部隠す」のがプロの流儀
Fortranのモジュールは、デフォルトでは全てが`PUBLIC`です。これは大きな建物で、全ての部屋の鍵が開いているようなもの。まずはこれを「デフォルトで閉鎖」する書き方から始めましょう。
MODULE physical_constants
! 最初に「デフォルトで全て隠す」と宣言する
IMPLICIT NONE
PRIVATE
! 外から使ってほしい定数だけをPUBLICにする
PUBLIC :: PI, GRAVITY
REAL(8), PARAMETER :: PI = 3.14159265358979323846d0
REAL(8), PARAMETER :: GRAVITY = 9.80665d0
! 外には見せたくない、計算途中の「中間定数」
! これをPRIVATEにすることで、名前の衝突を完全に防ぐ
REAL(8), PARAMETER :: INTERNAL_TOLERANCE = 1.0d-12
END MODULE physical_constants
この`PRIVATE`という記述が一本あるだけで、あなたのコードは格段に安全になります。特に大規模な数値解析プログラムでは、何千行にも及ぶコードの中で「うっかり同じ変数名を使ってしまう」という悲劇を未然に防ぐ生命線になります。
—
ステップ2：コンパイルの「待ち時間」を短縮する
コンパイル時間が長いと、デバッグのモチベーションは下がりますよね。`MODULE`を細かく分けて`PRIVATE`を適切に設定すると、コンパイラが「このモジュールに依存しているのはここだけだな」と即座に判断できるようになります。
特に、`USE`文を使う際は以下のように絞り込むのが鉄則です。
! 悪い例：全部取り込んでしまい、依存関係が複雑化する
USE physical_constants
! 良い例：必要なものだけを明示する（可読性も向上します）
USE physical_constants, ONLY: PI, GRAVITY
こうすることで、コンパイラは無駄なシンボルテーブルの検索を行う必要がなくなり、ビルド時間が劇的に短縮されます。大規模な解析コードであればあるほど、この塵も積もれば山となる最適化が効いてきます。
—
ステップ3：現場のデバッグで役立つ小話
最後に、現場で血を流した者からのアドバイスです。
C言語の`static`と似ていますが、Fortranの`PRIVATE`は「モジュールという単位」に対して非常に強力に働きます。特に数値計算でよく使う「作業用配列（ワークエリア）」を`PRIVATE`で隠しておくと、並列計算（OpenMPなど）を導入する際に、「どの変数が共有され、どの変数がスレッドごとに独立しているべきか」という設計ミスを減らすことができます。
今日から皆さんのコードに加えるべきルール：
1. すべての`MODULE`の冒頭には`IMPLICIT NONE`と`PRIVATE`を書く。
2. 外に見せたいものだけを`PUBLIC`にする。
3. `USE`するときは`ONLY`で必要なものだけを呼び出す。
これだけで、あなたの書くFortranコードは「モダンで堅牢なプロの設計図」に一歩近づきます。最初は面倒に感じるかもしれませんが、計算結果が怪しい時に「公開範囲」を制限しておくだけで、バグの探索範囲がモジュール単位に絞れる恩恵は、計り知れないものがありますよ。
さあ、まずは一つの定数モジュールを、`PRIVATE`の鍵で守ることから始めてみましょう！
コメント