| 查看: 478 | 回复: 3 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】arpack输出本征向量的问题
|
|||
|
请高手指点怎么在arpack里看本征向量？ 我用的是dsarpack，下面 call dsaupd ( ido, bmat, n, which, nev, tol, resid, & ncv, v, ldv, iparam, ipntr, workd, workl, & lworkl, info ) call HMATRIX(workd(ipntr(1)), workd(ipntr(2))) call dseupd ( rvec, 'All', select, d, v, ldv, sigma, & bmat, n, which, nev, tol, resid, ncv, v, ldv, & iparam, ipntr, workd, workl, lworkl, ierr ) call HMATRIX(v(1,j),ax) call daxpy(n, -d(j,1), v(1,j), 1, ax, 1) hmatrix是把哈密顿量作用在前一个变量上，输出后一个变量 根据我的肤浅理解本征向量应该在v(1,j)里面，但是又好像应该是ax，现在应经晕了  拜托拜托！ |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有12人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
happy_clp
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1618.9
- 帖子: 873
- 在线: 260.8小时
- 虫号: 435913
- 注册: 2007-08-25
- 专业: 金属材料的力学行为
2楼2009-04-09 21:02:49
3楼2009-04-10 19:35:05
4楼2009-05-19 22:13:57
20