Discussion :

[feynman]

Bonjour,
je suis un débutant en Fortran, j'essaye d'écrire un programme pour créer un tableau A tel que A(i)= i+2. Je suis sûr qu'on doit mettre une boucle de type do enddo mais je n'ai pas su quoi mettre dans cette boucle. Pouvez-vous me compléter mon programme ci-dessous ? Merci pour votre aide

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
program tableau 
implicit none
integer I
 
print* ,'entrer le nombre n:'
read*,n
REAL:: A(n) 
integer n
 
 
      do i=1,n
 
      enddo 
 
end

[cyrius24]

TOUJOURS LES DECLARATIONS EN PREMIER (rien avant)

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
program test 
implicit none
 
integer :: n,i
real,dimension(:),allocatable :: A
 
write(*,*) 'entrer le nombre n:'
read(*,*) n
 
allocate(A(n))
 
do i=1,n
A(i)=i+2.
end do 
 
end program test

[cyrius24]

c'est dimension( : )

[feynman]

TOUJOURS LES DECLARATIONS EN PREMIER (rien avant)

program test
implicit none

integer :: n,i
real,dimension(,allocatable :: A

write(*,*) 'entrer le nombre n:'
read(*,*) n

allocate(A(n))

do i=1,n
A(i)=i+2.
end do

end program test

Merci pour votre reponse, j ai esseye' le programme , mais il m a donne' que le nombre n, je veux faire apparaitre le tableau A, alors comment faire?

[Sylvain Bergeron]

Pour afficher le tableau A:

Code :

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

1
2
 
write(*,*)A

Ça affiche le tableau sur une ligne. Mais ce n'est pas très lisible, surtout si n est grand. Il est préférable d'afficher aussi les indices :

Code :

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

1
2
 
write(*,*) (i, A(i), i=1, n)

Ça s'affiche toujours sur une ligne, mais avec des indices. Il serait cependant préférable d'afficher un élément par ligne. Par ce faire, le plus simple est d'afficher un élément à la fois, soit 1 élément par write :

Code :

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

1
2
3
4
 
do i = 1, n
  write(*,*) i, A(i)
enddo

[feynman]

Pour afficher le tableau A:

Code :

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

1
2
 
write(*,*)A

Ça affiche le tableau sur une ligne. Mais ce n'est pas très lisible, surtout si n est grand. Il est préférable d'afficher aussi les indices :

Code :

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

1
2
 
write(*,*) (i, A(i), i=1, n)

Ça s'affiche toujours sur une ligne, mais avec des indices. Il serait cependant préférable d'afficher un élément par ligne. Par ce faire, le plus simple est d'afficher un élément à la fois, soit 1 élément par write :

Code :

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

1
2
3
4
 
do i = 1, n
  write(*,*) i, A(i)
enddo

Donc comment ca devient le program en totale si j utilise ton dernier code?

[Sylvain Bergeron]

Tu dois afficher le tableau après lui avoir assigné des valeurs.

[feynman]

C'est-à-dire comme ça :

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
program test 
implicit none
 
integer :: n,i
real,dimension(:),allocatable :: A
 
write(*,*) 'entrer le nombre n:'
read(*,*) n
 
allocate(A(n))
 
do i=1,n
A(i)=i+2.
enddo 
do i = 1, n
  write(*,*) i, A(i)
enddo
end

?
A quoi ça sert alors allocatable :: A? et est-ce qu'il y a une commande équivalente ?
Merci de nouveau

[Sylvain Bergeron]

C est a dire comme ca :
...
do i=1,n
A(i)=i+2.
enddo
do i = 1, n
write(*,*) i, A(i)
enddo
...

Tout à fait. Remarque que tu pourrais combiner les boucles comme suit :

Code :

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

1
2
3
4
5
6
...
do i = 1, n
  A(i)=i+2.
  write(*,*) i, A(i)
enddo
...

mais en pratique, dans un programme plus imposant, il faut presque toujours séparer les tâches de calculs et les tâches d'impression. La première forme est donc « meilleure ».

A quoi ca sert alors allocatable :: A? et est ce qu il y a une commande equivalente?

Ça sert à déclarer un tableau dont les dimensions ne seront déterminées qu'à l'exécution. Il n'y a pas vraiment de commande équivalente.

En F77, il n'y avait pas d'allocatable. On devait déclarer les tableaux assez grand pour supporter les plus gros cas prévus et valider si la situation était prévue:

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
program test77
implicit none
 
integer n,i,DimMaxA
parameter (DimMaxA=1000)
real A(DimMaxA)
 
write(*,*) 'entrer le nombre n:'
read(*,*) n
if (n .gt. DimMaxA) then
  write(*,*) 'n est trop grand'
  stop
endif
 
do i=1,n
A(i)=i+2.
enddo
do i = 1, n
write(*,*) i, A(i)
enddo
end

[feynman]

Ca marche pour les tableaux à une dimension, mais maintenant j'essaye de faire le tableau de 3 lignes et trois colonnes (matrice) tel que a(i,j)=i*j.
J'ai essayé mais ça n'a pas marché, pouvez-vous m'aider et dire où est mon erreur ?
Merci
Voilà mon programme :

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
program test 
implicit none
 
integer :: n,i,j
real,dimension(:,:),allocatable :: A
 
write(*,*) 'entrer le nombre n:'
read(*,*) n
 
allocate(A(n,n))
 
do i=1,n
  do j=1,n
A(i,j)=i*j
enddo
enddo 
do i = 1,n
  do j=1,n
  write(*,*) i,j, A(i,j)
enddo
enddo
end

[Sylvain Bergeron]

Je ne vois pas d'erreur. J'ai compilé et ça fonctionne.

Quelle erreur as-tu ?

[genteur slayer]

non plus je vois pas d'erreurs

une petite note supplémentaire: pour afficher ta matrice sous forme d'un tableau (que ce soit à l'écran ou dans un fichier) tu demande à WRITE de t'afficher la liste des valeur de la première ligne (FORTRAN stocke Toujours les donnée en collone):

