问题描述
我的第一个问题是:“可以将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)
这让我很困惑。
我的第二个问题与以下代码有关。
我尝试将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)
。