24小时热门版块排行榜    

查看: 594  |  回复: 3

superrice

木虫 (正式写手)

[求助] fortran菜鸟求帮助 已有1人参与

我想练习一下MKL库里的lapack95解本征方程,看了下manual半天没看懂,感觉代码写的不太对,麻烦专家帮忙看下哪里不对,谢谢。

program main
use lapack95
use blas95
implicit none

real::A(4,4)=(/ 1,2,3,4, &
                 2,3,4,5, &
                 3,4,5,6, &
                 4,5,6,7 /)
real::B(3),E(4)
integer::i

call sytrd(A,B,'U')
call orgtr(A,B,'U')
call rsteqr(E,B,A,'I')

do i=1,4
print '(4F8.4)', E(i)
end do

do i=1,4
print '(4F8.4)', A(i,
end do

end
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嗨王and娜娜

木虫 (正式写手)

博士

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
superrice: 金币+5 2014-03-28 21:37:38
定义要放着在最开始,你都给数组赋值了,才定义,这能行嘛?
努力学习!
2楼2014-03-28 08:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superrice

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 嗨王and娜娜 at 2014-03-28 08:00:50
定义要放着在最开始,你都给数组赋值了,才定义,这能行嘛?

谢谢专家回复,我想了解mkl的rsteqr怎么用~
3楼2014-03-28 21:38:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嗨王and娜娜

木虫 (正式写手)

博士

引用回帖:
3楼: Originally posted by superrice at 2014-03-28 21:38:43
谢谢专家回复,我想了解mkl的rsteqr怎么用~...

你直接调用就行了呗。http://www.netlib.org/lapack/explore-3.1.1-html/dsteqr.f.html,这个上面 有例子,你应该能看懂的。
努力学习!
4楼2014-03-28 22:16:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 superrice 的主题更新
信息提示
请填处理意见