Discussion :

[tourby]

Bonjour,

je débute avec le Fortran, donc ma question pourrait paraître assez simple pour certains

je dispose de 2 fichiers texte, comprenant plusieurs variables en colonnes et ce sur plusieurs lignes.
Je désire comparer les 2 fichiers afin de réaliser des opérations assez simples.
Ils se composent ainsi:

Id1 Id2 Composante
Id1 Id2 Composante
...

Je pose une condition if et .and. pour retrouver les mêmes variables d'un fichier à l'autre Id et Id2. Lorsque le programme trouve des 2 mêmes Id, il réalise une opération sur la Composante.

Cependant, mon programme ne marche pas. J'ai établi une boucle avec do while sur i pour les 2 fichiers. Ca ne marche pas, alors je pensais faire une boucle avec i pour le 1er fichier et une 2eme avec j pour le 2eme fichier afin qu'il parcourt d'abord 1 ligne pour le 1er fichier et toutes les lignes pour le 2e et ainsi de suite. Mais ça ne marche pas.

Soit j'ai fait un erreur dans mon programme genre dans la place de la boulce de j.
Soit c'est tout autre chose...

[Grame]

Bonjour,

je débute avec le Fortran, donc ma question pourrait paraître assez simple pour certains

Bah, on a tous débuté un jour

Cependant, mon programme ne marche pas. J'ai établi une boucle avec do while sur i pour les 2 fichiers. Ca ne marche pas,

Donne nous ton fichier source. On essaiera de voir ce qui ne va pas.
Précise l'OS et le compilateur (et avec quelles options) tu utilises.

[tourby]

Bonjour,


Je travaille et compile sous Linux

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
34
35
36
37
38
39
40
 
program clca
implicit none
 
real	:: T(99999),T09(99999),T10(99999)
integer :: Bdr_id(99999),Bd(99999),Bf(99999),S1(99999),S2(99999),S3(99999)
integer :: Bdr_id09(99999),Bd09(99999),Bf09(99999),S109(99999),S209(99999),S309(99999)
integer :: Bdr_id10(99999),Bd10(99999),Bf10(99999),S110(99999),S210(99999),S310(99999)
integer	:: i,j
 
i=1
j=1
 
open(unit=10,name="clca10.txt",status="old",err=1000)
open(unit=11,name="clca09.txt",status="old",err=1000)
open(unit=12,name="clca.txt",status="replace",err=1000)
 
do while(.true.)
 do while(.true.)
 
 read(10,*,end=1000) Bdr_id10(i),Bd10(i),Bf10(i),S110(i),S210(i),S310(i),T10(i) 
 
 read(11,*,end=1000) Bdr_id09(j),Bd09(j),Bf09(j),S109(j),S209(j),S309(j),T09(j)
 
   if(Bdr_id10(i)==Bdr_id09(j).and.Bd10(i)==Bd09(j))then
     T(i)=T10(i)-T09(j)
   endif
 
 
write (12,*) Bdr_id10(i),Bd10(i),Bf10(i),S110(i),S210(i),S310(i),T(i)
 
j=j+1
 enddo 
 
i=i+1
enddo
 
1000 continue
 
end program

merci

[genteur slayer]

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
34
35
36
37
38
39
 
program clca
implicit none
 
integer, parameter :: Taille=99999
 
real,dimension(Taille)       :: T09(Taille),T10(Taille)
integer,dimension(Taille,6):: B09,B10 
integer    :: i,j,N09,N10,ios
 
!lecture du fichier clca10.txt et stockage.
open(unit=10,file="clca10.txt",status="old",iostat=ios)
i=0
do while(ios==0)
  i=i+1
  read(10,*,iostat=ios) (B10(i,j),j=1,6),T10(i)
end do
N10=i
close(10)
!lecture du fichier clca09.txt et stockage
open(unit=10,file="clca09.txt",status="old",iostat=ios)
do while(ios==0)
  read(10,*,iostat=ios) (B09(i,j),j=1,6),T09(i) 
