大規模行列の計算

非零要素の位置だけを記録して解くのが大規模疎行列を扱う場合の基本です。 計算の過程で発生する非零要素（専門的にはfill-inと呼びます）の発生位置は予め計算可能なので、その分のエリアは予め空けておきます。 ウエーブ・フロント法は確かfill-inの発生をできるだけ抑制するための変換手法ですよね？ この変換手法を「オーダリング」と呼びますが、実は解く問題によってオーダリングの手法は善し悪しがあります。 ガウスであればより対角に近い場所に非零要素とfill-inが移動するすようにします。 んで、要素のうまい記憶方法ですが、ＬＤＵ分解を用いる解法であれば、 Ｄ：対角要素 d[10000] Ｕ：上△要素 U[??] Ｌ：下△要素 L[??] d0, U1, U2, ... Un L1, d1, Un+1, ... L2,Ln+1,d2, ... : : Ln こんな感じに記憶します。（分かりにく