Discussion :

[Yu Cloud]

Bonjour à tous,

Je suis en train d'essayer de faire un programme qui lit les fichiers txt où les points d'une courbe sont rentrées (nous avons donc 2 colonnes de points) et ces points sont affichés à l'écran.

J'aimerais donc savoir la manière de procéder en fortran afin de lire les points qui se trouvent dans le fichier txt ?

En vous remerciant par avance,

Yu Cloud

[Marlan]

Bonjour,

Ce genre de questions a sa réponse dans la FAQ. Voici un bout de code qui te permettras de faire les choses :

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
 
PROGRAM LECT
IMPLICIT NONE
REAL(KIND=8), ALLOCATABLE,DIMENSION(:,:) :: DONNEES
INTEGER :: ios,I
 
OPEN(11,FILE='MON_FICHIER.txt',STATUS='old',ACTION='read',IOSTAT=ios)
 
I=0
DO 
READ(11,*,IOSTAT=ios)
IF(ios/=0) EXIT
I=I+1
END DO
 
REWIND(11)
 
ALLOCATE(DONNEES(I,2))
READ(11,*)DONNEES
 
PRINT'(2f15.10)', DONNEES
 
CLOSE(11)
 
DEALLOCATE(DONNEES)
END PROGRAM LECT

Le code ci-dessus part du principe que tu ne connais pas la taille de tes données (d'où le tableau dynamique). De fait, on lit d'abord le fichier à l'aide d'une variable d'erreur (ios) pour déterminer le nombre de lignes à allouer au tableau DONNEES. On alloue ensuite DONNEES avant de le lire dans le fichier puis on l'imprime à l'écran en format réel.

En espérant avoir répondu à ta question,

Si cela convient, n'hésites pas à cocher la case Résolu en bas à gauche du sujet,

Marlan

[Yu Cloud]

Mon programme fontionne maintenant !

Par contre j'ai eu un petit soucis avec "real(kind=8)" car le compilateur afficher une erreur lors de la compilation du programme, j'ai donc du le remplacer par une déclaration "double precision".

Merci pour ton aide,

Yu Could

[Marlan]

Cela revient au même, mais tu peux essayer avec REAL*8 sinon (sans trop m'avancer, il me semble que la syntaxe REAL(KIND=8) est réservée au FORTRAN90+)

[Sylvain Bergeron]

(...) Par contre j'ai eu un petit soucis avec "real(kind=8)" car le compilateur afficher une erreur lors de la compilation du programme, j'ai donc du le remplacer par une déclaration "double precision"(...)

Tu n'indique pas le message d'erreur. Il y a 2 possibilités :

1. Tu utilises un compilateur Fortran 77 alors que les kinds sont introduits par la norme Fortran 90.
2. C'est le 8 qui est refusé car ton compilateur utilise les kind numbers 1, 2, 3, ... plutôt que 4, 8, 16... Il y a au moins 1 compilateur qui a adopté cette stratégie. Les chiffres 4 et 8 ne sont pas dans la norme. Pour avoir un programme respectant la norme, il faut utiliser selected_real_kind(), mais peu de personnes le font.

Ta solution "double précision" est tout à fait correcte. La solution "real*8" est une extension Fortran 77 assez généralement supportée (mais probablement pas par le compilateur ayant adopté la numérotation des kinds 1, 2, 3...)

[Yu Cloud]

J'écris en Fortran 95 et j'utilise Silverfrost ftn95 comme compilateur Fortran.

En effet le compilateur affichait la deuxième possibilité que tu énonçais :

C'est le 8 qui est refusé car ton compilateur utilise les kind numbers 1, 2, 3, ... plutôt que 4, 8, 16... Il y a au moins 1 compilateur qui a adopté cette stratégie. Les chiffres 4 et 8 ne sont pas dans la norme. Pour avoir un programme respectant la norme, il faut utiliser selected_real_kind(), mais peu de personnes le font.

Cordialement,

Yu Cloud

[Sylvain Bergeron]

Tu devrais donc développer le réflexe d'utiliser la stratégie suivante pour tous tes développements :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
module MyKinds
   integer, parameter :: rks = selected_real_kind(6)  ! rks = real kind single
   integer, parameter :: rkd = selected_real_kind(14) ! rkd = real kind double
end module
 
program Exemple
  use MyKinds
  implicit none
  ...
  real(rks) :: x,y,z
  real(rkd) :: A,B,C
  real(rkd), parameter :: pi = 4.0_rkd * atan(1.0_rkd)
  ...

Ça semble lourd, mais le bénéfice sera appréciable à long terme.

[Yu Cloud]

Ok pour mes prochains développement j'essaierai d'utiliser ce module. Ce module permet bien de jouer sur la précision des réels que j'utilise ?

Yu Cloud

[Yu Cloud]

Bonjour,

Il y a quelques temps vous m'avez conseillé d'utiliser un module (Module MyKinds) dans mes programmes.

Est ce qu'il serait possible que vous m'expliquiez son utilité ?
Quel est la différence entre déclarer mes variables en real ou double precision et déclarer mes variables avec le module ?

En vous remerciant d'avance,

Yu Cloud

[Sylvain Bergeron]

Mieux vaut tard que jamais...

Pourquoi les kinds plutôt que real / double precision ?

Pour les raisons suivantes :

• Pour accéder aux autres précisions, quads par exemple.
• Pour développer du code générique.
• Pour simplifier le changement éventuel de précision.
• Pour contrôler la précision par « module » ou « objet ».

Quads : Je ne pense pas qu'on puisse déclarer un quad autrement qu'avec les kinds.

Code générique : Tu peux écrire une fonction de tri par exemple pour des kinds "MyKind". Au début, tu n'as qu'à définir integer, parameter :: MyKind = "mes besoins", et tout s'ajuste.

Changement éventuel de précision : Comme le code générique, on change à un endroit et tout s'ajuste.

Contrôle de précision : Tu peux prévoir une précision différente pour différentes sections d'un programme, et changer la précision de chaque section indépendamment. Par exemple, je travaille en finance. En pratique, on veut beaucoup de précision sur les données et calculs antérieurs à aujourd'hui (accumulation avec des données réelles) et une faible précision pour les projections / simulations futures. J'aurai donc des précisions définies comme suit :

Code :

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

1
2
3
4
5
6
7
8
 
   integer, parameter :: Single = selected_real_kind(6)
   integer, parameter :: Double = selected_real_kind(14)
   integer, parameter :: Quad = selected_real_kind(30)
...
   integer, parameter :: Accu_k = Double
   integer, parameter :: Proj_k = Single
...

Finalement, c'est souvent une question de complexité et d'espérance de vie du code : pour un petit programme jetable, les kinds, surtout par selected_real_kind, c'est inutilement lourd. Pour un système dont on s'attend à ce qu'il nous survive, c'est plus justifié.