2番目の派生型の配列へのポインターを持つ派生型があります
TYPE vertex
REAL :: x, y, z
END TYPE
TYPE path
TYPE(vertex), DIMENSION(:), POINTER :: vertices => NULL()
END TYPE
目的は、頂点配列のサイズを変更できるようにして、任意の数の頂点ポイントを配列に追加できるようにすることです。そのポインタに頂点を追加するコードを作成しました。
SUBROUTINE path_append_vertex(this, x, y, z)
IMPLICIT NONE
!-------------------------------------------------------------------------------
! Variable declarations.
!-------------------------------------------------------------------------------
TYPE(path), INTENT(inout) :: this
REAL, INTENT(in) :: x, y, z
!-------------------------------------------------------------------------------
! Local Variable declarations.
!-------------------------------------------------------------------------------
INTEGER :: status
TYPE(vertex), DIMENSION(:), ALLOCATABLE :: vertices
!-------------------------------------------------------------------------------
! Start of executable code
!-------------------------------------------------------------------------------
IF (ASSOCIATED(this%vertices)) THEN
! Create a temporary array the same size as current number of vertices. Copy the
! contents of the old array to the new array then delete the old array.
ALLOCATE(vertices(SIZE(this%vertices)), STAT = status)
CALL check_status(status)
vertices = this%vertices
DEALLOCATE(this%vertices)
! Create a new array with one extra element. Copy the contents of the temporary
! array to the new one the delete the temporary array.
ALLOCATE(this%vertices(SIZE(vertices) + 1), STAT = status)
CALL check_status(status)
this%vertices(1:SIZE(vertices)) = vertices
DEALLOCATE(vertices)
this%vertices(SIZE(this%vertices))%x = x
this%vertices(SIZE(this%vertices))%y = y
this%vertices(SIZE(this%vertices))%z = z
ELSE
ALLOCATE(this%vertices(1), STAT = status)
CALL check_status(status)
this%vertices(1)%x = x
this%vertices(1)%y = y
this%vertices(1)%z = z
ENDIF
END SUBROUTINE
いくつかのパスオブジェクトを作成します。
TYPE ipch_desc
...
TYPE(path) :: chordPath
END TYPE
SUBROUTINE ipch_desc_construct(this, ...)
...
TYPE (ipch_desc), INTENT(inout) :: this
...
! Must NULL out the vertices array or else it will point to the last
! integration_path created in memory. Not sure why these are defaulting
! to NULL
this%chordPath%vertices => NULL()
CALL path_append_vertex(this%chordPath, xcart_i(1), xcart_i(2), xcart_i(3))
CALL path_append_vertex(this%chordPath, xcart_f(1), xcart_f(2), xcart_f(3))
! Check the value of the path vertices.
write(*,*) this%chordPath%vertices
END SUBROUTINE
すべてがうまく機能し、各頂点の正しい値を取得します。たとえば、作成された3つのパスオブジェクトの場合、
-0.33808113528699218 1.0467574437103653 0.10713720000000000 -0.16057879084545851 0.49717960298733294 0.10713720000000000
-0.33322243268266594 1.0483142707971911 1.42240000000000010E-003 -0.14945358419461796 0.47017940500485894 1.42240000000000010E-003
-0.33656460666251325 1.0472460386853264 -0.10629900000000000 -0.15821659220752302 0.49230280357365630 -0.10629900000000000
コードの後半でこれらのパスオブジェクトを使用する場合、
SUBROUTINE ipch_mc_model_compute(a_ipch, ...)
...
TYPE (ipch_desc), INTENT (inout) :: a_ipch
...
! Check the value of the path vertices again.
write(*,*) a_ipch%chordPath%vertices
...
END SUBROUTINE
最初のN-1値のみが正しいままです。上で作成したものと同じ値の場合、
-0.33808113528699218 1.0467574437103653 0.10713720000000000 -0.16057879084545851 0.49717960298733294 0.10713720000000000
-0.33322243268266594 1.0483142707971911 1.42240000000000010E-003 -0.14945358419461796 0.47017940500485894 1.42240000000000010E-003
0.15094203233057696 6.94277920927416864E-310 -0.10629900000000000 1.63041663127611360E-322 3.01884064661153912E-003 6.94277920927179713E-310
作成するオブジェクトの数に関係なくpath、N番目は常に間違った値になります。これを引き起こしている可能性がありますか？