création de matrice

Bonjour,

j'ai écrit un programme pour créer des matrices en donnant les éléments de la matrice.

C.....Déclarations
implicit none
integer i,j,k,nl1,nc1
integer nl3,nc3
integer mat1(10,10)
character*3 test
C.....Saisie des matrices
open(7,file='mat1.dat')
1000 call lect (nl1,nc1,mat1)
print*, 'donner les valeurs des éléments de la matrcice 1 '
do i=1,nl1
print*, (mat1(i,j),j=1,nc1)
end do
write(7,*) mat1
close(7)
end
C.....Subroutine de lecture de matrices
subroutine lect(NL,NC,MAT)
integer mat(10,10)
print*, 'nombres de lignes et de colonnes de la matrice'
read*, NL,NC
print*, 'donner les valeurs des ‚l‚ments de la matrice'
do i=1,NL
read*, (mat(i,j),j=1,NC)
end do
return
end

Je voudrais maintenant que les éléments de la matrice s'inscrivent automatiquement à partir d'un unique nombre k. Ce même nombre k défini le nombre de ligne et de colonne de la matrice tel que:
nl=k+1
nc=k

La matrice a cette forme:

0 0 0 0 ....
a b b b ...
b a b b ...
b b a b ....
b b b a ....

Les éléments a et b sont calculés à partir du nombre demandé.
J'ai essayé avec des boucles DO WHILE, DO et IF mais rien ne fonctionne.
Petite précision qui peut avoir son utilité, je travail en FORTRAN 77.

Si vous avez des idées, je suis preneur.

Yvan

Salut,

Il semble que tu ais besoin d'allocation dynamique pour ton programme, puisque tu ne connais pas la taille à la compilation du coup

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
 
    ! déclaration de notre matrice, on peut aussi utiliser "pointer" à la place d'allocatable
    ! 
    REAL, DIMENSION(:,:), ALLOCATABLE :: matrix
 
    ! on alloue de l'espace pour notre matrice en fonction des 
    !*paramètres lus
    ALLOCATE(matrix(nl, nc))
 
    ! on utilise la matrice
 
    ! on libère la mémoire quand on en a plus besoin
    DEALLOCATE(matrix)

En espérant que ça puisse aider,

Bon we ;)

allocate

Le problème c'est que la fonction "allocate" ne fonctionne pas en FORTRAN 77.

Mais j'ai réussi à me dépatouiller.

Yvan

salut Yvan,

la question que tu as posé concernant l'écriture de la matrice était vraiment intéressante, apparemment tu as trouvé le moyen de le faire, pourrais-tu m'en dire un peu plus.

merci d'avance

En faite, ça ne me présent pas les résultats sous forme de matrice mais de blocs de données qui correspondent à chaque ligne de la matrice.
Vu l'utilisation que j'en fait derrière, ça me suffit largement.
Voici le bloc de programmation:

C.....Cr‚ation de la matrice d'essais
print*, 'cr‚ation de la matrice d''essai'
OPEN(3,FILE='matessai.dat')
DO i=1,nbessai
WRITE(3,'(a)')' Valeurs des variables pour l''essai'
WRITE(3,*) i
DO v=1,k
j=i-1
IF(i.EQ.1)THEN
var=val_init(v)
ELSE IF(v.EQ.j)THEN
var=val_init(v)+pas(v)*p
ELSE
var=val_init(v)+pas(v)*q
END IF
WRITE(3,6) var
WRITE(*,'(2(a,X,I2,2X),2X,F8.3)') 'essai',i,'var',v,var
END DO
END DO

Tu t'en sors avec le simplex? Je présente le 26 mars...ça arrive trop vite.