如何使用function重载运算符运算后返回一个数组

306908677

RT
求高手指点

306908677

补充：数组运算得到数组

百事可乐

我这里有一个重载 matmul 的代码,给你当例子吧

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	module matmul_int

02	INTERFACE OPERATOR(.mm.)

03	module procedure fc

04	END INTERFACE

05

contains

06

07	function fc(a,b)

08	real , intent( IN ) :: a(:,:) , b(:,:)

09	real :: fc(size(a,dim=1),size(b,dim=2))

10	fc = matmul(a,b)

11	end function fc

12	end module matmul_int

13

14	program main

15	use matmul_int

16	real :: x(3,2) , y(2,4)

17	x = 1 ; y = 2

18	write(*,*) x .mm. y

19

end

li913

普通的 + - * / 都是可以运用于数组的。

real a(3,4), b(3,4), c(3,4)
c = a + b

306908677

百事可乐 发表于 2014-12-6 16:13
我这里有一个重载 matmul 的代码,给你当例子吧

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	module matmul_int

02

[/quote]

03	[mw_shl_code=fortran,true]  program main

04	implicit none

05

06

interface

07	subroutine random_arr(a)

08	implicit none

09	integer(4),intent(out),allocatable ::a(:)

10

end

11	end interface

12

13	interface operator(.U.)

14	Function addy(c1, c2)

15	implicit none

16	integer, intent(in) :: c1(:), c2(:)

17	integer, intent(out):: addy(20)

错误：
Function addy(c1, c2)
1
错误： (1)处的符号不是一个 DUMMY 变量
/Users/StarkLee/Desktop/EQU2.F90:33.1:

function addy(c1,c2)
1
错误： (1)处的符号不是一个 DUMMY 变量
Jarvis:~ StarkLee$

什么原因啊＝＝

306908677

li913 发表于 2014-12-6 19:10
普通的 + - * / 都是可以运用于数组的。

real a(3,4), b(3,4), c(3,4)

嗯嗯，thanx～

fcode

照猫画虎还不会？

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01

Module S

02	Interface Operator (.U.)

03	Module Procedure Addy

04	End Interface Operator (.U.)

05

Contains

06	Function Addy(C1, C2)

07	Implicit None

08	Integer, Intent (In) :: C1(:), C2(:)

09	Integer :: Addy(20)

10	Integer :: Arry1(20), Arry2(20)

11	Integer I, Len

12	Logical :: Flg(20)

13	Len = Size(C1)

14	Arry1 = -1

15	Do I = 1, Len

16	Arry1(C1(I)) = C1(I)

17

End Do

18	Len = Size(C2)

19	Arry2 = -2

20	Do I = 1, Len

21	Arry2(C2(I)) = C2(I)

22

End Do

23

Addy = -3

24	Flg = All(Arry1.Eq.Arry2)

25	Print *, Flg

26	Addy = Merge(Arry1, Addy, Flg)

27	End Function Addy

28	End Module S

29	Program Main

30

Use S

31	Implicit None

32

33

Interface

34	Subroutine Random_arr(A)

35	Implicit None

36	Integer (Kind=4), Intent (Out), Allocatable :: A(:)

37	End Subroutine Random_arr

38	End Interface

39

40

41	Integer, Allocatable :: A(:)

42	Integer, Allocatable :: B(:)

43	Integer Lena, Lenb

44

45	Call Init_random_seed()

46	Call Random_arr(A)

47	Call Random_arr(B)

48	Lena = Size(A)

49	Lenb = Size(B)

50	Print *, 'A'

51	Write (*, 10) A

52	Print *, 'B'

53	Write (*, 10) B

54	10 Format (20I3)

55

56	End Program Main

306908677

fcode 发表于 2014-12-7 00:31
照猫画虎还不会？

[mw_shl_code=fortran,true]Module S

多谢大家的耐心指点～

仔细看了一遍原来function不能定义为 out～
问题已解决

vvt

是的，in,out只能对虚参进行。返回值不是虚参