chkginv_JPH {lazy.symbolic}R Documentation

第二の行列が第一の行列の一般化逆行列であるか、またはその逆かをチェック

Description

英語のヘルプファイルは: chkginv

Arguments

A

行列

Ainv

A の一般化逆行列である可能性のある行列

eps

ゼロ判定のための許容誤差

print

= 1 の場合、結果を出力

Details

以下の4つの条件が満たされているかを確認します。 g1) A - A g2) Ainv - Ainv g3) A g4) Ainv さらに、入力行列が正方行列の場合は以下の2条件も追加されます。 g01) I - A g02) I - Ainv

各条件について、絶対値の最大値を eps と比較し、ゼロと見なせるかを判定します。 デフォルト値は max(abs(c(diag(Ainv), diag(A)))) * sqrt(.Machine$double.eps) です。

Value

以下のリストを返します: check=list(g1,g2,g3,g4,g01,g02) 各 g1, g2,... は OK または NG maxad=list(g1,g2,g3,g4,g01,g02) 各 g1, g2,... は条件ごとの絶対値の最大値

Examples

chkginv(demomat(3,3), matSwp(demomat(3,3)))

# Ginv や matSwp から呼び出す場合
A=demomat(9,4)/100; A=t(A)%*%A; Ainv=Ginv(A,chkginv=2)


[Package lazy.symbolic version 1.0.0.20250830 ]