用于 Scipy solve_banded 矩阵的求解方法

Question

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

Answer 1

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

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