用于Scipy solve_banded矩阵的求解方法

Question

有谁知道scipy.linalg.solve_banded用什么方法来解决方程组？ documentation未说明函数使用的解法。通常Thomas algorithm，a.k.a。TDMA用于这些类型的系统，但我想知道这个Scipy函数是否使用了其他解决方法。

Answer 1

Github code表明scipy使用lapack例程gbsv()来解决此问题。您可以阅读gbsv() here和here。

我不确定这是否与托马斯算法相同。看起来两者都使用LU分解。