可分配的Intent(inout)参数是否可选?

时间:2017-08-15 02:04:02

标签: fortran optional allocatable-array

我在尝试定义子例程时遇到问题,子例程的参数包含一个可分配的,可选的,intent(inout)变量,如下所示。代码编译正常,但得到运行时错误"分段错误 - 无效的内存引用"。

子程序test_routine.f90

SUBROUTINE test_routine(A,B)

IMPLICIT NONE

    REAL,ALLOCATABLE,INTENT(IN)               :: A(:,:)
    REAL,ALLOCATABLE,OPTIONAL,INTENT(INOUT)   :: B(:,:)

    B = A
    B(:,:) = 1

END SUBROUTINE

此子程序打包在一个模块中,并在Main中调用。

模块test_module.f90

MODULE test_module

IMPLICIT NONE

    INTERFACE test_routine
        MODULE PROCEDURE test_routine
    END INTERFACE

END MODULE test_module

主要test_main.f90

PROGRAM main

    USE test_module

IMPLICIT NONE

    REAL,ALLOCATABLE    :: A(:,:),B(:,:)

    ALLOCATE(A(6,6))
    ALLOCATE(B(6,6))
    A(:,:) = 0
    CALL test_routine(A,B)    ! WORKS FINE
    CALL test_routine(A)      ! DOESN'T WORK!

END PROGRAM main

然后我尝试分配另一个变量op_B来弥补B,如果主程序没有传入它,它实际上并不存在。但是下面的代码仍然无法工作

SUBROUTINE test_routine(A,B)

IMPLICIT NONE

    REAL,ALLOCATABLE,INTENT(IN)    :: A(:,:)
    REAL,ALLOCATABLE,OPTIONAL,INTENT(INOUT)   :: B(:,:)
    REAL,ALLOCATABLE               :: op_B(:,:)

    IF(.NOT. PRESENT(B)) THEN
        ALLOCATE(op_B(SIZE(A,1),SIZE(A,2)))
        B = op_B
    END IF
    B = A
    B(:,:) = 1

END SUBROUTINE

顺便说一句,我也尝试使用固定大小的阵列,但仍然没有帮助。我想知道做这些事情是不可能的吗?

2 个答案:

答案 0 :(得分:2)

虚拟参数可以是可选的,可分配的和意图的(inout)。

但是不允许定义或引用可选的伪参数,bar将其传递给PRESENT或将其与另一个可选参数相关联。如果“B”不存在,则无法执行“B = A”。

答案 1 :(得分:0)

基于@lanH的回答,我目前唯一的解决方法是为每个可选的伪参数赋予别名。这个版本现在有效。

DismissViewController