LAPACK

LAPACK (Linear Algebra PACKage) 数値線形代数のための数値解析ソフトウェアライブラリで、線型方程式や線型最小二乗問題、固有値問題、特異値問題等を数値的に解くために利用される。本ライブラリは複素数または実数を成分とする行列を扱うことが可能であり、LU分解やコレスキー分解、QR分解、シュア分解等の行列の分解を行うためのサブルーチンを含む。サブルーチンは単精度版と倍精度版が提供される。1992年 (1992)のLAPACKの初版はFORTRAN 77 で実装されていたが、現在はFortran 90が用いられている。LAPACK 3.4.0からはC言語インターフェースであるLAPACKEが統合され、C言語やC++からの利用が容易になった。 LAPACKはLINPACKおよびEISPACKの後継と見做されている。ただし、LINPACKの設計が開発当時近代的であった共有メモリ型ベクトルコンピュータを意識したものであるのに対して、本ライブラリの設計はキャッシュを用いたアーキテクチャを有する、より近代的なコンピュータを意識したものである。LAPACKはLINPACK同様にBLAS（Basic Linear Algebra Subprograms、基本線型代数サブプログラム群）ライブラリ上に構築されている。LAPACKは後に分散メモリ型のコンピュータ向けにScaLAPACK（英語版）やPLAPACK（英語版）へと拡張された。 LAPACKはBSDライセンスで提供されるオープンソースソフトウェアである。 LAPACKのサブルーチンは以下の三種類に大別される。 ドライバルーチン（driver routines） LAPACKが扱うことが可能な問題を解くためのルーチン。 問題の例として線型方程式系を解く問題や対称行列の固有値問題などが挙げられる。利用者の要請に合致する機能のドライバルーチンが存在する場合にはそのルーチンの利用が推奨される。 計算ルーチン（computational routines） 問題を解くために必要な計算タスクを実行するためのルーチン。LAPACKのドライバルーチンは計算ルーチンを連続的に呼び出すことで問題を解く。計算タスクの例として行列をLU分解することや対称行列を三重対角行列に変換することなどが挙げられる。前者は線型方程式系を解くために、そして後者は対称行列の固有値問題を解くために必要である。利用者の要請に合致するドライバルーチンが存在しない場合は計算ルーチンを組み合わせて問題を解くことになる。 補助ルーチン（auxiliary routines） 補助的に利用されるルーチン。ブロックアルゴリズム内部で利用される計算タスクの一部を実行するものや、BLASの機能をわずかに拡張したものが含まれる。 LAPACKとBLASのサブルーチンの名称は機能の判別が平易である範囲で短くなるような規則で命名されている。 これは初期のFORTRANにおける関数の名称に関する仕様上の制限を受けたものである。 サブルーチンはpmmaaaの規則で命名される。 以下、LAPACKのDGESV（倍精度一般行列の方程式系の求解）とBLASのDGEMM（倍精度行列の積の計算）を例に挙げる。 pは通常は一文字の英字で数値データの型を表現するために利用される。SとDはそれぞれ単精度と倍精度の実数を意味し、CとZはそれぞれ単精度と倍精度の複素数を意味する。なお、一文字目を別の文字に置き換えてサブルーチンの名称を記述することがある。ただし、サブルーチンが採用するアルゴリズムによっては例外が存在することに注意が必要である。 mmは二文字の英字でサブルーチンが採用するアルゴリズムが想定する行列の型を意味する。行列の型を現す略号は以下に表で示す。サブルーチンが想定する行列データの格納方法は型ごとに異なる。例えば、対角行列を意味するDIが与えられた場合は対角要素が格納された長さnの配列を、一般行列を意味するGEが与えられた場合は行列の要素が格納されたn×nの配列という具合いである。 aaaは一文字から三文字の英数字でサブルーチンの処理内容を表現する。例えばSVは線型方程式系の求解を意味し、MMは行列の積を意味する。 命名規則の詳細はLAPACK Users' Guideの当該項目を参照。 LAPACKにはサブルーチンの機能とインターフェースに互換性のある実装が多く存在する。 以下に例を示す。それぞれの動作環境についてはリンク先を参照。 ACML (AMD Core Math Library): AMDによる実装。 Accelerateフレームワーク: Appleによる実装。 HP MLIB: ヒューレット・パッカードによる実装。 ESSL (Engineering and Scientific Subroutine Library): IBMによる実装。 Intel Math Kernel Library: インテルによる実装。 MathKeisan: 日本電気による実装。 SCSL (Scientific Computing Software Library): SGIによる実装。 Sun Performance Linaray: サン・マイクロシステムズによる実装。 LAPACKはFortran 90以外のプログラミング言語から利用することが可能であり、これを目的とした言語バインディングのためのライブラリも開発されている。LAPACK 3.4.0よりC言語インターフェースであるLAPACKEが統合された。 以下に例を示す。 LAPACK95: Fortran 95用。インターフェースが簡略化される。 CLAPACK: C言語用。Fortran版をf2c（英語版）によって変換したもの。 LAPACKE: C言語用。インテルのMath Kernel Library互換のC言語インターフェース。 LAPACK++: C++用。 Armadillo: C++用。 CPPLAPACK: C++用。 Boost Numeric Bindings: C++用。 SciPy: Python用。 jlapack: Java用。 CSLapack: C#用。 Linalg: Ruby用。 LACAML: OCaml用。 hmatrix: Haskell用。 Accelerateフレームワーク: Objective-C用インターフェースが提供される。 Gonum: LAPACK and BLAS in Native Go 多くの処理系ではライブラリのC言語バインディングが可能であるため、いくらかの制約があるにせよLAPACKのサブルーチンをC言語の関数のように利用できる。ただし、Fortranコンパイラが存在しない処理系ではCLAPACKも有効な選択肢となる。なお、Automatically Tuned Linear Algebra Software（英語版）が提供するC言語インターフェースも存在するが、これはf2cのCLAPACKとは互換性が無い。 C言語では慣用のrow-major orderingな行列を計算させる場合、LAPACKEでは指定オプションが用意されている。その場合、内部で行列転置が行われオーバヘッドが発生する。例えば行列積は、内部の計算は下記のように行われる。ここで「 ⋅ {\displaystyle \cdot } 」はBLAS/LAPACK本来のcolumn-major orderingな行列積を意味している。 C = ( A T ⋅ B T ) T {\displaystyle C=\left(A^{\textrm {T}}\cdot B^{\textrm {T}}\right)^{\textrm {T}}}