06 2014 档案

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