这篇文章关注的是非常短,很宽的数组(#列可能比行数大几个数量级)。
由于行数/列数的差异以及我使用的矩阵的大尺寸,将LU分解的U部分保留在内存中通常是不可行的。 Eigen是否具有仅计算L的功能?等价地,使用行操作将输入矩阵放置为梯形形式?
一般说明
(1)我在这里看到了一个相关问题
https://forum.kde.org/viewtopic.php?f=74&t=138686&p=371097&hilit=echelon#p371097
答案建议查看FullPivLU下的image()方法,但我无法在文档中找到必要的信息。特别地,实际上,获得矩阵L通常很重要。矩阵的列空间的任意基础是不够的。
(2)这里还有一个问题
https://forum.kde.org/viewtopic.php?f=74&t=130430&p=348923&hilit=echelon#p348923
但它似乎没有得到回应。
(3)在引起这个问题的(相当专业的)应用程序领域中,稳定性问题较少受到关注,因为我们通常在有限域中工作。
谢谢!