継承とポリモーフィズムの罠：数値計算屋がCLASSキーワードと正しく付き合う方法
Fortranの`CLASS`キーワードによるポリモーフィズムは、一見するとオブジェクト指向言語の甘い果実のように見えます。しかし、大規模数値シミュレーションの現場において、何も考えずにこれを乱用することは「性能に対する背信行為」と同義です。
なぜなら、`TYPE`（静的型）がコンパイル時にメモリレイアウトを確定させ、コンパイラに強力な最適化のヒントを与えるのに対し、`CLASS`（動的型）は実行時まで型を確定させないため、コンパイラから「最適化の機会」を奪い去る可能性があるからです。
今回は、数値計算屋が現場で「安全」かつ「高速」にポリモーフィズムを実装するための流儀を伝授します。
—
1. なぜ「ポリモーフィズム」が性能を殺すのか
`CLASS`オブジェクトは、実行時に型情報を参照するための「ディスパッチテーブル（VTable）」を背負っています。ホットループの中で`CLASS`型変数を通じてメソッドを呼び出すと、コンパイラはインライン展開を躊躇し、分岐予測ミスやキャッシュミスを誘発します。
鉄則：ホットループ（計算の核心部）で`CLASS`を使わない。
ポリモーフィズムは、計算の前処理や後処理、あるいは物理モデルの切り替えといった「抽象度が高いレイヤー」に限定し、数億回のループを回す演算カーネル内では、静的な`TYPE`にキャストするか、継承を使わない設計を優先してください。
—
2. 実装例：セキュアかつ高速なインタフェース設計
以下は、流体解析における「境界条件クラス」を想定した設計例です。`CLASS`を使って抽象化しつつ、計算負荷を最小限に抑えるためのポイントを盛り込んでいます。
module m_boundary_base
implicit none
private
public :: boundary_t, apply_boundary
! 抽象基底クラス
type, abstract :: boundary_t
contains
procedure(apply_if), deferred, public :: apply
end type
! インターフェース定義
abstract interface
subroutine apply_if(this, field)
import :: boundary_t
class(boundary_t), intent(in) :: this
real(8), intent(inout) :: field(:,:)
end subroutine
end interface
contains
! ポリモーフィックな呼び出し口
subroutine apply_boundary(obj, field)
class(boundary_t), intent(in) :: obj
real(8), intent(inout) :: field(:,:)
! ここで動的ディスパッチが発生する
call obj%apply(field)
end subroutine
end module
3. 最適化を殺さないための「静的への引き戻し」
もしポリモーフィズムを使って計算の切り替えを行いたい場合、ループの内側で`CLASS`を判定させるのではなく、「ポリモーフィズムをループの外側で行う」のが定石です。
! NG: ループ内でCLASS判定（分岐が重すぎる）
do j = 1, ny
do i = 1, nx
call obj%compute(field(i,j))
end do
end do
! OK: 具象型にキャストしてループを回す（コンパイラがベクトル化できる）
select type(obj)
type is (concrete_a_t)
call solve_kernel(obj, field) ! ここでコンパイラは型を確信し、SIMD展開できる
type is (concrete_b_t)
call solve_kernel(obj, field)
end select
このように`select type`で型を確定させた後に、引数として具象型を渡すことで、コンパイラは初めて「メモリレイアウトが静的である」と認識し、AVX-512等のベクトル命令をフル活用できるようになります。
—
4. シニアエンジニアからの警告：デバッグとビルド設定
`CLASS`を用いた設計において、最も恐ろしいのは「実行時型エラー」です。以下のデバッグフラグは、現場で必ず有効にしてください（gfortran/ifort共通）。
- `-fcheck=all` (gfortran) / `-check all` (ifort):
`select type`の不整合や、未割り当ての`CLASS`ポインタへのアクセスを即座に捕捉します。性能は落ちますが、開発段階では必須です。
- `-O3 -march=native -fno-trapping-math`:
数値計算のパフォーマンスを追求する際、`CLASS`を使うとどうしても最適化が鈍ります。その分、浮動小数点演算の堅牢性を担保しつつ、ベクトル化を阻害する境界チェック等を外す設定を併用してください。
まとめ
1. `CLASS`は「設定」に使う。 物理モデルのスイッチや入出力などの抽象化に留める。
2. 計算の「核心（カーネル）」は`TYPE`に引き戻す。 `select type`で静的に型を確定させ、コンパイラにベクトル化の余地を渡す。
3. データ構造は「列優先（Column-major）」を意識せよ。 `CLASS`で包んだとしても、その内部の配列アクセスがメモリの連続性を守っているかを確認し続ける。
Fortranの真の強みは、オブジェクト指向という「近代的な道具」を使いながら、同時に「ハードウェアの限界まで演算性能を絞り出せる」点にあります。`CLASS`を単なる抽象化の道具としてではなく、ハードウェアへの命令変換を制御するための「スイッチ」として扱うことができれば、貴方のコードは一段上のレベルに達するはずです。
コメント