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