摘要： numericalrecipe里一共讲了两种实数对称矩阵的对角化，jacobi法tred2生成上三角阵以后用tqli对角化前者稳定但慢易并行，后者较快但疑似不稳定，串行。花了一下午，一点点调试终于知道了第二种方法不稳定的原因在哪里 1 SUBROUTINE tred2(a,d,e,n... 阅读全文