end do
N09=i
close(10)
 
open(unit=12,file="clca.txt",status="replace",iostat=ios)
do i=1,N10
 do j=1,N09
! si il y a correspondance, alors on fait l'opération
   if((B10(i,1)==B09(j,1)).and.(B10(i,2)==B09(j,2)))then
     T10(i)=T10(i)-T09(j)
   endif
 end do
 write (12,*) (B10(i,j),j=1,6),T10(i)
end do
close(12)
end program

moi je le ferai plutôt de cette manière mais cela te demande beaucoup de mémoire et ce n'est pas forcément utile:
on lit entièrement le fichier 10 et on stocke ses valeurs le compteur i permet en fin de boucle de connaitre le nombre de valeurs on fait de même pour le fichier 9
on fait ensuite une boucle pour trouver la corespondance et si oui, on éfectue l'opération, sinonon laisse le tableau T01 tranquille... puis on écrit le fichier...
cependant tu n'a pas besoin de stocker entièrement le fichier 10:

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
34
35
36
37
38
39
program clca
implicit none
 
integer, parameter :: Taille=99999
 
real,dimension(Taille)       :: T09(Taille),T10
real T10
integer,dimension(Taille,6):: B09,B10 
integer,dimension(6) :: B01
integer    :: i,j,N09,N10,ios
 
 
!lecture du fichier clca09.txt et stockage
open(unit=10,file="clca09.txt",status="old",iostat=ios)
do while(ios==0)
  read(10,*,iostat=ios) (B09(i,j),j=1,6),T09(i) 
end do
N09=i
close(10)
 
 
!lecture du fichier clca10.txt et stockage.
open(unit=10,file="clca10.txt",status="old",iostat=ios)
open(unit=12,file="clca.txt",status="replace",iostat=ios)
 
 do while(ios==0)
   read(10,*,iostat=ios) (B10(j),j=1,6),T10
  do j=1,N09
! si il y a correspondance, alors on fait l'opération
   if((B10(1)==B09(j,1)).and.(B10(2)==B09(j,2)))then
     T10=T10-T09(j)
   endif
 end do
 write (12,*) (B10(j),j=1,6),T10
end do
close(10)
close(12)
 
end program

dans ce second exemple, tu lit le fichier 9 entièrement pour sauvegarder ses valeurs en mémoire ensuite, tu lit ligne par ligne le fichier 10 que tu compare à chacune des valeurs de 9, effectue l'opération si correspondance enfin écrit la ligne dans le fichier résultat... tu économise donc pas mal de mémoire.

enfin dans ton programme à toi, ce qui est très dérangeant, c'est qu'à chaque fois que tu lit une ligne de 10, tu lit TOUTES les lignes de 9 ce qui veut dire que tu lit et relit les données de 9 en permanence ce qui est très lent puisque ton fichier est sur le disque dur... il vaut mieux le charger en mémoire c'est beaucoup plus rapide en traitement!!!

p.s. l'attribut IOSTAT renvoie un entier te disant si la lecteur ou l'ouverture (dans le open) s'est bien dérouler... il devient non nul également lorsque tu arrive à la fin du fichier... donc tant que ios est nul, on peut continuer à lire...

[tourby]

Bonjour, merci pour ta réponse précise!

Cependant, lorsque j'applique ta première solution, le programme ne semble pas calculer ma condition. Il me réécrit exactement leT10 du fichier clca10

J'ai aussi changé les open unit, tu écrivais 2 fois unit=10 pour les 2 fichiers à lire

EDIT :

Oups,

en fait il manquait 2, 3 lignes que tu devais sans doute avoir oublié, ou fait exprès pour me faire cogiter mais ça marche

Merci beaucoup

[genteur slayer]

non c 'était pas fait exprès...
mais pour les deux open unit=10, vu que tu ouvre puis ferme la 10, elle est de nouveau libre pour ouvrir un autre fichier!