R求解多个线性系统,其中每个RHS是矩阵的一列

R solving multiple linear systems where each RHS is a column of a matrix

我知道我可以使用 R 中的 backsolve 求解系统 Lx = b 其中 L 是上三角。但是,现在我有一个矩阵 B 其中每个B 的列是 Lx = b 的右侧。基本上我想解决 Lx = b,其中 b 是矩阵 B 的每一列。有没有比循环更有效的方法?

我没有发布 MWE,因为我很确定应该有一个函数可以做到这一点。

或许你可以试试下面的方法

x <- ginv(L)%*%b