Matlabを使用して最小ピボットだけのLU分解を見つける可能性はありますか？

Question

私は、非線形ソルバーの各繰り返しで方程式の線形システム（Ax = b）を解こうとしています。 Lに単位対角要素がある場合、Uの最小ピボットに対してのみ、Aの不完全LUに関する情報が必要です。私のプログラムはMatlabで書かれています。

私は完全なLUを実行し、最小のピボットに対して対応するLおよびUベクトルを見つけることができますが、行列の次元が非常に大きいのでコストがかかります。

Answer 1

thisコードが役立つと思います。私はあなたが探しているものは完全にはわかりませんが、Aの各列の最小値をピボットとして使用する場合は、上記のコードの13行目のmaxをminに変更してください。