Discussion :

[Bizounett]

Bonjour à tous,

Je dois transcrire des subroutines de fortran en Delphi. Tous va bien sauf une que je ne comprends pas, En voici une partie que je ne comprends pas :

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
 
*
*      ************************
*      * minver : inverse of a matrix *
*      ************************
*
       SUBROUTINE MINVER(A,N,D,L,M)
C
       IMPLICIT DOUBLE PRECISION (Z-H, O-Z)
       DIMENSION A(N*N),L(N),M(N)
C
      D=1.
      NK=-N
      DO 80 K=1,N
      NK=NK+N
      L(K)=K
      M(K)=K
      KK=NK+K
      BIGA=A(KK)
...

la variable A est un tableau de deux dimensions passer un paramètre, donc c'est quoi BIGA=A(KK), il manque une dimension ou elle est traité autrement. Si quelqu'un à une idéee ? Je dois transcrire cette subroutine assez rapidement, c'est le département d'ingénierie qui attends après cette réponse pour un calcule important.

Merci

[Fortran90]

Ici A(N*N) est un tableau unidimensionnel.

En fortran, un tableau bidimensionnel NxN peut se voir aussi comme un tableau unidimensionnel N*N

En fortran le stockage est par colonne ( à l'inverse du C ), cad que les N premiers chiffres correspondent à A(:,1), la première colonne.

Dans ton code, KK vaut successivement 1, N+2, 2N+3...cad que tu acccèdes les éléments diagonaux de la représentation bidimensionnelle de A.

[Bizounett]

OK, mais la tu me perd pas a peut près. Dans mon call de la subroutine, le A est bien un tableau a deux dimensions. Suposons que mon tableau est comme suit :

tableau(1,1) = 1
tableau(1,2) = 2
tableau(1,3) = 3

tableau(2,1) = 4
tableau(2,2) = 5
tableau(2,3) = 6

tableau(3,1) = 7
tableau(3,2) = 8
tableau(3,3) = 9

dans ma fonction minver, que vaut A(2) ? Ça pourrais m'aider a comprendre, je suis un visuel !

En fait voici la partie du code qui déclare et affecte le tableau a deux dimensions

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
 
integer*4     rank
parameter   (rank = 9)
real*8        matrix(rank,rank)
real*8        vector_w1(rank), vector_w2(rank)
real*8        determinant
 
matrix(1,1) = 0.d+00
matrix(1,2) = 0.d+00
matrix(1,3) = 0.d+00
...
 
*
* le call pour la fonction minver
*
 
call minver(matrix,rank,determinant,vector_w1,vector_w2)
...

Merci

[Fortran90]

Je me répète: fortran stocke par colonne.

En mémoire, ton tableau sera la première colonne, puis la 2e et la 3e soit

1 4 7 2 5 8 3 6 9

cad

COLONNE 1
tableau(1,1) = 1
tableau(2,1) = 4
tableau(3,1) = 7
COLONNE 2
tableau(1,2) = 2
tableau(2,2) = 5
tableau(3,2) = 8
COLONNE 3
tableau(1,3) = 3
tableau(2,3) = 6
tableau(3,3) = 9

donc A( 2 )...c'est 4.

[Fortran90]

Pour le côté "visuel", regarde les touches de ton téléphone.

En fortran, tu regardes de haut en bas puis de gauche à droite.
En C , tu regardes de gauche à droite puis de haut en bas.

Sauf si tu regardes les touches de ton clavier numérique et alors là ça serait plus de haut en bas mais de bas en haut.

Enfin. On se comprend.

[Bizounett]

Merci beaucoup, ça répond très bien à ma question maintenant.

il est possible que j'ai d'autres questions plus tard sur ce même projet, mais je partirais un autre thread.

Merci encore