Discussion :

[brahim20]

Code :

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

1
2
3
4
5
6
7
8
9
10
 
---------programme principal--------
do i=1,n
 
   do j=1,m
        !génération d'une matrice m(i,j)
   end do
 
end do
---------------------------------------

je veux utiliser cette matrice dans une subroutine, comment faire pour la stocker et l'appeler à partir de cette dernière

[Ladgalen]

Bonjour

Est ce que tu sais transférer des variables à un sous programmes ? Que ce soit des matrices ou des variables plus simples, scalaires ?

Regardes ici :

http://fortran.developpez.com/faq/?page=section8

Tu as plusieurs solutions :
1) passages dans les arguments
2) utilisation d'un sous programmes interne au programme principal (partage des variables)
3) utilisation d'un module

Le plus simple si tu débutes c'est d'utiliser des arguments :

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
 
subroutine maroutine( ndim, matrice )
 
  implicit none
 
  ! arguments
  integer, intent(in) :: ndim
  double precision, dimension( ndim, ndim), intent(inout) :: matrice
 
  ! variables locales
  integer :: i,j
 
  do i = 1 , ndim
    do j = 1, ndim
      matrice(i,j) = dble(i + j)
    end do
  end do
end subroutine

[Marlan]

Bonjour,

Pour que la matrice créée dans la subroutine puisse être appelée dans le programme principal, il est nécessaire de la définir comme variable sortante (à l'aide de l'instruction INTENT(OUT)). Soit, lors de la déclaration de variables dans ta subroutine, mettre quelque chose sous la forme:

SUBROUTINE Matrice(A) !A est la matrice souhaitée
REAL,DIMENSION(:,: ), INTENT(OUT):: A !le type OUT définit une variable sortante, IN une entrante, INOUT une entrante ET sortante

...
END SUBROUTINE Matrice

Il faut cela étant que la matrice A soit définit dans le programme principal, au même titre que sa taille. Par exemple, dans le programme principal:
REAL, DIMENSION(3,3):: A

Dans ton cas, il serait nécessaire de définir un tableau dynamique pour A (du fait que tu sembles faire intervenir des indices n et m variables), soit dans la déclaration:

REAL, ALLOCATABLE, DIMENSION(:,: )::A

Puis dans la suite du programme:
ALLOCATE(A(n,m)) !n: indice de colonne, m indice de ligne

Et à la fin du programme (ou quand tu n'as plus besoin de A), mettre la ligne:
DEALLOCATE(A)
En espérant que cela t'ai aidé,