用 FLAME 方法系统推导斜对称矩阵 LTLT 分解算法,融合 BLAS 类运算后性能较 PFAPACK 与 Pfaffine 大幅提升
该工作用 FLAME 形式化方法系统推导斜对称矩阵 X = LTLT(L 为单位下三角、T 为斜对称三对角)三角三对角化的一族算法,给出无主元与带主元的多种分块右看、左看及融合变体,识别出新的 level-2 与 level-3 BLAS 类运算,并借助 BLIS 2.0 的打包机制与 OpenMP 并行实现,实验表明其最佳实现性能大幅超过此前唯一已知软件 PFAPACK 与 Pfaffine,同时与相关对称矩阵分解软件相当或更优。