| 查看: 937 | 回复: 11 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】FORTRAN 中有没有计算稀疏矩阵本征值的函数？
|
|||
|
如题，IMSL库中只有用稀疏矩阵解线性方程的函数！有没有解本征值的啊？ [ Last edited by hedge on 2009-3-9 at 16:28 ] |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有7人回复
-
请教限项目规定
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有15人回复
-
Materials Today Chemistry审稿周期
已经有6人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
推荐一本书
已经有13人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
9楼2009-03-10 13:09:31
★ ★ ★ ★
lvzhu2007(金币+4,VIP+0):谢谢交流 3-8 20:32
lvzhu2007(金币+4,VIP+0):谢谢交流 3-8 20:32
|
有的。我用过。找了一个当时联系的程序 program LearnToUse !学习使用内部函数求随机矩阵的本征值 use imsl implicit none integer,parameter::n=4,lda=n real::a(n,n) complex::eval(n) real::ds integer::iseed,i,j open(7,file="eigenvalue.dat"  data a/-2.0,-3.0,-2.0,-1.0,& 2.0,3.0,0.0,0.0,& 2.0,2.0,4.0,0.0,& 2.0,2.0,2.0,5.0/ do i=1,n write(*,*)(a(i,j),j=1,n) end do write(*,*)" " call evlrg(n,a,lda,eval) call wrcrn('eigenvalues',1,n,eval,1,0) write(*,*)" " do i=1,n write(*,10)eval(i) write(7,20)i,eval(i) end do 10 format(4x,"real",4x,f4.2,4x,"imig",4x,f4.2) write(*,*)" " 20 format(i4,8x,f4.2,8x,f4.2) close(7) stop end program |
2楼2009-03-07 19:00:41
3楼2009-03-07 19:01:31
argo
铁杆木虫 (著名写手)
上善若水
- 应助: 4 (幼儿园)
- 金币: 9672.8
- 散金: 200
- 红花: 34
- 帖子: 1940
- 在线: 486.4小时
- 虫号: 507533
- 注册: 2008-02-19
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
call wrcrn('eigenvalues',1,n,eval,1,0) wrcrn是自已编的吗？编译时出现以下提示： Error: If the actual argument is scalar, the corresponding dummy argument shall be scalar unless the actual argument is an element of an array that is not an assumed-shape or pointer array, or a subs tring of such an element. [TITLE] call wrcrn('eigenvalues',1,n,eval,1,0) |
4楼2009-03-07 19:36:49