24小时热门版块排行榜    

查看: 2705  |  回复: 8

quantumage

金虫 (小有名气)

[求助] 怎么找出lapack函数及函数的使用方法

我的电脑上已经安装好了lapack3.2.1版程序包,用的是linuix系统gfortran编译器!由于我刚开始接触,不知道怎么在lapack中找出我需要的函数及函数的使用方法,有那位兄台知道给指点一下啊!我主要是想求一个矩阵的本征值,是一般的矩阵,但是矩阵的本征值里含有复数!我查到一个是dgeev,但不知道是怎么用,格式是什么样!要是知道的能否讲的详细一点啊!谢谢了!
回复此楼

» 收录本帖的淘帖专辑推荐

matlab和maple专业交流

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

study hard every day!
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chembetsey

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
参看LAPACK的手册,官网上有详细的说明和实例。
2楼2012-12-06 07:57:23
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

quantumage

金虫 (小有名气)

引用回帖:
2楼: Originally posted by chembetsey at 2012-12-06 07:57:23
参看LAPACK的手册,官网上有详细的说明和实例。

我已经在官网上查过了,还是没有搞清楚!知道的给点具体的意见吧!
study hard every day!
3楼2012-12-06 11:16:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chembetsey

木虫 (小有名气)

你去看看Intel的MKL手册,上面对每个函数的使用参数都有解释。只要按要求写好矩阵,正确的调用就可以了。难道要别人帮你写出来?
4楼2012-12-06 11:23:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

quantumage

金虫 (小有名气)

我找到了dgeev函数的说明,但是有几个参数看不明白!
我是这样理解的
fortran77  call dgeev(jobvl, jobvr,n,a,lda,wr,wi,vl,ldvl,vr,lvdr,work,lwork,info)
jobvl 和jobvr 分别表示是否计算左、右本征矢量
n是矩阵a的阶
a是要计算的矩阵
lda是a的第一维数(这里不是很清楚,如果矩阵a(2,3)那么这里的lda就是2,是这样的吗)?
wr、wi分别是本征值的实部和虚部
vl、vr分别是左右本征矢量,ldvl、ldvr分别表示左右本征矢量的维数
work是计算用的空间数组,lwork是空间数组的维数
info是用来表示计算结果参数,0,i,-i可以用来表示计算成功与否(这里也不是很清楚,我看后面的解释是一般可以选-1)
不知道我这样的理解对不对,我试着按这个方法算了一下(输出的时候wr(j),我就想看一下本征值),但是结果输出的时候总是段错误,不知道为什么?
另外让我迷惑的是lapack的解释分别用了输入参数和输出参数,但我怎么知道这个函数括号里的那个是输入参数,那个是输出参数呢?
study hard every day!
5楼2012-12-07 17:29:03
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

quantumage

金虫 (小有名气)

引用回帖:
4楼: Originally posted by chembetsey at 2012-12-06 11:23:02
你去看看Intel的MKL手册,上面对每个函数的使用参数都有解释。只要按要求写好矩阵,正确的调用就可以了。难道要别人帮你写出来?

谢谢你,按照你的线索我找到了lapack的使用手册!
study hard every day!
6楼2012-12-07 17:31:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

quantumage

金虫 (小有名气)

program main
integer i,j,k
real a(3,3),y(3,3),wr(3),wi(3),vl(3,3),vr(3,3)
do i=1,3
  do j=1,3
  a(i,j)=i+j
  end do
end do
do k=1,3
  write(*,*) a(k,
end do
call dgeev('N','V',3,a(3,3),3,wr(3),wi(3),vl(3,3),3,vr(3,3),3,y,-1,x)
write(*,*) wr
end program main
这一次调试成功了,但是为什么算出来的本征值不对啊!我用matlab对比了一下,怎么不一样啊,原因出在那里啊?
study hard every day!
7楼2012-12-07 18:40:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chembetsey

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jjdg: 金币+1, 感谢提供建议 2012-12-08 01:08:06
quantumage: 金币+25, ★★★★★最佳答案, 按照你的指示,最终我算出来了!太感谢了!好心人,谢谢你! 2012-12-10 11:25:28
引用回帖:
7楼: Originally posted by quantumage at 2012-12-07 18:40:36
program main
integer i,j,k
real a(3,3),y(3,3),wr(3),wi(3),vl(3,3),vr(3,3)
do i=1,3
  do j=1,3
  a(i,j)=i+j
  end do
end do
do k=1,3
  write(*,*) a(k,
end do
call dgeev('N','V',3,a(3,3),3 ...

你的调用方式不对。你需要先调用一次,得到work数组的大小,然后再将此大小的数组传进函数里。这是MKL中LAPACK函数的使用方式。
8楼2012-12-07 22:46:45
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

quantumage

金虫 (小有名气)

引用回帖:
8楼: Originally posted by chembetsey at 2012-12-07 22:46:45
你的调用方式不对。你需要先调用一次,得到work数组的大小,然后再将此大小的数组传进函数里。这是MKL中LAPACK函数的使用方式。...

太完美了,不但实根解出来了,连虚根都一起解出来了!真的好感谢你!
study hard every day!
9楼2012-12-10 12:21:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 quantumage 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +8 ziyangfang 2026-08-25 11/550 2026-09-01 09:51 by zzuzxg
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见