Fortran Ifort-mkl-lapack特征值

我正在使用ifort 14.0、mkl和f95 lapack接口进行特征值分解

我当前使用的函数是DSYEV_F95，它被称为：

        call DSYEV_F95(A=RMEigenVectors,W=RVEigenValues,Info=ISTmp&
      &,JobZ="N")

仅产生本征值，或

        call DSYEV_F95(A=RMEigenVectors,W=RVEigenValues,Info=ISTmp&
      &,JobZ="V")

它还返回特征向量

然而，使用call one，给定100x100矩阵的前三个特征值为：

-0.00000000000000022204 1.01690323906836455059 1.0405135339583818602

然而，如果使用调用2，结果会发生变化：

0.00000000000000246374 1.01690323906836477263 1.0405135339583885216

关于第一个值符号变化的任何想法这个开关完全是调用的函数，因此使用调用二，到目前为止还没有出现负特征值

---

谢谢

你检查信息的值吗？但是第一个特征值的值太小了，我认为它只是0，考虑到数值精度，因为它比其他的小16个数量级——64位双精度的范围。阅读有关浮点数的内容。问题不在于数字太小，而在于函数调用时值的差异。我可以用完全相同的数字和完全相同的差异一次又一次地重现结果