2011-08-08 22 views
3

私はBLAS/Lapackの新しいユーザーです。ガウス消去やガウス・ジョーダン消去を行うルーチンがあるのでしょうか?私はグーグルで調べて、彼らの文書を見ましたが、まだそれらを見つけることができませんでした。BLAS/LAPACKガウス消去のルーチン

ありがとうございました!

+0

Gaussian/Gauss-Jordanでなければならないのですか、それとも線形方程式の解を求めるアルゴリズムであればよいですか? –

+0

@Matt Ball私はガウス/ガウスジョーダンを好むかもしれないと思うのは、得られたマトリックスがエシェロン形式になるからです。他のアルゴリズムでも行のエシェロン形式が生成されますか? –

答えて

6

ガウス消去は基本的にLU分解と同じです。ルーチンxGETRFはLU分解を計算する(例えば、実際の倍精度行列についてはDGETRF)。 U因子は、ガウス消去後の行列に対応する。 U係数は、出口の行列Aの上三角部分(対角を含む)に格納されます。

LU分解/ガウス消去は、一般に線形連立方程式を解くために使用されます。 LU分解を計算したら、xGETRSルーチンを使用して線形システムを解くことができます。