有沒有可能找到使用Matlab的最小樞軸LU分解？

Question

我想在非線性求解器的每次迭代中求解一個線性方程組（Ax = b）。我需要關於A的不完全LU的信息，只有U中的最小樞軸，如果L有單位對角線條目。我的程序是用Matlab編寫的。

我可以做完整的LU，並找到相應的L和U向量爲最小的樞軸，但它是昂貴的矩陣尺寸非常大。

Answer 1

我想你可能會發現this代碼很有幫助。我不完全確定你在找什麼，但是如果你只想用A的每一列作爲支點使用最小值，你可以在上面代碼的第13行上將max更改爲min。