作者bmw530li ( )
看板Fortran
标题[问题] function如何回传阵列(矩阵)
时间Fri Jun 1 00:14:48 2012
想使用function回传一整个阵列而不是一个数值，不知到与法该怎麽用
program main
real*8, allocatable :: a(:,:),a1(:,:))
integer*4 :: b=3
allocate (a(b,b),a1(b,b))
!allocate (a1(b,b))
allocate (a2(b,b))
a=0.d0; a1=0.d0;
a1(1,1:3)=1.5d0;
a(2,1:3)=2d0*a1(1,1:3)+5d0
write(*,*) a;
write(*,*) a1;
print * , '\'
write(*,*) add(b,a,a1)
pause
stop
end program main
function add(b,a,a1)
integer*4 :: b
real*8,allocatable :: a(:,:),a1(:,:)
allocate( a(b,b),a1(b,b) )
real*8, dimension(b,b) :: add
add=a+a1
return
end function add
跑出两个错误讯息，不知道与法该怎麽用QQ
Error: A specification statement cannot appear in the executable section.
real*8, dimension(b,b) :: add
--------^
Error: The shapes of the array expressions do not conform.
add=a+a1
--------^
除此之外，顺便问一下可以自动加入行号吗....难道真的要一行一行打QQ?
我用的是Compaq Visual Fortran
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.121.23.113
1F：→ biggun:第一个错误 请把宣告放在 allocate 前面 06/01 11:03
2F：→ biggun:第二个错误 属於基本问题 请看有关阵列运算的章节 06/01 11:04
3F：推 terryys:第二个错误其实就是第一个引起的，解决第一个就可以了 06/01 11:20
4F：→ terryys:但是副程式要用allocatable的话要加interface，或者把它 06/01 11:21
5F：→ terryys:用contains包在主程式内 06/01 11:21
6F：→ terryys:对了，你这个程式里还有另外的错误，但都容易解决 06/01 11:22
7F：→ bmw530li:所以要用到interface就是了？我再看看，感谢 06/01 15:10
8F：推 kermomo:我的经验是把这种函数写在MODULE再USE该MODULE就行了 06/12 03:52