24小时热门版块排行榜    

查看: 2690  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 02:26 by xxniao123
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见