zpbcon zpbtrf zpbtrs zptcon zpttrf zpttrs

zpbcon() void zpbcon ( char  uplo, int  n, int  kd, int  ldab, doublecomplex  ab[], double  anorm, double *  rcond, doublecomplex  work[], double  rwork[], int *  info  ) 行列の条件数 (正定値エルミート帯行列) 目的 本ルーチンはzpbtrfにより計算されたコレスキー分解 A = U^H*U または A = L*L^H を用いて正定値エルミート帯行列の(1-ノルムによる)条件数の推定を行う. norm(inv(A))の推定値を求め, 条件数の逆数を次のように計算する. rcond = 1 / (anorm * norm(inv(A))) 引数 [in] uplo = 'U': 上三角行列Uをab[][]に格納. = 'L': 下三角行列Lをab[][]に格納. [in] n 行列Aの行および列数. (n >= 0) (n = 0 の場合, rcond = 1 を返す) [in] kd 上帯幅(uplo = 'U'の場合)あるいは下帯幅(uplo = 'L'の場合). (kd >= 0) [in] ldab 二次元配列ab[][]の整合寸法. (ldab >= kd + 1) [in] ab[][] 配列 ab[lab][ldab] (lab >= n) 正定値エルミート帯行列Aのコレスキー分解 A = U^H*U または A = L*L^Hの三角行列 U または L (kd+1×n対称帯行列形式). [in] anorm 正定値エルミート帯行列Aの1ノルム(または無限大ノルム). [out] rcond 行列Aの条件数の逆数. 次のように計算する.   rcond = 1/(anorm * ainvnm) ここで, ainvnmは本ルーチン内で求められたinv(A)の1-ノルムの推定値である. [out] work[] 配列 work[lwork] (lwork >= 2*n) 作業領域. [out] rwork[] 配列 rwork[lrwork] (lrwork >= n) 作業領域. [out] info = 0: 正常終了 = -1: 入力パラメータ uplo の誤り (uplo != 'U'および'L') = -2: 入力パラメータ n の誤り (n < 0) = -3: 入力パラメータ kd の誤り (kd < 0) = -4: 入力パラメータ ldab の誤り (ldab < kd+1) = -6: 入力パラメータ anorm の誤り (anorm < 0) 出典 LAPACK