是否有可能使用Matlab找到最小枢轴的LU分解？

我试图在非线性解算器的每次迭代中解一个线性方程组（Ax=b）。我需要关于A的不完全LU的信息，只针对U中最小的轴，如果L有单位对角线条目。我的程序是用Matlab编写的

我可以进行全LU运算，并为最小的轴找到相应的L和U向量，但由于矩阵维数非常大，因此成本很高。

我认为您可能会发现代码很有用。我不完全确定您要查找的是什么，但如果您只想使用

A

每列中的最小值作为枢轴，您可以在上面代码的第13行将

max

更改为

min