小对称矩阵的特征值

计算科学 数值分析 特征值 高性能计算
2021-12-14 06:06:28

如果我要解决对称特征值系统A=QDQT, 在哪里ARn×nn很小(在 4 - 64 范围内);我想要所有的特征向量和特征值;

我的设计中有两个主要考虑因素;我希望我的实施尽可能快;我需要Q尽可能正交,但是我可以容忍特征值的误差。什么是解决它的好方法?目标平台是单核 x86 处理器;

所有 MKL 例程,如 sYEV、sYEVD 和 sYEVR 似乎都非常慢(根据失败率,大约是预期的 10 倍)。因此需要自己实现它;我希望设计空间不仅仅是三对角归约类的方法;

我想对所有三个方面提出建议,即此类问题的算法、实现和现有实现;

1个回答

以我的经验,主要的 BLAS 库没有针对小型n,但您不应该期望获得 10 的全部因子。您可以查看Eigen以了解它在比较中的表现。

您将被包含在较大尺寸的 L2 缓存中,而 L1 将包含在较小尺寸的缓存中。您可能需要考虑到这一点才能获得良好的性能。