作者sam92791432 (sammy)
看板Fortran
标题Compaq fortran可执行，但Intel Fortran却不行
时间Thu Jun 5 08:42:03 2014
版上各位大大好：
小弟刚从Compaq fortran 转用 Intel Fortran
不过原来在Compaq可以使用的执行的Code
在Intel Fortran却不行
Code如下：
************************************************************************
* *
* Inverse matrix Calculating module (Including some subroutines) *
* *
************************************************************************
module LinearAlgebra
implicit none
contains
************************************************************************
* *
* Main Inverse Matrix subroutines *
* *
************************************************************************
subroutine main_inverse (A,IA)
implicit none
real*8 :: A(:,:),IA(:,:)
real*8,allocatable ::B(:,:)
integer ::i,j,N
N = size(A,1)
allocate(B(N,N))
!先把IA设定成单位矩阵
forall(i=1:N,j=1:N,i==j) IA(i,j)=1.0
forall(i=1:N,j=1:N,i/=j) IA(i,j)=0.0
!保存原本的矩阵A,使用B来计算
B=A
!把B化成对角线矩阵(除了对角线外,都为0)
call Upper(B,IA,N) !先把B化成上三角矩阵
call Lower(B,IA,N) !再把B化成下三角矩阵
!求解
forall(i=1:N) IA(i,:)=IA(i,:)/B(i,i)
return
end subroutine
************************************************************************
* *
* UPPER Matrix subroutines *
* *
************************************************************************
subroutine Upper(M,S,N)
implicit none
integer :: N
real*8 :: M(N,N)
real*8 :: S(N,N)
integer :: I,J
real*8 :: E
do I=1,N-1
do J=I+1,N
E=M(J,I)/M(I,I)
M(J,I:N)=M(J,I:N)-M(I,I:N)*E
S(J,:)=S(J,:)-S(I,:)*E
end do
end do
return
end subroutine Upper
************************************************************************
* *
* LOWER Matrix subroutines *
* *
************************************************************************
subroutine Lower(M,S,N)
implicit none
integer :: N
real*8 :: M(N,N)
real*8 :: S(N,N)
integer :: I,J
real*8 :: E
do I=N,2,-1
do J=I-1,1,-1
E=M(J,I)/M(I,I)
M(J,1:N)=M(J,1:N)-M(I,1:N)*E
S(J,:)=S(J,:)-S(I,:)*E
end do
end do
return
end subroutine Lower
end module
************************************************************************
* *
* The connecting-subroutine *
* *
************************************************************************
subroutine inverse(B,BI,M)
use LinearAlgebra
implicit real*8 (a-h,o-z)
real*8 :: A(M,M)
real*8 :: IA(M,M)
integer :: i
real*8 :: B(M,M),BI(M,M)
do imn = 1,M
do jmn = 1,M
A(imn,jmn) = B(imn,jmn)
end do
end do
call main_inverse(A,IA)
do imn = 1,M
do jmn = 1,M
BI(imn,jmn) = IA(imn,jmn)
end do
end do
return
end
执行时，会出现下列状况，
allocatable array is already allocated
在这应该是指我的B矩阵。
不晓得是不是语法上出了点问题，还请各位大大指教。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.115.62.41
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/Fortran/M.1401928926.A.71E.html
1F：推 aix0209:module B最後要解掉? 06/05 10:09
2F：→ darrenmm:你是不是开错档了? 06/05 12:34