| 查看: 2272 | 回复: 8 | ||||
quantumage金虫 (小有名气)
|
[求助]
怎么找出lapack函数及函数的使用方法
|
| 我的电脑上已经安装好了lapack3.2.1版程序包，用的是linuix系统gfortran编译器！由于我刚开始接触，不知道怎么在lapack中找出我需要的函数及函数的使用方法，有那位兄台知道给指点一下啊！我主要是想求一个矩阵的本征值，是一般的矩阵，但是矩阵的本征值里含有复数！我查到一个是dgeev，但不知道是怎么用，格式是什么样！要是知道的能否讲的详细一点啊！谢谢了！ |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab和maple专业交流 |
» 猜你喜欢
垃圾破二本职称评审标准
已经有18人回复
-
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
EST投稿状态问题
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab如何读取中括号内的数据
已经有8人回复
- MATLAB实用中文手册
已经有507人回复
- 特征函数的性质
已经有8人回复
- 矩阵的本征矢量
已经有6人回复
- MATLAB2010多元函数的泰勒展开怎么计算？
已经有4人回复
- word中如何输入矩阵的中括号？
已经有5人回复
- fortran中如何调用matlab中求解特征值函数？
已经有9人回复
- 高性能计算函数库:LAPACK;BLAS;GOTO;Atlas;FFTW 这些能安装到WIN7系统上吗？
已经有6人回复
- 卡方分布的特征函数怎么求啊？
已经有6人回复
- 如何寻找最逼近的曲线
已经有3人回复
- LDA的使用方法
已经有10人回复
- 如何分离函数的实部和虚部？
已经有4人回复
- 安装LAPACK出错，该怎么办
已经有3人回复
- 矩阵求特征值和特征向量用arpack和lapack哪个好些？
已经有11人回复
- 泰勒公式的问题！多项式怎么敢逼近一个函数 ，多项式与函数的关系
已经有7人回复
- 求助ifort 调用lapack的问题！
已经有4人回复
- 上传一本《固体物理中的格林函数方法》希望对大家有用
已经有168人回复
- 想调用fortran的BLAS函数库。。。忙活半天还是不行
已经有5人回复
- fortran程序怎样调用lapack库进行矩阵对角化？
已经有7人回复
- 为什么画出的对关联函数第三峰刚开始,之后的函数值就全变为零了,怎么解决这个问题
已经有5人回复
- 行列式的值，矩阵相乘
已经有6人回复
- 【分享】MKL程序编译与连接：Lapack篇
已经有6人回复
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
2楼2012-12-06 07:57:23
quantumage
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 900.8
- 帖子: 292
- 在线: 136.8小时
- 虫号: 1855704
- 注册: 2012-06-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2012-12-06 11:16:30
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
4楼2012-12-06 11:23:02
quantumage
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 900.8
- 帖子: 292
- 在线: 136.8小时
- 虫号: 1855704
- 注册: 2012-06-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
我找到了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的解释分别用了输入参数和输出参数，但我怎么知道这个函数括号里的那个是输入参数，那个是输出参数呢？ |
5楼2012-12-07 17:29:03
quantumage
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 900.8
- 帖子: 292
- 在线: 136.8小时
- 虫号: 1855704
- 注册: 2012-06-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
6楼2012-12-07 17:31:56
quantumage
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 900.8
- 帖子: 292
- 在线: 136.8小时
- 虫号: 1855704
- 注册: 2012-06-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
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对比了一下，怎么不一样啊，原因出在那里啊？ |
7楼2012-12-07 18:40:36
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
8楼2012-12-07 22:46:45
quantumage
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 900.8
- 帖子: 292
- 在线: 136.8小时
- 虫号: 1855704
- 注册: 2012-06-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
9楼2012-12-10 12:21:21