Thomas算法是否适用于含复数元素的非厄米非正定三对角矩阵？

Thomas算法（追赶法）本质是三对角矩阵LU分解的高效实现，完全适用于元素为复数的三对角矩阵，其核心要求仅为：

• 矩阵是**非奇异（可逆）**的
• LU分解过程中不会出现主元为0的情况（无选主元版本的限制）

针对你提到的场景补充说明：

• 矩阵是否为厄米矩阵、是否正定，不影响Thomas算法的适用性——只要是结构正确的可逆三对角矩阵，不管元素是实数还是复数，算法逻辑都成立。
• 带完美匹配层（PML）的微分方程时间步求解中，矩阵出现复数元素是正常情况，Thomas算法本身可以处理这类场景。

如果你实现后结果不符合预期，建议从以下方向排查：

• 检查复数运算逻辑：确保复数的加减乘除、赋值等操作没有错误（比如除法的分子分母顺序，复数乘法的实部虚部计算）
• 验证矩阵结构：确认输入的矩阵严格符合三对角结构，没有非三对角位置的元素被误纳入计算
• 做基准测试：构造一个简单的可逆复数三对角矩阵，用你的实现求解，再与scipy.linalg.solve_banded等成熟工具的结果对比，定位差异点
• 排查数值稳定性：如果矩阵主元过小，无选主元的Thomas算法可能会出现数值误差，可检查分解过程中的主元大小，必要时改用带选主元的三对角求解方法

内容的提问来源于stack exchange，提问作者Tom