1

http://fortranwiki.org/fortran/show/Linked+listにある Fortran で書かれた一般的なリンク リスト (genII.f90) を実装しました 。

私はそれをテストし、LI_Remove_Head 関数がメモリを解放していないように見えるという事実を除いて、すべて問題ありません。

元のモジュールに関数 LI_Destruct (以下を参照) を追加すると、同じ結果でメモリが解放されません。

SUBROUTINE LI_Destruct(List)  
implicit none
TYPE(List_Type),INTENT(INOUT),TARGET :: List
TYPE(Link_Ptr_Type) :: Link_current, Link_next

Link_next%P =>List%Head%next


do while (associated(Link_next%P))
    Link_current%P => Link_next%P

    Link_next%P => Link_next%P%next
    deallocate(Link_current%P)
end do
end subroutine LI_destruct

私は確かに何かを見逃しているので、私の質問は2つあります.1-コードにエラーはありますか? 「割り当て解除」によってメモリが空にならない理由は何ですか?

2- Fortran 用のより優れた、ほぼ標準的な汎用リンク リストが存在しますか?

テストを行うために使用する以下の簡単なコードを追加します。

PROGRAM test_list

! Defines data and other list(s) and arrays for particles. 

USE Generic_List, ONLY : Link_Ptr_Type,Link_Type,List_Type
USE Generic_List, ONLY : LI_Init_List,LI_Add_To_Head,LI_Add_To_Tail,LI_Get_Head,&
     LI_Remove_Head,LI_Get_Next,LI_Associated,LI_Get_Len, LI_destruct

IMPLICIT NONE

TYPE:: Particle_data
    REAL, dimension(2)   :: pos       !! Coordinate dimensionali 
END TYPE Particle_data 

! Definition of the types necessary for the list
TYPE Particle_Node
  TYPE(Link_Type) :: Link
  TYPE(Particle_data), pointer :: Data  
END TYPE Particle_Node

TYPE Particle_Node_ptr
  TYPE(Particle_Node), pointer :: P
END TYPE Particle_Node_ptr

! Create array of lists in order to allow classify the particles
TYPE(List_Type), allocatable      :: ao_Particle_List(:)
TYPE(Link_Ptr_Type)  :: Link
TYPE(Particle_Node_ptr)  :: Particle_elem

!-------------------------------------------------------------!
!-------------------------------------------------------------!

INTEGER, parameter :: Npart_test = 1000000  ! , nPart
INTEGER :: i,iter,j,item,nBuffer
REAL :: pos(2)

    nBuffer = 5

    IF (ALLOCATED(ao_Particle_List)) DEALLOCATE(ao_Particle_List)
    ALLOCATE(ao_Particle_List(0:nBuffer))  

    ! Init list used for temporary construction
    DO iter=0,nBuffer
        CALL LI_Init_List(ao_Particle_List(iter))
    ENDDO

    DO j=1,NBuffer
        DO i=1,Npart_test
            pos(1)=i*1.0;  pos(2)=j*i

            ALLOCATE(Particle_elem%P); ALLOCATE(Particle_elem%P%Data)   ! Allocate data before store
            Particle_elem%P%Data%pos = pos

            ! Elem is treated and should be put at head of the list ao_Particle_List(item)
            item=j
            Link = TRANSFER(Particle_elem,Link); CALL LI_Add_To_Head(Link,ao_Particle_List(item)) ! STORAGE
        END DO
    END DO

    WRITE(*,*) "List is full, see RAM"; READ(*,*)


    ! Write(*,*) "Destruct list"
    DO iter=0,nBuffer
        CALL LI_Destruct(ao_Particle_List(iter))  
    ENDDO

    IF (ALLOCATED(ao_Particle_List)) DEALLOCATE(ao_Particle_List)
    WRITE(*,*) "List is empty, see RAM";     READ(*,*)

END PROGRAM

みんなのおかげで、ジョン

4

1 に答える 1

0

問題は解決しました。転送操作があるため、完了した機能でデータに到達できないようLI_DESTRUCTです。

リンクされたリストの割り当てを解除する正しい方法は次のとおりです (fortran wiki の例とは少し異なります)。

DO
    Link = LI_Remove_Head(ao_Particle_List)
    IF(.NOT.LI_Associated(Link))EXIT
    User = TRANSFER(Link,User)
!~                 WRITE(6,*)User%P%Data%Index,User%P%Data%User_Stuff
    DEALLOCATE(User%P%Data)
    DEALLOCATE(User%P)
ENDDO
于 2012-09-04T14:08:25.053 に答える