appeler une matrice dans une subroutine

Code:

---

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

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:

---

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

---

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é,