Code :

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

1
2
3
4
 
do i=1,n
  write(*,*)(A(i,j),j=1,n)
enddo

de plus je te suggère si tu as des matrices bien plus grande de trouvé un cours sur les FORMAT en fortran permettant d'afficher que les chiffres significatifs par exemple et surtout dans l'exemple que je t'es donné d'afficher tous les nombres de la même manière et donc de parfaitement alligner les colones....

[feynman]

Merci beaucoup pour votre reponse j'ai essaye aujourd'hui avec do i=1,n
write(*,*)(A(i,j),j=1,n)
enddo
et ca marche correctement mais est-ce que je peux mettre ma matrice entre colonne, point d'exclamation (comme scilab)?
merci de nouveau

[genteur slayer]

que veux-tu dire par là.... la syntaxe que je t'es donné te l'écris en tableau donc trois ligne de trois colones, qu'appele-tu affichage en colonne?
dans le style:
1 1 5.2266
1 2 7.336
...

???
parce que dans ce cas, ta première syntaxe est la bonne

[feynman]

en effet je veux ma matrice sous cette forme(avec les points d'exlamation !) par exemple
!1 2 3!
!9 2 4!
!3 1 2!

[genteur slayer]

dans ce cas:

Code :

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

1
2
3
do i=1,n
write(*,*)'!',(A(i,j),j=1,n),'!'
enddo

après, je te dis, au lieu du "(*,*)" tu peux mettre un indicateur de format à la place de la second étoile, pour le cas avec les point d'exclamation:

Code :

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

'(a1,3i1,a1)'

seulement celui-là ne marche que pour une matrice de trois colonnes avec des entiers à un chiffre....

[feynman]

Bizare , votre premier code ne marche pas, ca me donne erreur apres la compilation

[feynman]

maintenant ca marche merci beaucoup

[feynman]

Autre question sur les tableaux,
Je vient de voire sur un cite(http://www.u-cergy.fr/sir/data/cours.../wilkinson.f90) cette expression en fortran :vect(2:n*n:n+1) = 1
j arrive pas a la comprendre. Est ce quelqu un pourrai me l expliquer?
merci de nouveau

[Ehouarn]

Bonjour,

L'expression

Code :

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

1
2
 
vect(i:j:k)

correspond à une boucle implicite sur les indices du vecteur 'vect', partant de l'indice 'i', allant jusqu'à l'indice 'j', par pas de 'k'.

L'instruction

Code :

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

1
2
 
vect(2:n*n:n+1)=1

correspond donc à la boucle suivante:

Code :

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

1
2
3
4
 
do ind=2,n*n,n+1
  vect(ind)=1
enddo