应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 统计、非线性&软物质 » 【求助】FORTRAN 中有没有计算稀疏矩阵本征值的函数？
北京石油化工学院2026年研究生招生接收调剂公告
121/2返回列表12下一页
查看: 1009  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

hedge

银虫 (小有名气)

[交流] 【求助】FORTRAN 中有没有计算稀疏矩阵本征值的函数？

如题，IMSL库中只有用稀疏矩阵解线性方程的函数！有没有解本征值的啊？

[ Last edited by hedge on 2009-3-9 at 16:28 ]
回复此楼

» 猜你喜欢
1楼 2009-03-06 19:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happymmin

银虫 (正式写手)
★ ★ ★ ★
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
赞一下(10人) 回复此楼
The wisdom of God gives us balance, insight and stability.
2楼2009-03-07 19:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happymmin

银虫 (正式写手)
好久没写程序了，我自己都忘得差不多了。
赞一下 回复此楼
The wisdom of God gives us balance, insight and stability.
3楼2009-03-07 19:01:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allawwp

金虫 (正式写手)
★
zt970831(金币+1,VIP+0):感谢您的交流 3-8 21:53
用intel 的mkl 数学库，或者 arpack
赞一下(10人) 回复此楼
5楼2009-03-08 13:07:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happymmin

银虫 (正式写手)
引用回帖:
Originally posted by argo at 2009-3-7 19:36:
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 argume ...

不是啊，调用imsl的那个求本征值的函数时，说明用那个作格式输出。你可以打开imsl文件看看。
赞一下 回复此楼
The wisdom of God gives us balance, insight and stability.
6楼2009-03-09 09:23:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hedge

银虫 (小有名气)
★ ★ ★
zt970831(金币+3,VIP+0):感谢您的交流 7-13 23:30
CALL EVLRG (N, A, LDA, EVAL)
Compute all of the eigenvalues of a real matrix.
但这不是解稀疏矩阵的本征值！！
不合我题目的要求！！
但还是感谢您的回答！
引用回帖:
Originally posted by happymmin at 2009-3-7 19:00:
有的。我用过。找了一个当时联系的程序

program LearnToUse   !学习使用内部函数求随机矩阵的本征值
use imsl
implicit none
integer,parameter::n=4,lda=n
real::a(n,n)
complex::eval(n)
real::ds
in ...

[ Last edited by hedge on 2009-3-9 at 15:35 ]
赞一下(20人) 回复此楼
7楼2009-03-09 15:31:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hedaors

禁虫 (正式写手)

隐清堂主
用 mkl 的数学库很简单的，你看一看 mkl  的 useguide 和 manual 就应该知道了，一般不推荐用 imsl , 没有免费版的，效率也不高
赞一下 回复此楼
8楼2009-03-10 08:10:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hedge

银虫 (小有名气)
但mkl没有免费版本。
不知有没有试用版呢？
赞一下 回复此楼
9楼2009-03-10 13:09:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hedaors

禁虫 (正式写手)

隐清堂主
引用回帖:
Originally posted by hedge at 2009-3-10 13:09:
但mkl没有免费版本。
不知有没有试用版呢？

mkl 有试用版，在 linux 下有非商业使用免费版。
赞一下 回复此楼
10楼2009-03-10 15:37:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hedge 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320调剂 +4 农业工程与信息� 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 一志愿南昌大学324求调剂 +13 hanamiko 2026-04-01 13/650 2026-04-03 18:30 by ls刘帅
[考研] 土木水利328分求调剂 +6 疾风知劲草666 2026-04-02 6/300 2026-04-03 11:38 by znian
[考研] 一志愿深大085601材料工程专业（专硕）300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 282求调剂 +3 aaa车辆 2026-04-02 3/150 2026-04-02 21:55 by zllcz
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[考研] 301求调剂 +4 细胞相关蛋白 2026-04-02 8/400 2026-04-02 16:30 by 271179835
[考研] 初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666！
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 材料调剂 +14 一样YWY 2026-04-01 14/700 2026-04-01 21:07 by lijunpoly
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 求调剂 +4 DADA怪 2026-03-31 4/200 2026-04-01 14:30 by ZXlzxl0425
[考研] 318一志愿吉林大学生物与医药 求调剂 +6 笃行致远. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 调剂考研 +3 王杰一 2026-03-29 3/150 2026-03-29 08:09 by fmesaito
信息提示
关闭
请填处理意见
关闭