Jak przekazywać alokowane tablice do podprogramów w Fortranie

Question

Następujący kod zwraca błąd segmentacji, ponieważ przydzielana tablica, którą próbuję przekazać, nie jest poprawnie rozpoznawana (rozmiar zwraca 1, kiedy powinien być 3). Na tej stronie (http://www.eng-tips.com/viewthread.cfm?qid=170599) podobny przykład wydaje się wskazywać, że powinien działać poprawnie w F95; mój plik kodu ma rozszerzenie .F90, ale próbowałem go zmienić na F95 i używam gfortran do kompilacji.

Zgaduję, że problem powinien polegać na tym, że przekazuję tablicę alokacji do podprogramu; Co ja robię źle?

!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 
PROGRAM test 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 
IMPLICIT NONE 
DOUBLE PRECISION,ALLOCATABLE :: Array(:,:) 
INTEGER      :: iii,jjj 

ALLOCATE(Array(3,3)) 
DO iii=1,3 
DO jjj=1,3 
    Array(iii,jjj)=iii+jjj 
    PRINT*,Array(iii,jjj) 
ENDDO 
ENDDO 
CALL Subtest(Array) 

END PROGRAM 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 
SUBROUTINE Subtest(Array) 
DOUBLE PRECISION,ALLOCATABLE,INTENT(IN) :: Array(:,:) 
INTEGER         :: iii,jjj 

PRINT*,SIZE(Array,1),SIZE(Array,2) 
DO iii=1,SIZE(Array,1) 
DO jjj=1,SIZE(Array,2) 
    PRINT*,Array(iii,jjj) 
ENDDO 
ENDDO 

END SUBROUTINE 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 

Answer 1

Jeśli dana procedura ma fałszywy argument, który można alokować, wymagany jest jawny interfejs w dowolnym zasięgu wywoływania.

(Istnieje wiele rzeczy, które wymagają wyraźnego interfejs, allocatable manekin jest tylko jeden).

Można przewidywać, że wyraźny sprzęgowi się poprzez umieszczenie bloku interfejsu dla podprogramu wewnątrz głównego programu. Alternatywą i daleko, o wiele lepszym rozwiązaniem jest umieszczenie podprogramu wewnątrz modułu, a następnie użycie tego modułu w programie głównym - automatycznie tworzony jest jawny interfejs. Jest taki przykład na stronie Eng-tips, do której podałeś link - patrz post przez xwb.

Zauważ, że to ma sens tylko do manekina argumentu mieć allocatable atrybut Jeśli zamierzasz robić coś związanego z jego statusu alokacji - sprawdzanie stanu, przesunięciu go, zwalnianie go itp

Answer 2

Proszę również zwróć uwagę, że Twój przydzielający atrapowy argument array jest zadeklarowany jako intent(in), co oznacza, że ​​jego status przydziału będzie zgodny z faktycznym argumentem (i nie może zostać zmieniony podczas procedury). Faktyczny argument przekazany do twojego podprogramu może być nieprzydzielony, a zatem niedozwolony do odniesienia, nawet z jawnym interfejsem. Kompilator nie będzie o tym wiedział, a zachowanie zapytań takich jak size jest w takich przypadkach niezdefiniowane.

W związku z tym należy najpierw sprawdzić stan alokacji array za pomocą allocated(array) przed odniesieniem do jego zawartości. Chciałbym dalej proponuję wdrożyć pętle na całej tablicy z lbound i ubound, ponieważ w ogóle nie mogą być pewni o granicach array Add

subroutine subtest(array) 
    double precision, allocatable, intent(in) :: array(:,:) 
    integer         :: iii, jjj 

    if(allocated(array)) then 
    print*, size(array, 1), size(array, 2) 
    do iii = lbound(array, 1), ubound(array, 1) 
     do jjj = lbound(array, 2), ubound(array, 2) 
     print*, array(iii,jjj) 
     enddo 
    enddo 
    endif 
end subroutine 

Answer 3

Jest to prosty przykład, który używa allocatable manekina z argumentów moduł.