Discussion :

[manaliac]

Bonjour encore !

j'ai un problème d'allocation de mémoire d'un tableau que je dois lire dans un fichier, voici mon code tout simple :

num_zo = nombre de zone lue.
Ni*Nj*Nk = nombre de lignes du tableau T
num_var = nombre de colonnes du tableau t

le tableau est lu sur Ni*Nj*Nk lignes et sur num_var colonnes.

Pour chaque num_zo, Ni, Nj,Nk varient :

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
integer, dimension(:),allocatable :: Ni,Nj,Nk,ifin
real*8, dimension(:,:),allocatable :: T
 
allocate(Ni(num_zo))
allocate(Nj(num_zo))
allocate(Nk(num_zo))
 
do i=1,num_zo
read(7,*) i1,i2,i3
Ni(i)=i1  ! On récupère les dimensions pour chaque zone
Nj(i)=i2
Nk(i)=i3
ifin(i)=Ni(i)*Nj(i)*Nk(i) ! ifin(i) determine le nombre de lignes lues dans le fichier)
!allocate(T(num_var,ifin(i)) ! c'est ici mon problème
do j=1,ifin(i)
read(7,*) T(1:num_var,j)
enddo
enddo

le problème est la taille du tableau T ou je stocke des valeurs pour chaque zone.
car quand je mets : allocate(T(num_var,ifin(i)), j'ai une erreur : attempting to allocate already allocated array

avez vous une astuce pour eviter cette erreur?

Merci

[genteur slayer]

c'est normal:
tu alloue DANS la boucle Do donc à la première itération, pas de problème, il alloue le tableau, mais dès la seconde, problème: il cherche à re-allouer ton tableau!
de plus tu peux simplifer Nj,ni,nk:

Code :

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

1
2
3
 
integer, dimension(:,:),allocatable :: N
allocate N(num_zo,3)

par exemple.

tu n'alloue pas non plus ifin
je te suggère une lecture en deux fois:
- la première pour déterminer la taille totale
- la second pour la lecture des données

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
 
integer, dimension(:,:),allocatable :: N
integer, dimension(:),allocatable :: ifin
real(kind=8), dimension(:,:,:),allocatable :: T
 
allocate(N(num_zo,3))
allocate(ifin(num_zo))
 
do i=1,num_zo
   read(7,*) N(i,1),N(i,2),N(i,3)
   ifin(i)=N(i,1)*N(i,2)*N(i,3) 
   do j=1,ifin(i)
     read(7,*)
   enddo
end do
!allocation du tableau global:
allocate(T(num_zo,numvar,maxval(ifin)))
rewind(7)
do i=1,num_zo
   read(7,*)! on saute la ligne
   do j=1,ifin(i)
     read(7,*)T(i,1:Numvar,j)
   enddo
end do

cependant cette approche consomme beaucoup de mémoire inutile surtout si tes zones sont de taille très différente: on est obligé d'alloué comme si toutes les zones avaient la même taille et donc obligé d'allouer le tableau à la taille de la plus grande zone!

il faudrait chercher le moyen de ranger plus intelligemment les données dans le tableau, genre une lecture de données brute adjoint à un tableau d'adressage!

[manaliac]

Je vais essayer de faire un rewind, ça devrait marcher, cela étant, mon programme est super moche, j'ai plein de rewind partout mais la fin justifie le moyen !!

Merci encore!

[manaliac]

en ce qui concerne l'allocation du tableau T :

allocate(T(num_zo,numvar,maxval(ifin))

pour le maxval(ifin), est ce qu'il ne vaudrait mieux pas définir une nouvelle variable, qui serait le résultat de la somme de tout les ifin pour chaque zone :

Code :

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

1
2
3
4
somme =0 ! initialisation de la variable
do i=1,ifin(num_zo)
somme= somme+ifin(i)
enddo

et puis de faire :

allocate(T(num_zo,numvar,somme)

car je voudrais avoir la taille exacte du tableau.

[genteur slayer]

je n'y avais pas pensé, tu as num_zo zones, de tailles respective ifin(
tu dois donc lire Sum(ifin) lignes de numvar variables

d'où une alloc:

Code :

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

1
2
 
allocate(T(numvar,sum(ifin))

et là, tu as exactement la bonne taille de tableau

en revanche pour acceder à tes données:

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
 
function accesdata(izone,iligne,numvar,ifin,T)
  implicit none
  integer, intent(in) :: izone,iligne,numvar
  integer,dimension(:),intent(in)::ifin
  real(kind=8),dimension(:,:),intent(in)::T
  real(kind=8),dimension(numvar)::accesdata
!
  integer debut
  if (izone==1) then
    debut=1
  else
    debut = sum(ifin(1:i-1))
  end if
  accesdata= T(:,debut+iligne)
  return
end function accesdata

afin de pouvoir accèdé à la ligne 'iligne' de la zone 'izone' dans le bloc de données T

[manaliac]

Merci de ton idée, j'avais oublié la fonction sum(valeur(i)) :arf:

Je repose la même question pour une autre lecture :



seulement le tableau que je lis, doit etre lu de la manière suivante :

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
27
28
29
30
31
32
33
 
read(7,*) (((T(num_zo,num_var_,Ni,Nj,Nk),i=1, Ni), j=1,Nk), k=1,Nk)
 
 
!mon premier problème se pose quand je veux récupérer les Ni,Nj,Nk pour !chaque num_zo (chaque zone):
! num_var designe le nombre de variable qui est pareil pour n'importe quel !num_zo
 
integer, dimension(:),allocatable :: Ni,Nj,Nk
real*8, dimension(:,:),allocatable :: T
 
allocate(Ni(num_zo))
allocate(Nj(num_zo))
allocate(Nk(num_zo))
 
do i=1,num_zo
read(7,*) i1,i2,i3
Ni(i)=i1  ! On récupère les dimensions pour chaque zone
Nj(i)=i2
Nk(i)=i3
 
! il faut lire ensuite le tableau T
 
! il faudrait un do j=1, n_lignes
!read(7,*) seulement je ne connais pas sa taille, et donc combien de lignes il !faut lire
! enfin, on peut trouver le nombre de ligne par :
!  n_ligne= (Ni(i)*Nj(i)*nk(i) )/5  ! il existe 5 valeurs de T par ligne
! quoique  c'est problematique car il faut vérifier que n_ligne est divisible par 5
! ( cas de la dernière ligne avec moins de 5 valeurs)
! enddo
 
enddo
 
allocate(T(num_zo,num_var,produit(Nbi),produit(Nbj)),produit(Nbk))

en ce qui concerne la taille des produit(Ni) et produit(Nj),je suis entrain de me casser la tête dessus !



est ce que vous auriez une idée?

[genteur slayer]

En fait c'est axactement le même tableau que précédement!!!!

Code :

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

1
2
 
allocate(T(numvar,sum(ifin))

avec la même lecture!!!

seulement, on va changer la fonction d'accès (je te laisserai faire l'éventuelle fonction d'écritue: c presque pareil )

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
function accesdata2(T,izone,ivar,i,j,k,Ni,Nj,Nk,)
  implicit none
  integer, intent(in) :: izone,ivar,i,j,k,Ni,Nj,Nk
  real(kind=8),dimension(:,:),intent(in)::T
  real(kind=8) accesdata
!
  integer debut,iligne
  integer ii,jj,kk
!
  if (izone==1) then
    debut=1
  else
    debut = sum(Ni(1:izone-1)*Nj(1:izone-1)*Nk(1:izone-1))
  end if
  iligne=(i-1)*Ni(izone)+(j-1)*Nj(izone)+k
  accesdata= T(ivar,debut+iligne)
  return
end function accesdata2