「なぜFortranは行列計算が速いのか？」メモリ配置の深淵を覗く
数値計算の世界へようこそ。C言語やPythonの経験がある皆さんなら、「Fortranはなぜ古いのに、スパコンのトップシェアを維持しているのか？」と不思議に思ったことがあるかもしれません。
その答えの大部分は、「メモリの扱い方」への徹底したこだわりにあります。今日は、Fortranの心臓部とも言える「列優先（Column-Major Order）」という概念について、現場で血を流した経験を交えながらお話ししましょう。
—
1. 「列優先」とは何か？——図書館の本棚に例えてみよう
C言語やPython（NumPy）のユーザーは、メモリ上にデータが「行」の順に並んでいる世界（行優先）に慣れています。しかし、Fortranは違います。
Fortranは、「縦の列」がメモリ上で隣り合っているというルールを採用しています。
これを図書館の本棚で例えてみましょう。
- 行優先（C/Python）： 本を左から右へ、1冊ずつ並べる。
- 列優先（Fortran）： 本を「上から下へ」積み上げ、それが終わったら次の列へ移動する。
Fortranで `A(i, j)` と書いたとき、`i`（行）をインデックスとして回すことは、メモリ上では「隣り合った場所」を次々と参照することを意味します。これが、CPUのキャッシュメモリと相性が抜群に良い理由なのです。
—
2. 最適化の鉄則：ループの回し方ひとつで世界が変わる
では、具体的に見てみましょう。以下のコードは、行列の全要素を合計するだけの単純な例です。
! 悪い例：列優先の原則を無視している
do j = 1, n
do i = 1, m
sum = sum + A(i, j) ! ここが重要！
end do
end do
なぜこの書き方が「正解」なのか。もしこれが `A(j, i)` だった場合、メモリ上では大きく離れた場所へ次々とジャンプすることになります。CPUは「次に使うデータ」を事前にキャッシュに読み込んでおこうとしますが、飛び飛びのアクセスだとキャッシュが役に立たず、メモリへわざわざ値を取りに行く「キャッシュミス」という高コストな現象が発生します。
現場の教訓：
「一番左のインデックス（`i`）を、一番内側のループで回せ」。これがFortranで計算速度を稼ぐための黄金律です。
—
3. 実践！キャッシュを味方につけるコード設計
単にコードを書くだけではなく、コンパイラに「ここを最適化してくれ！」と背中を押すことも大切です。以下のコードは、モダンFortranで配列演算を行う際の推奨パターンです。
program matrix_sum
implicit none
integer, parameter :: m = 1000, n = 1000
real(8), allocatable :: A(:, 🙂
real(8) :: total_sum
integer :: i, j
allocate(A(m, n))
A = 1.0d0 ! 配列全体に値を代入（Fortranの得意技）
total_sum = 0.0d0
! 【最適化のポイント】
! 内側のループで第一添字(i)を動かすことで、
! メモリ上の連続アクセスを保証する
do j = 1, n
do i = 1, m
total_sum = total_sum + A(i, j)
end do
end do
print , “合計値: “, total_sum
deallocate(A)
end program matrix_sum
—
4. コンパイラへの「魔法の呪文」
モダンな開発環境では、コードの書き方だけでなく、コンパイルオプション（フラグ）も重要です。例えば `gfortran` を使う場合、以下のようなフラグを検討してみてください。
現場でのビルド例
gfortran -O3 -march=native -ffast-math -o solver main.f90
- `-O3`: 最も攻撃的な最適化を指示します。
- `-march=native`: あなたのCPUの機能をフル活用する設定です。
- `-ffast-math`: 厳密な数学的順序よりも速度を優先します（数値的安定性に注意が必要ですが、爆速になります）。
—
最後に：なぜ今、Fortranなのか
皆さんが他言語から持ち込んだ「データ構造の考え方」は素晴らしい武器です。しかし、Fortranの世界では「データはメモリにどう並んでいるか？」を想像するだけで、計算速度が10倍、100倍と変わることがあります。
最初は戸惑うかもしれませんが、「第一添字を内側で回す」というたった一つのルールを守るだけで、あなたのプログラムはスパコンの性能を最大限に引き出す準備が整います。
ぜひ、自分の書いたコードがメモリの上をどう走っているか、脳内でシミュレーションしてみてください。そこには、数値計算のプロだけが知る「静かなる興奮」が待っていますよ！
何か詰まったら、いつでも聞いてくださいね。応援しています！
コメント