連立一次方程式を解く

遅くなりましたが、LAPACKを使用してみました。
恐ろしく計算速度が速いですね。
環境や、コードの最適化がかなりされているようですね。
1000×1000の係数行列を数秒で計算していました。
ちなみにLinux PIII 500MHz 512MBです。

で、係数行列を増やしたときの計算速度を見てみましたが、
7000×7000ぐらいからスワップを使用し始めたためか、
極端に速度が落ちましたね。

しかし、環境毎にLAPACKをインストールする必要があるため、
プログラムの移植が面倒ですね。
で、ライブラリ(***.a)をコンパイル時にリンクするという方法を取ったら、上手く動作しました。
こんな方法でもいいのかな？

LAPACKの移植性の悪さから見て、TNTと言うのは良いですね。
ヘッダーをインクルードするだけでいいんですかね。
しかも、LAPACKの作成者と同じ人が作っている？
けど、処理速度はLAPACKの方が速いみたいですね。
昔から最適化されつづけた差でしょうか？
ちなみにTNT試しましたが、コンパイルできませんでした。
原因がよくわからないので、ver1.00のリリースを待つかな。

そこで、私は、このLAPACKに勝てるものは無いような気がするので、このLAPACKを使用しようと思います。
ちょっと見てみるとバンドマトリクスも解くことが出来るようなので、それを利用しようと思います。
しかし、ドキュメントを見てみると、引数のサブ対角線？スーパー対角線？という意味がわかんない。
これは、バンドマトリクスを理解する必要がありますね。
自分なりに理解してみるつもりですが、分からなかったらまた質問しちゃうかも。ちょっと、他力本願すぎますかな？

後本ですが、古いバージョンの日本語訳があるようなので、余裕があれば見てみようと思います。

ありがとうございます。
コレスキー法は対称行列についてかなり力を発揮する解法のようですね。

今回は使用しないかもしれませんが、対称連立一次方程式を扱うときは是非使用しようと思います。