BLAS,LAPACKの使い方BLAS,LAPACKの使い方 LAPACK(Linear Algebra PACKages)はC/Fortran用の高速線形代数演算パッケージで、行列とベクトル(行列)の積、行列方程式(Ax=b)の数値解法などを高速に行います。 例1:連立一次方程式の解法 以下の方程式の数値解を求めるとする。 x-y+2w = 0 y+2z+3w = 2 2x-z+3w = 1 2x-3y+z+4w = 1 この方程式をLapackで解くソースファイル(gausselim.c)はここ。 ※(注意点その1) 行列も１次元配列で定義することと、列ごとに定義していくこと(一行一列目の次は二行一列目を定義する。CとFortranとの二次元配列のアドレスの違いによる) ※(注意点その2) コンパイル時は > gcc -llapack -lblas -lg2c gausselim.c
