問題固有のロジックを実行できるように、(スパース A の場合) A\b がいつ失敗したかをプログラムで調べたいと考えています。バックスラッシュ演算子の使用
A\b
コンソールに警告が出力されますが、これらの状態 (特異またはほぼ特異) についてプログラムで知りたいので、問題固有の処理を行うことができます。
高密度システムの場合、私はできる
[soln, cond_recip] = linsolve(A,b);
if cond_recip < 1e-15, ..., end
しかし、linsolve は疎行列では機能せず、行列を高密度化したくありません。