随机数组 & 接口冲突

hanshikai

GFortran编译通过，但每次运行结果一样，无变化
IVF无法通过编译，并提示【there is a conflict beetween local interface and external interface [GET_RANDOM]】

文本为输入文件及运行结果[random.dat]
RANDOM.DAT (113 Bytes, 下载次数: 1)

2014-4-4 17:30 上传

点击文件名下载附件

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

PROGRAM RANDOMNUM
IMPLICIT NONE
INTERFACE !为函数GET_RANDOM(N,L,U)定义接口
!==================================================
FUNCTION GET_RANDOM(N,L,U)
IMPLICIT NONE
INTEGER N,I
REAL L,U,LENGTH,t
REAL::GET_RANDOM(N)!除了N,L,U外应该不许要声明其他吧？？
!RETURN
END FUNCTION 
!==================================================
END INTERFACE
INTEGER N
REAL L,U
REAL,ALLOCATABLE::A(:)!因为我想让这个函数具有通用性，所以用了动态数组

CALL RANDOM_SEED()!随机数种子
!石头大人说可能是没有种子的原因，可是我这里一直都有哇。
!我用的GFortran，后来在IVF上试验，无法编译
OPEN(9,file='random.dat')
READ(9,*)N,L,U
ALLOCATE(A(N))
A=GET_RANDOM(N,L,U)
WRITE(9,90)A
90 FORMAT(/,F7.3)
DEALLOCATE (A)
END

!==================================================
FUNCTION GET_RANDOM(N,L,U)
IMPLICIT NONE
INTEGER N,I
REAL L,U,LENGTH,t
REAL::GET_RANDOM(N)
!L,U分别代表了数组的上下限
LENGTH=U-L
DO I=1,N
CALL RANDOM_NUMBER(t)
GET_RANDOM(I)=L+LENGTH*t
END DO

RETURN
END FUNCTION 
!==================================================

hanshikai

不好意思，刚刚给代码加上了注释。已更新。

hanshikai

chuxf 发表于 2014-4-4 18:06
1. GFortran 的 seed 函数需要费点劲，你需要手动获得系统时间并当做种子。它的帮助里给了一个例子，可以复 ...

非常感谢

PS:我的ivf是11.*的

hanshikai

chuxf 发表于 2014-4-4 18:06
1. GFortran 的 seed 函数需要费点劲，你需要手动获得系统时间并当做种子。它的帮助里给了一个例子，可以复 ...

谢谢你指导性的意见【返回数组的话，通常用subroutine，避免使用function。如果非要用 function，可以封装到 module 里】

还有个问题：如果不用interface、module，是不是function不能返回数组？

hanshikai

chuxf 发表于 2014-4-4 18:40
是的，不能。

你把一切都否定了，自然就不能了。

谢谢，明白了。

hanshikai

pasuka 发表于 2014-4-4 18:38
可以预先定义结构体，里面设置一个动态分配数组
然后function里面分配动态数组并返回该结构体，只是有必 ...

谢谢。
结构体是什么？ 英文？

hanshikai

chuxf 发表于 2014-4-4 18:51
C 里面叫 struct
Fortran 里面叫 type（中文通常叫派生类型）

3Q.又长知识了。