getrsnp_batch (バッファー・ストライド・バージョン)

getrsnp_batch (バッファー・ストライド・バージョン)#

複数の右辺を持つ LU 分解された正方係数行列の線形方程式系を解きます。このルーチンは、oneapi::mkl::lapack 名前空間に属します。

説明#

このルーチンは、Xi について次の線形方程式を解きます。

  • trans = mkl::transpose::notrans の場合、Ai * Xi = Bi

  • trans = mkl::transpose::trans の場合、AiT * Xi = Bi

  • trans = mkl::transpose::conjtrans の場合、AiH * Xi = Bi

このルーチンを呼び出す前に、getrfnp_batch (USM バージョン) (ピボットなし) を呼び出して Ai の LU 分解を計算します。

API#

構文#

 namespace oneapi::mkl::lapack { 
  void getrsnp_batch(sycl::queue &queue, 
  mkl::transpose trans, 
  int64_t n, 
  int64_t nrhs, 
  sycl::buffer<T> &a, 
  int64_t lda, 
  int64_t stride_a, 
  sycl::buffer<T> &b, 
  int64_t ldb, 
  int64_t stride_b, 
  int64_t batch_size, 
  sycl::buffer<T> &scratchpad, 
  int64_t scratchpad_size) 
}

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

T

デバイスサポート

float

CPU および GPU

double

CPU および GPU

std::complex<float>

CPU および GPU

std::complex<double>

CPU および GPU

入力パラメーター#

queue

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

trans

方程式の形式を示します。

trans = mkl::transpose::nontrans の場合、Xi について Ai*Xi = Bi が解かれます。

trans = mkl::transpose::trans の場合、Xi について AiT*Xi = Bi が解かれます。

trans = mkl::transpose::conjtrans の場合、Xi について AiH*Xi = Bi が解かれます。

n

行列 Ai の次数と行列 Bi の行数 (n ≥ 0)。

nrhs

右辺の数 (nrhs ≥ 0)。

a

getrfnp_batch (バッファー・ストライド・バージョン) によって返される、行列 Ai の因数分解を含む配列。

lda

Ai の先頭次元 i(lda ≥ max(1, n))。

stride_a

バッチ配列 a 内の行列 Ai の先頭間のストライド (stride_a ≥ max(1, lda * n))。

b

連立方程式の右辺の列を持つ行列 Bi を含む配列。

ldb

Bi の先頭次元 i(ldb ≥ max(1, n))。

stride_b

バッチ配列 b 内の行列 Bi の先頭間のストライド (stride_b ≥ max(1, ldb * nrhs))。

batch_size

バッチ内の問題数 (batch_size ≥ 0)。

scratchpad

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

scratchpad_size

T タイプの浮動小数点要素数のスクラッチパッド・メモリーのサイズ。サイズは、getrsnp_batch_scratchpad_size (ストライドバージョン) が返す値よりも小さくてはなりません。

出力パラメーター#

b

解行列 Xi によって上書きされます。

例外#

例外

説明

mkl::lapack::batch_exception

この例外は、計算中に問題が発生した場合にスローされます。例外オブジェクトの info() メソッドを使用して、問題の情報コードを取得できます。

info = -n の場合、n 番目のパラメーターの値が不正です。

info がスクラッチパッドのサイズとして渡された値と等しく、detail() がゼロ以外を返す場合、渡されたスクラッチパッドのサイズが十分ではなく、必要なサイズは例外オブジェクトの detail() メソッドによって返される値以上である必要があります。

info がゼロの場合、Ui の一部の対角要素がゼロになり、解決を完了できません。バッチ内のこのような行列のインデックスは、例外オブジェクトの ids() メソッドを使用して取得できます。例外オブジェクトの infos() メソッドを使用して、これらの Ui 行列の最初のゼロ対角要素のインデックスを取得できます。