直接法 † 大抵の行列に対して使用可。密な行列に対して使う。疎行列、帯行列のようなものに対しては反復法を利用した方が計算効率・精度の面でよい。 ガウスの消去法(部分ピボット付) LU分解 LU分解(部分ピボット付) LAPACK95を使った方法 Driver routineを使う la_gesv(EXAMPLES1にla_gesv_exampleがあるのでそちらを参照) Computational routineを使う la_getrf でLU分解し、la_getrsで解を求める。(la_getrsは後退代入、前進代入を行うroutine) ex1 ex2 ※例題は「計算による線形代数(共立出版)第5章練習問題」より ↑ 反復法 † ※対角優位行列( , ) に対して解が収束することが保証されているが、他はダメ。 ※安定性、収束の速さ、適用範囲の広さでBiCGStab法が全て