递归调用出错，error #7121

bfjc

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	type octree_node

02	type(point), dimension(:), pointer :: points

03	type(octree_node), dimension(:), pointer :: children

04	real :: x_min, x_max, y_min, y_max, z_min, z_max

05	integer :: num_points

06	end type octree_node

07

08

Contains

09	recursive subroutine build_octree(node, points, x_min, x_max, &

10	y_min, y_max, z_min, z_max)

11	type(octree_node), pointer :: node

12	type(point), dimension(:), pointer :: points

13	real, intent(in) :: x_min, x_max, y_min, y_max, z_min, z_max

14	integer :: i, j, k, num_points, num_children

15	real :: x_mid, y_mid, z_mid

16

17	do k = 0, 1

18	do j = 0, 1

19	do i = 0, 1

20	num_children = num_children + 1

21	allocate(node%children(num_children))

22	call build_octree(node%children(num_children), points,&                     !这里报错，error #7121: A ptr dummy may only be argument associated with a ptr, and this array element or section does not inherit the POINTER attr from its parent array.   [CHILDREN]

23	x_min + i*(x_mid - x_min), x_min + (i+1)*(x_mid - x_min),&

24	y_min + j*(y_mid - y_min), y_min + (j+1)*(y_mid - y_min),&

25	z_min + k*(z_mid - z_min), z_min + (k+1)*(z_mid - z_min))

26

end do

27

end do

28

end do

29

30	end subroutine build_octree

这里截取了一部分问题代码，在递归调用函数时第一个变量出现了问题，如果第一个变量改为node则可以运行，是因为node%children(num_children)不是指针吗，不太明白，希望大佬们帮忙看看。

fcode

node%children 的类型虽然为  type(octree_node), dimension(:), pointer :: children 是指向数组的指针。
但 node%children(num_children) 表示 node%children 指针指向目标的第 num_children 个元素。
因为 node%children 指向的目标类型是 type(octree_node), dimension(:)，所以它的第 num_children 个元素，类型为 type(octree_node) 而不是 type(octree_node), pointer

有两种修改方法，我推荐第一种
一，第11行，由 type(octree_node), pointer :: node 修改为 type(octree_node) :: node
二，build_octree 函数中增加一个指针变量 type(octree_node), pointer :: pNode，在 call build_octree(pNode 时传入，传入之前先 pNode => node%children(num_children)

由于Fortran是默认传地址的，所以除非你要在子程序中分配指针，否则没有什么必要把虚参定义为指针。

bfjc

好的好的，谢谢刘涛。