问题描述
我想写一个返回FORTRAN的可分配数组的函数
I want to write a function that returns an allocatable array in fortran
program test
implicit none
real a(3)
real, allocatable :: F18(:)
a = (/1,2,3/)
print *, F18(a)
end program test
function F18(A)
implicit none
real A(:) ! An assumed shape array
real F18(size(A,1)) ! The function result itself is
! the second dimension of A.
F18 =A !
end function F18
预计打印1 2 3屏幕上,但我得到了一个错误:
It is expected to print "1 2 3" on the screen but I got an error:
forrtl:严重(157):程序异常 - 访问冲突
这是什么问题?
另外,我已经试过code是这样的:
Also, I have tried code like this:
program test
implicit none
real a(3)
real, allocatable :: F18(:)
a = (/1,2,3/)
print *, F18(a,3)
end program test
function F18(A,n)
implicit none
integer n
real A(:) ! An assumed shape array
real F18(size(A,1)) ! The function result itself is
! the second dimension of A.
F18 =A !
end function F18
在编译期间我得到了:
During compiling I got:
Intel(R) Visual Fortran Intel(R) 64 Compiler XE for applications running on Intel(R) 64, Version 14.0.4.237 Build 20140805
Copyright (C) 1985-2014 Intel Corporation. All rights reserved.
D:\Fortran\Elephant.f90(6): error #6351: The number of subscripts is incorrect. [F18]
print *, F18(a,3)
-------------^
compilation aborted for D:\Fortran\Elephant.f90 (code 1)
我真的很困惑与FORTRAN的功能。
I am really confused with fortran's function.
什么是调用FORTRAN数组值函数的正确方法?
What is the right way to call Array-valued Functions in fortran?
@Fortranner
@Fortranner
推荐答案
您需要让别人知道给调用者函数的性质。最简单的方法就是把它变成一个模块,并使用该模块。在你的例子,在你的主程序你声明数组'F18',这是没有的功能。
You need to make the properties of the function known to the caller. The easiest way is to put it into a module and 'use' that module. In your examples, in your main program you are declaring an array 'F18', which is not the function.
module mystuff
contains
function F18(A,n)
implicit none
integer n
real A(:) ! An assumed shape array
real F18(size(A,1)) ! The function result itself is
! the second dimension of A.
F18 =A !
end function F18
end module mystuff
program test
use mystuff
implicit none
real a(3)
a = (/1,2,3/)
print *, F18(a,3)
end program test
这篇关于如何调用数组值函数在FORTRAN?的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持!