Discussion :

[naruneo]

Bonjour,

Je dois exporter dans un fichier texte une matrice GAMMA dont la taille peut varier de 252x252 à 5313x5313. La dimension de la matrice est connu avant la compilation (en fait on créé 10 programmes, et pour le premier programme la matrice à une dimension 252 et pour le dernier 5313). On note dans le programme NR=252.

Mon problème est que je ne sais pas écrire la fonction WRITE de telle sorte que tous les éléments de la matrice GAMMA soit listée de telle sorte dans mon fichier texte :

a11
a12
...
a1n
a21
...
ann

Code :

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

1
2
 
OPEN(11,FILE='Gamma.dat',STATUS='UNKNOWN')

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
        do 20 i=1,NR
	do 21 j=1,i
 
	GAMMA(i,j)=0
 
	Gamma(i,j)= **Je vous épargne la formule de Gamma**
 
	Gamma(j,i)=Gamma(i,j) **C'est une matrice symétrique**
 
21	continue
20	continue
 
	WRITE(11,'(F8.5)') ,Gamma

Je précise que ce n'est tout simplement que la ligne WRITE(11,'(F8.5)'),Gamma qui me pose problème, puisque que je me sert de la matrice gamma pour calculer les valeurs propres et ça marche parfaitement. Par contre le fichier gamma.dat contient plus de 250 000 valeurs alors qu'il devrait en contenir 252*252=63 504

Merci de votre aide

[Ehouarn]

Bonjour,

Si tu écris

Code :

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

1
2
 
	WRITE(11,'(F8.5)') Gamma

(tu ne devrais pas mettre de virgule entre le 'write(...)' et 'Gamma') et que 'Gamma' est un tableau, tu demande implicitement à écrire tous les éléments, à la suite les uns des autres, et en principe sur une seule ligne. Mais sur bon nombre de systèmes (tous?) les lignes ont des tailles maximales et pour de grandes matrices le résultats se retrouverait possiblement tronqué et sur plusieurs lignes.

Si tu souhaites n'écrire qu'un élément par ligne, il faut le faire explicitement, en passant par des boucles:

Code :

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

1
2
3
4
5
6
 
       do i=1,NR
          do j=1,NR
            write(11,'(F8.5)') Gamma(i,j)
          enddo
       enddo

Bonne continuation.

[naruneo]

Bien, je te remercie beaucoup, je vais essayer ça tout de suite.