Python 如何在sympy中重用LU分解来求解常数M和变化b的Mx=b

求解辛矩阵Mx=b的最简单方法是。如果我有很多常数M和变化b的方程，有没有办法重用LU分解？我看不到。

使用LUdecomposition获得L和U矩阵，然后依次用Ly=b和Ux=y解Mx=b，速度够快吗？这就是我在高层次上计划做的，但是我没有看到任何内置函数可以用Ly=b和Ux=y对上下三角矩阵进行消元求解。我希望Symphy能够识别三角矩阵并更有效地求解它们，但是看起来U.LUsolveL.LUsolveb实际上比较慢。有/可能会有帮助。卡尔哈特：把它变成一个答案，我会接受的。不知怎的，我错过了那些。