我正在使用 LAPACK zgeev 例程来获取 C++ 中对称矩阵的特征值和特征向量。问题是 zgeev 在循环中被调用,但它有时对特征值(和特征向量)进行不同的排序。
例如,这是第一轮循环的特征值:
(-1.29007e-5 - 5.207e-6*i) (1.28782e-5 + 7.40505e-6*i)
这是第二次的结果:
(1.28782e-5 + 7.40505e-6*i) (-1.29007e-5 - 5.207e-6*i)
我需要将这些特征值和向量的演变绘制为循环变量的函数,但它们每次都会被交换,并为我提供我需要的图的组合。
如何解决这个问题?