hetrf#
複素エルミート行列のバンチ・カウフマン分解を計算します。このルーチンは、oneapi::mkl::lapack 名前空間に属します。
説明#
このルーチンは、バンチ・カウフマン対角ピボット法を使用して複素エルミート行列 A の因数分解を計算します。
uplo='U'の場合、A = U*D*UHuplo='L'の場合、A = L*D*LH,
ここで、A は入力行列、U と L は順列行列と単位対角の三角行列の積 (U の場合は上三角行列、L の場合は下三角行列)、D は 1 行 1 列と 2 行 2 列の対角ブロックを持つエルミートブロック対角行列です。U と L には、D の 2 行 2 列のブロックに対応する 2 行 2 列の単位対角ブロックがあります。
注
このルーチンは、進捗ルーチン機能をサポートします。
API#
構文#
namespace oneapi::mkl::lapack {
void hetrf(sycl::queue &queue,
mkl::uplo uplo,
int64_t n,
sycl::buffer<T> &a,
int64_t lda,
sycl::buffer<T> &ipiv,
sycl::buffer<T> &scratchpad, int64_t scratchpad_size)
}hetrf は次の精度とデバイスをサポートしています。
T |
デバイスサポート |
|---|---|
|
CPU |
|
CPU |
入力パラメーター#
- queue
この計算が実行されるデバイスキュー。
- uplo
Aの上三角部分と下三角部分のどちらが格納されているか、またAがどのように因数分解されるかを示します。uplo = uplo::upperの場合、配列にはAの上三角部分が格納され、AはU*D*UH として因数分解されます。uplo = uplo::lowerの場合、配列には行列Aの下三角部分が格納され、AはL*D*LH として因数分解されます。- n
行列
A(0 ≤ n)の次数。- a
行列
Aの係数を保持するバッファー。サイズはmax(1,lda*n)で、行列Aの上三角部分または下三角部分が含まれます (uplo を参照)。a の 2 番目の次元はmax(1,n)以上でなければなりません。- lda
a の先頭次元。
- scratchpad
ルーチンが中間結果を維持するスクラッチパッド・メモリーを保持するバッファー。
- scratchpad_size
Tタイプの浮動小数点要素の数であるスクラッチパッド・メモリーのサイズ。サイズは、hetrf_scratchpad_size 関数が返す値以下であってはなりません。
出力パラメーター#
- a
aの上三角部分または下三角部分は、ブロック対角行列Dの詳細と、係数U(またはL) を取得する乗数によって上書きされます。- ipiv
max(1, n)以上のサイズの配列を保持するバッファー。Dのインターチェンジとブロック構造の詳細が含まれます。ipiv(i) = k >0の場合、diiは 1 行 1 列のブロックとなり、Aのi行目と列がk行目と列と交換されます。uplo= mkl::uplo::upper かつipiv(i) =ipiv(i-1) = -m< 0 の場合、Dにはi行目とi-1列目に 2 行 2 列のブロックがあり、Aの (i-1) 行目と列はm行目と列と交換されます。uplo= mkl::uplo::lower かつipiv(i) =ipiv(i+1) = -m< 0 の場合、Dにはi行目とi+1列目に 2 行 2 列のブロックがあり、Aの (i+1) 行目と列はm行目と列と交換されます。
例外#
例外 |
説明 |
|---|---|
|
この例外は、計算中に問題が発生した場合にスローされます。例外オブジェクトの info() メソッドを使用して、問題の情報コードを取得できます。
|