hetrd#
複素エルミート行列を三角対角形式に縮小します。このルーチンは、oneapi::mkl::lapack 名前空間に属します。
説明#
このルーチンは、ユニタリー相似変換によって複素エルミート行列 A を対称三角形式 T に縮小します: A = Q*T*QH。ユニタリー行列 Q は明示的には形成されませんが、n-1 個の基本鏡映子の積として表されます。この表現で Q を操作するルーチンが提供されています。
API#
構文#
namespace oneapi::mkl::lapack {
sycl::event hetrd(sycl::queue &queue,
mkl::uplo uplo,
int64_t n,
sycl::buffer<T> &a,
int64_t lda,
sycl::buffer<T> &d,
sycl::buffer<T> &e,
sycl::buffer<T> &tau,
sycl::buffer<T> &scratchpad,
int64_t scratchpad_size)
}hetrd は次の精度とデバイスをサポートしています。
T |
デバイスサポート |
|---|---|
|
CPU および GPU* |
|
CPU および GPU* |
*インターフェイスのサポートのみ。すべての計算は CPU 上で実行されます。
入力パラメーター#
- queue
この計算が実行されるデバイスキュー。
- uplo
uplo::upperまたはuplo::lowerである必要があります。uplo = uplo::upperの場合、a はAの上三角部分を格納します。uplo = uplo::lowerの場合、a はAの下三角部分を格納します。- n
行列
A(0 ≤ n)の順序。- a
行列
Aを保持するサイズ(lda,*)のバッファー。Uplo で指定された上部または下部の三角形が含まれます。- lda
a の先頭次元;
max(1, n)以上- scratchpad
ルーチンが中間結果を維持するスクラッチパッド・メモリーを保持するバッファー。
- scratchpad_size
Tタイプの浮動小数点要素の数であるスクラッチパッド・メモリーのサイズ。サイズは、hetrd_scratchpad_size 関数が返す値以下であってはなりません。
出力パラメーター#
- a
終了時、
uplo = uplo::upperの場合、Aの対角成分と最初の上対角成分は三重対角行列Tの対応する要素で上書きされ、最初の上対角成分より上の要素はバッファー tau とともに、基本鏡映子の積としてユニタリー行列Qを表します;uplo = uplo::lowerの場合、Aの対角成分と最初の副対角成分は三重対角行列Tの対応する要素で上書きされ、最初の副対角成分成分より下の要素は配列 tau とともに、基本鏡映子の積としてユニタリー行列Qを表します。- d
行列
Tの対角要素を保持するバッファー。d の次元はmax(1, n)以上でなければなりません。- e
行列
Tの対角要素を保持するバッファー。e の次元はmax(1, n-1)以上でなければなりません。- tau
max(1, n)以上のサイズの配列を保持するバッファー。ユニタリー行列Qの分解における基本鏡映子を定義する(n-1)個のスカラーをn-1個の基本鏡映子の積に格納します。tau(n)はワークスペースとして使用されます。
例外#
例外 |
説明 |
|---|---|
|
この例外は、計算中に問題が発生した場合にスローされます。例外オブジェクトの info() メソッドを使用して、問題の情報コードを取得できます。
|