对于主对角线,请小心,您可以执行以下操作:
PROGRAM diagonal
IMPLICIT NONE
REAL, TARGET :: array(4,4)
REAL, POINTER :: ptr(:)
INTEGER :: i
array = RESHAPE([(i,i=1,SIZE(array))], SHAPE(array))
CALL get_diagonal_pointer(array, SIZE(array, 1), ptr)
PRINT "(*(G0,:,','))", ptr
CONTAINS
SUBROUTINE get_diagonal_pointer(arr, n, ptr)
REAL, INTENT(IN), TARGET :: arr(*)
INTEGER, INTENT(IN) :: n
REAL, INTENT(OUT), POINTER :: ptr(:)
!****
ptr => arr(1:n*n:n+1)
END SUBROUTINE get_diagonal_pointer
END PROGRAM diagonal
但请注意array
主程序中是简单连续并具有 TARGET 属性。如果array
不只是连续的,那么事情就变得……复杂了。
您可以使用中间等级一指针和指针边界重新映射来完成相同的任务。适用与假设尺寸方法相同的要求。
REAL, POINTER :: tmp(:)
tmp(1:SIZE(array)) => array
ptr => tmp(::SIZE(array,1)+1)
矩阵的上三角形不是“规则的”(您想要指向的元素之间的间距不同),因此无法指向它。