Matlab 线性系统中矩阵的部分消去

我想知道在matlab中消除（在线性系统求解的意义上）矩阵的某些行最有效的方法是什么。基本上，给定行“s”和“l”（例如s和l是不同的，它们的并集是1:n），如何有效地计算

A是大的nxn稀疏SPD矩阵

R = chol(A(s,s))
Atemp = R' \ A(s,l)
All = A(l,l) - Atemp' * Atemp

基本上我是在计算部分高斯消去法的尾随矩阵。 上面的代码，对于大而稀疏的A，几乎所有的时间都花在第二行，尤其是第三行。关于如何改进这一点有什么提示吗？ s的尺寸通常为l尺寸的10%

非常感谢