将Fortran POINTER作为模块过程传递给SUBROUTINES

如何解决将Fortran POINTER作为模块过程传递给SUBROUTINES

我的第一个问题是：“可以将POINTER传递到INTENT旁边的子例程吗？”

尽管用INTENT（IN）或INTENT（INOUT）传递Fortran POINTER似乎是一种常见的做法，但是Intel for Fortran POINTER的文档说（https://software.intel.com/content/www/us/en/develop/documentation/fortran-compiler-developer-guide-and-reference/top/language-reference/a-to-z-reference/o-to-p/pointer-fortran.html）

不能使用INTENT或PARAMETER属性声明指针

这让我很困惑。

我的第二个问题与以下代码有关。

我尝试将write_matrix重载到write_matrix_2d（用于2d数组）和write_matrix_2d_ptr（重载用于2d Fortran POINTER）。

原来的代码无法编译

Ambiguous generic interface WRITE_MATRIX: previously declared specific procedure WRITE_MATRIX_2D is not distinguishable from this declaration. [WRITE_MATRIX_2D_PTR]

如果我评论(1)，则它可以执行并运行良好（但是我认为这是一种不好的做法，因为POINTER参数进入了需要数组输入的subroutien write_matrix_2d中）。

有没有办法像C / C ++一样通过其参数的数据类型来重载Fortran子例程？

module mod_write

    implicit none

    interface write_matrix
        module procedure write_matrix_2d
(1)     !module procedure write_matrix_2d_ptr
    end interface

contains

    subroutine write_matrix_2d(a)

        implicit none
        real,intent(in) :: a(:,:)
  
        ! local variables
        integer :: i,j

        print *,"in ARRAY printing subroutine"
        do i = lbound(a,1),ubound(a,1)
            write(*,*) ( a(i,j),j = lbound(a,2),2) )
        enddo
  
    end subroutine

    subroutine write_matrix_2d_ptr(a)

        implicit none
        real,pointer,"in POINTER printing subroutine"
        do i = lbound(a,2) )
        enddo
  
    end subroutine

end module

program test

    use mod_write

    implicit none

    real,dimension(9)            :: b = (/21,22,23,24,25,26,27,28,29/)
    real,dimension(3,3)          :: c
    real,3),target  :: ct(3,3)
    real,dimension(:,:),pointer :: cptr(:,:)

    c = reshape( b,(/3,3/) )

    write(*,*)
    write(*,*) 'c'
    call write_matrix(c)
(2) !call write_matrix_2d(c)

    ct = c
    cptr => ct(2:3,2:3)

    write(*,*) 'ct'
    call write_matrix(ct)
(2,3)!call write_matrix_2d_ptr(ct)

    write(*,*) 'cptr'
    call write_matrix(cptr)
(2) !call write_matrix_2d_ptr(cptr)

end program

我的第三个问题再次与所附代码有关。重载并不是什么大问题，因为我可以仅将write_matrix_2d和write_matrix_2d_ptr用作代码中的(2)。但是，正如您在(3)上看到的那样，call write_matrix_2d_ptr(ct)在target输入ct的情况下运行良好，而write_matrix_2d_ptr期望pointer的输入。这样可以吗，还是我应该再做一个write_matrix_2d_target？

谢谢。

解决方法

编译器文档错误（或严重过期）。作为指针的伪参数也可以具有intent属性。该属性适用于参数的指针关联状态，而不适用于目标。

如果通用接口中的过程的伪参数是可区分的，则它们可以被重载。如果伪参数具有不同的类型，种类或等级，则可以区分。某物是否为指针不是其类型的一部分。（为了完整起见，如果一个是没有意图的指针，而另一个是可分配的，则参数也是可区分的。）对通用接口中不同过程的可区分伪参数的要求意味着，特定泛型围绕特定过程的规则引用解析非常简单。

您的重载示例将是模棱两可的-可以将（n个关联的）指针实际参数传递给非指针虚拟对象-实际参数指定指针的目标。还可以将具有target属性的非指针实际值传递给intent（in）指针虚拟对象。

某物是否为指针与它是否为数组完全正交。

Fortran中的指针主要用于在程序执行期间引用不同的事物（目标）。

仅当您需要在与指针的关联状态相关的过程中执行某些操作时，虚拟参数才需要是指针-也许您想测试关联状态或将指针与其他东西关联。

（1）意图限制是某些复制和粘贴错误，或者是在某些我不理解的奇怪上下文中。 Fortran 90中有这样的限制。

（2）“但我想这是不好的做法，因为POINTER参数进入需要数组输入的subroutien write_matrix_2d中。”

将Fortran指针变量作为普通变量传递并没有什么不好。默认的参数传递机制是使用普通的非ointer机制传递thee指针的目标。只有在有特定原因的情况下，才使用伪虚拟参数。大多数情况下，tbat意味着在子例程中使用一种关联状态进行操作。还有其他用途，但是您需要特定的理由将参数声明为指针，否则只会使事情复杂化。

（3）target可以在Fortran 2008中传递给pointer,intent(in)。

将Fortran POINTER作为模块过程传递给SUBROUTINES

如何解决将Fortran POINTER作为模块过程传递给SUBROUTINES

解决方法

相关推荐