Zmienne tablice rozmiarów w Fortranie bez przydziału()

Question

Czy istnieje sposób tworzenia tablic o zmiennych rozmiarach w Fortranie na stosie? Allocate() nie działa dla mnie, ponieważ umieszcza tablicę na stercie. Może to prowadzić do problemów z równoległością (zobacz moje inne pytanie: OpenMP: poor performance of heap arrays (stack arrays work fine)). Oczywiście, pewne inteligentne zarządzanie pamięcią pozwoliłoby obejść ten problem, ale zarządzanie pamięcią w Fortranie brzmi głupio.

Zasadniczo szukam dla Fortran równowartość co następuje: C:

scanf("%d", N); 
int myarray[N]; 

do ponownego iteracyjne: Nie chcę

Integer, PARAMETER :: N=100 
Integer, Dimension(N) :: myarray 

bo ten określa rozmiar tablicy przy kompilacji czas. Nie chcę też, ponieważ umieszcza tablicę na stercie.

Pomoc bardzo doceniana. Byłem bardzo szczęśliwy z tablicami Allocatable aż do ostatniego spotkania z problemem w powyższym pytaniu. Jeśli istnieje negatywna odpowiedź na to pytanie, bardzo bym docenił link do źródła.

Edycja: patrz komentarze do odpowiedzi M.S.B. Elegancki sposób robienia tego tylko stał się możliwy w Fortran 2008, i odbywa się to w konstrukcji block.

Answer 1

Fortran może automatycznie tworzyć tablice tylko z deklaracjami przy wejściu do podprogramów, o ile wymiary są znane w czasie wykonywania ... nie wymaga to deklarowania wymiarów jako atrybutów parametru, mogą to być argumenty, np. ,

subroutine MySub (N) 

integer, intent (in) :: N 
real, dimension (N) :: array 

jest ważna. Dlaczego więc nie zdecydować o rozmiarze "N" w programie głównym lub podprogramie, a następnie wywołać inny podprogram, aby kontynuować. Prawdopodobnie przy takim podejściu tablica będzie na stosie. Jak napisał @eriktous, standard języka Fortran nie określa tego zachowania. Niektóre kompilatory przełączają lokalne tablice na stos po określonym rozmiarze. Niektóre kompilatory udostępniają opcje kontrolowania tego zachowania. Umieszczenie dużych tablic na stercie prawdopodobnie zostanie przesłonięte rekursywnie lub OpenMP.

Można również przekazać tablicę alokacji jako rzeczywisty argument podprogramu, bez podania fałszywego argumentu podprogramu jako możliwego do przydzielenia. Co może nie pomóc w twoich obawach, ponieważ pierwotna tablica nadal będzie prawdopodobnie na kupce.

Answer 2

Standard Fortran nie ma pojęcia stosu i sterty, więc będzie to zależne od implementacji (tj. Kompilatora). Patrząc na dokumenty dla np. gfortran, ma opcję

 
-frecursive 
    Allow indirect recursion by forcing all local arrays to be allocated on the stack. 

Inne kompilatory mogą mieć podobne opcje. Być może robi to, co chcesz.