SciPy solve_banded 矩阵的解法方法

Question

有没有人知道 scipy.linalg.solve_banded 这个方法是怎么解决方程组的？它的 文档 没有说明这个函数使用了什么解法。通常来说，处理这类方程组时会用到 托马斯算法，也叫做TDMA，但我想知道这个Scipy的函数是否用了其他的解法。

Answer 1

这个Github代码显示，scipy使用了lapack里的gbsv()函数来解决这个问题。你可以在这里和这里了解更多关于gbsv()的信息。

我不太确定这是不是和托马斯算法一样。看起来这两者都使用了LU分解的方法。