24小时热门版块排行榜    

查看: 396  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ilovexiaomu

金虫 (小有名气)

[求助] 麻烦看哈这个很简单的求线性方程的问题? 已有1人参与

各位,遇到一个程序,在网上直接下载的,编译没错,但是运行时报错,不知什么原因,还望各位解答
!**********************************************
! THIS EXAMPLE USES THE LAPACK ROUTINE DGESV
! TO SOLVE A SYSTEM OF LINEAR EQUATIONS AX=B
! A = [1, 2, 3; 4, 5, 6; 7, 8 10]            
! B = [1, 0; 0, 1; 0, 0]
! *********************************************
program lapack_test
  integer ipiv(3), info, i, j
  double precision A(3,3), B(3,2)
  A(1,1)=1
  A(1,2)=2
  A(1,3)=3
  A(2,1)=4
  A(2,2)=5
  A(2,3)=6
  A(3,1)=7
  A(3,2)=8
  A(3,3)=10
  
  B(1,1)=1
  B(2,1)=0
  B(3,1)=0
  B(1,2)=0
  B(2,2)=1
  B(3,2)=0

  call dgesv (3,2,A,3,ipiv,B,3,info)
  if(info .EQ. 0) then
     do i=1,3
        write(*,'(2F8.3)') (B(i,j), j=1,2)
     enddo
  endif
end program lapack_test
不知道是自己的安装软件有问题还是那里出了什么问题,还望各位解答
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖


感谢参与,应助指数 +1
ilovexiaomu: 金币+1, 有帮助, 感谢 2015-05-21 16:05:55
原因很简单,Lapack没装成功。。。。。。。。。。这东西装起来还是有点麻烦的,好好整几天吧

如果不是特别密集的矩阵计算而且要求计算特别高效,可以不装这东西,另外找其他的子函数调用即可
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
3楼2015-05-21 16:01:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

ilovexiaomu

金虫 (小有名气)

问题在于call dgesv()这个地方,按照道理他是一个函数库,怎么没有其具体形式呢?
2楼2015-05-21 15:40:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见