我想知道在matlab中消除(在线性系统求解意义上)矩阵的某些行的最有效方法是什么。基本上,给定行“s”和“l”(这样s和l是不同的,它们的联合是1:n),如何有效地计算
A是一个很大的n x n稀疏SPD矩阵。
R = chol(A(s,s))
Atemp = R' \ A(s,l)
All = A(l,l) - Atemp' * Atemp
基本上我正在计算部分高斯消元的尾随矩阵。 上面的代码,对于大而稀疏的A,几乎所有的时间花在第二行,特别是第三行。有关如何改善这一点的任何提示? s的大小通常是l的大小的10%。
非常感谢