Mathematica的“精确解”和“数值解”得到的本征矢为何不一样？(4)

Mathematica的“精确解”和“数值解”得到的本征矢为何不一样

再往下进行，可以发现wk[GPi0]和wk[GPiPi]不是反对称矩阵了，所以它们的Pfaffian无法计算，如下

这样，我们看到“数值解”出现了问题！为了进一步查找原因，我们选择“精确解”的GPiPi点进行测试。

测试

“精确解”的GPiPi点为

它的“精确解”和“数值解”对应的矩阵分别定义为t1和t2，可以看到在Chop函数下，它们是相等的（用==判断），但是在精度上它们是不相等的（用===判断）。