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

m 行 n 列の一般行列バッチの LU 分解バッチを計算します。このルーチンは、oneapi::mkl::lapack 名前空間に属します。

説明

このルーチンは、m 行 n 列の一般行列 A_i の LU 分解を次のように計算します。

\[A_i = P_i \times L_i \times U_i\]

ここで、P_i は順列行列、L_i は単位対角要素を持つ下三角行列 (m > n の場合は下台形行列) であり、U_i は上三角行列 m < n の場合は上台形行列) です。このルーチンは、行の交換を伴う部分ピボットを使用します。

API

構文

namespace oneapi::mkl::lapack { 
  void getrf_batch( sycl::queue &queue, 
  int64_t m, 
  int64_t n, 
  sycl::buffer<T> &a, 
  int64_t lda, 
  int64_t stride_a, 
  sycl::buffer<int64_t> &ipiv, 
  int64_t stride_ipiv, 
  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

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

m

行列 A_i の行数 (m ≥ 0)。

n

行列 A_i の列数 (n ≥ 0)。

a

入力行列 A_i を保持する配列。

lda

A_i の先頭次元 i(lda ≥ max(1, m))。

stride_a

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

stride_ipiv

配列 ipiv 内の行列 ipiv_i の先頭間のストライド (stride_ipiv ≥ max(1, min(m,n)))。

batch_size

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

scratchpad

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

scratchpad_size

T タイプの浮動小数点要素の数であるスクラッチパッド・メモリーのサイズ。サイズは、getrf_batch_scratchpad_size (ストライドバージョン) によって返される値以下であってはなりません。

出力パラメーター

a

L_i と U_i によって上書きされます。L_i の単位対角要素は保存されません。

ipiv

max(1,min(m,n)) 以上のサイズのピボット・インデックス ipiv_i のバッチを含む配列。1 ≤ k ≤ min(m,n) の場合、A_i の行 k は行 ipiv_i(k) と交換されます。

例外

例外	説明
mkl::lapack::batch_exception	この例外は、計算中に問題が発生した場合にスローされます。例外オブジェクトの info() メソッドを使用して、問題の情報コードを取得できます。 info = -n の場合、n 番目のパラメーターの値が不正です。 info がスクラッチパッドのサイズとして渡された値と等しく、detail() がゼロ以外を返す場合、渡されたスクラッチパッドのサイズが十分ではなく、必要なサイズは例外オブジェクトの detail() メソッドによって返される値以上である必要があります。 info がゼロの場合、Ui の一部の対角要素がゼロになり、解決を完了できません。バッチ内のこのような行列のインデックスは、例外オブジェクトの ids() メソッドを使用して取得できます。例外オブジェクトの infos() メソッドを使用して、これらの Ui 行列の最初のゼロ対角要素のインデックスを取得できます。