getrs

getrs#

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

説明#

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

方程式

状態

A*X = B

trans=mkl::transpose::nontrans の場合

AT*X = B

trans=mkl::transpose::trans の場合

AH*X = B

trans=mkl::transpose::conjtrans の場合

このルーチンを呼び出す前に、getrf を呼び出して A の LU 分解を計算します。

API#

構文#

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

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

T

デバイスサポート

float

CPU および GPU

double

CPU および GPU

std::complex<float>

CPU および GPU

std::complex<double>

CPU および GPU

入力パラメーター#

queue

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

trans

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

trans=mkl::transpose::nontrans の場合、X について A*X = B が解かれます。

trans=mkl::transpose::trans の場合、X について AT*X = B が解かれます。

trans=mkl::transpose::conjtrans の場合、X について AH*X = B が解かれます。

n

行列 A の次数と行列 B(0 ≤ n) の行数。

nrhs

右辺の数 (0 ≤ nrhs)。

a

getrf によって返される、行列 A の因数分解を含む配列を保持するバッファー。a の 2 番目の次元は max(1, n) 以上でなければなりません。

lda

a の先頭次元。

ipiv

getrf によって返される、max(1, n) 以上のサイズの配列を保持するバッファー。

b

配列 b には、連立方程式の右辺の列を持つ行列 B が含まれます。b の 2 番目の次元は max(1,nrhs) 以上でなければなりません。

ldb

b の先頭次元。

scratchpad

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

scratchpad_size

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

出力パラメーター#

b

バッファー b は、解行列 X によって上書きされます。

例外#

例外

説明

mkl::lapack::exception

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

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

info = i の場合、U の i 番目の対角要素はゼロであり、解決は完了しませんでした。

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