ormrq#

gerqf によって形成された RQ 分解の直交行列 Q を実行列に乗算します。

説明 #

このルーチンは、実数 m 行 n 列の行列 C を Q または Q^T で乗算します。ここで、Q は、k 個の基本鏡映子 H_i の積として定義される実直交行列です: Q = H1H2 ... Hk は RQ 因数分解ルーチン gerqf によって返されます。

パラメーター side と trans に応じて、ルーチンは行列積 Q*C、QT*C、C*Q、または C*QT (結果を C に上書き) のいずれかを形成できます。

API #

構文#

namespace oneapi::mkl::lapack { 
  void ormrq(sycl::queue &queue, 
  mkl::side side, 
  mkl::transpose trans, 
  int64_t m, 
  int64_t n, 
  int64_t k, 
  sycl::buffer<T> &a, 
  int64_t lda, 
  sycl::buffer<T> &tau, 
  sycl::buffer<T> &c, 
  int64_t ldc, 
  sycl::buffer<T> &scratchpad, 
  int64_t scratchpad_size) 
}

ormrq は次の精度とデバイスをサポートしています。

T	デバイスサポート
`float`	CPU
`double`	CPU

入力パラメーター#

queue

計算が実行されるデバイスキュー。

side

side = mkl::side::left の場合、Q または Q^T が左から C に適用されます。side = mkl::side::right の場合、Q または Q^T が右から C に適用されます。

trans

trans=mkl::transpose::trans の場合、ルーチンは C に Q を乗算します。

trans=mkl::transpose::nontrans の場合、ルーチンは C に Q^T を乗算しますT。

m

行列 A の行数 (0 ≤ m)。

n

行列 A の列数 (0 ≤ n ≤ m)。

k

積が行列 Q を定義する基本反射体の数 (0 ≤ k ≤ n)。

a

gerqf 関数の結果を保持するバッファー。a の 2 番目の次元は max(1,k) 以上でなければなりません。

lda

a の先頭次元。

tau

gerqf 関数によって返された tau を保持するバッファー。

c

行列 C を保持するバッファー。c の 2 番目の次元は max(1,n) 以上でなければなりません。

ldc

c の先頭次元。

scratchpad

ルーチンが中間結果を維持するスクラッチパッド・メモリーを保持するバッファー。

scratchpad_size

T タイプの浮動小数点要素の数であるスクラッチパッド・メモリーのサイズ。サイズは、ormrq_scratchpad_size 関数が返す値以下であってはなりません。

出力パラメーター#

c: 積 Q*C、Q^T*C、C*Q、または C*Q^T (side および trans で指定) によって上書きされます。

例外#

例外	説明
`mkl::lapack::exception`	この例外は、計算中に問題が発生した場合にスローされます。例外オブジェクトの info() メソッドを使用して、問題の情報コードを取得できます: `info = -i` の場合、`i` 番目のパラメーターの値が不正です。`info` がスクラッチパッドのサイズとして渡された値と等しく、detail() がゼロ以外を返す場合、渡されたスクラッチパッドのサイズが十分ではなく、必要なサイズは例外オブジェクトの detail() メソッドによって返される値以上である必要があります。

例外

説明

mkl::lapack::exception

この例外は、計算中に問題が発生した場合にスローされます。例外オブジェクトの info() メソッドを使用して、問題の情報コードを取得できます: info = -i の場合、i 番目のパラメーターの値が不正です。info がスクラッチパッドのサイズとして渡された値と等しく、detail() がゼロ以外を返す場合、渡されたスクラッチパッドのサイズが十分ではなく、必要なサイズは例外オブジェクトの detail() メソッドによって返される値以上である必要があります。

アプリケーション・ノート#

このルーチンの複素数対極は unmrq です。

ormrq

目次