复杂对称矩阵的Autonne-Takagi分解 该代码最初是为工具箱编写的,但最终不需要。 因此,我根据LGPL许可证发布了它。 我还没有看到任何混合方法的实现,但是它似乎大大优于其他方法。 通过使用矢量化,可以进一步优化takagi_jacobi脚本。 但是,请注意,对Jacobi旋转的特定实现是为了使舍入误差最小。 因此我们没有试图改变它。 (C)丹尼斯·罗塞(Denis Rosset),2021年