Discussion :

[einrish]

Bonjour,

Je coince sur un petit problème en fortran !
Je ne suis pas informaticien mais physicien donc vous allez surement pouvoir m'aider
J'aurais besoin d'une structure de tableaux déchiquetés (jagged array) en fortran 95.
C'est à dire un tableau de N lignes où chaque ligne peut avoir un nombre de colonnes différent.

Est ce possible ? Efficace ?

Merci pour votre aide

Bonne journée

[Sylvain Bergeron]

Ça se fait au moyen d'un vecteur de pointeurs représentant les lignes, et où chaque pointeur contient un vecteur pour les éléments de la ligne, soit les colonnes.

La difficulté est que Fortran ne permet pas de définir un vecteur de pointeurs ! Il faut donc utiliser une astuce, au moyen d'un type utilisateur (type). Et tu ne retrouveras pas la syntaxe A(i,j) comme il est possible de faire en C par exemple...

Voici un exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
 
type TCol
   real, pointer :: Col(:)
end type
type(TCol), allocatable(:) :: Ligne
...
allocate(Ligne(1:N))
do k=1,N
   allocate(Ligne(k)%Col(1:NbColRequisePourLaLigne(k)))
enddo
...
a = f(Ligne(i)%Col(j))
...

Ce code est valide (sous réserve de coquilles) en F90 et F95. À compter de F95 + TR allocatables, donc pour les compilateurs récents, tu peux changer la ligne real, pointer :: Col(:) par real, allocatable(:) :: Col pour avoir un code un peu plus robuste (moins de risque de perte de mémoire (memory holes)) et plus efficace (il est plus facile pour un compilateur d'optimiser des allocatables que des pointeurs).

Tu pourrais aussi pousser un peu plus loin et créer des fonctions et routines pour accéder au tableau et alléger la syntaxe des références :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
type TCol
   real, pointer :: Col(:)
end type
type(TCol), allocatable(:) :: A
...
allocate(Ligne(1:N))
do k=1,N
   allocate(A(k)%Col(1:NbColRequisePourLaLigne(k)))
enddo
...
call Set(A,i,j,expression)
b = f(Get(A,i,j))
...
contains
  function Get(A,i,j)
    type(TCol), intent(in) :: A(:)
    integer, intent(in) :: i,j
    Get = A(i)%Col(j)
  end function Get
  subroutine Set(A,i,j,V)
    type(TCol), intent(inout) :: A(:)
    integer, intent(in) :: i,j
    real, intent(in) :: V
    A(i)%Col(j) = V
  end subroutine Set