potri

potri#

コレスキー分解を使用して対称 (エルミート) 正定値行列の逆行列を計算します。このルーチンは、oneapi::mkl::lapack 名前空間に属します。

説明#

このルーチンは、対称正定値行列、または複素数の場合はエルミート正定値行列 A の逆行列 inv(A) を計算します。このルーチンを呼び出す前に、potrf を呼び出して A を因数分解します。

API#

構文#

namespace oneapi::mkl::lapack { 
  void potri(sycl::queue &queue, 
  mkl::uplo uplo, 
  int64_t n, 
  sycl::buffer<T> &a, 
  int64_t lda, 
  sycl::buffer<T> &scratchpad, 
  int64_t scratchpad_size) 
}

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

T

サポートされるデバイス

float

CPU および GPU

double

CPU および GPU

std::complex<float>

CPU および GPU

std::complex<double>

CPU および GPU

入力パラメーター#

queue

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

uplo

入力行列 A がどのように因数分解されたかを示します。

uplo = mkl::uplo::upper の場合、A の上三角が格納されます。

uplo = mkl::uplo::lower の場合、A の下三角が格納されます。

n

行列 A の次数を指定します (0 ≤ n)。

a

potrf (USM バージョン) によって返される行列 A の因数分解が含まれます。a の 2 番目の次元は max(1, n) 以上でなければなりません。

lda

a の先頭次元 。

scratchpad

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

scratchpad_size

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

出力パラメーター#

a

A の逆数の上三角形または下三角形によって上書きされます。uplo で指定されます。

例外#

例外

説明

mkl::lapack::exception

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

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

info = i の場合、コレスキー因子の i 番目の対角要素 (つまり因子自体) はゼロであり、反転は完了できません。

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