MLDIVIDE提供了一个优秀的稀疏矩阵三对角矩阵求解器,以及Thomas算法的其他实现(详见下文),但我需要更快的方法来解决复杂数据的三对角系统。在我的系统(R2018b)上,这种方法比MLDIVIDE或MATLAB中的直接实现快四倍。它使用AVX指令处理复数交错操作,编译后只需将其添加到路径上即可使用:“mex -R2018a 'CFLAGS=-mavx' tdma.c”。