Matrice triangulaire dans un tableau 1D

Version imprimable

21/08/2008, 13h21
lguignar

Matrice triangulaire dans un tableau 1D

Bonjour,
J'ai une matrice triangulaire que j'aimerai stocker dans un tableau à une dimension.
Je me souvient que c'est possible mais je n'arrive pas à me souvenir de "l'équation" qui lie les coordonnées en 2D de ma matrice (par exemple i et j) et les coordonnées 1D de mon tableau (par exemple k)
Sachant que je peux stocker la partie inférieure ou supérieure de la matrice et que je ne souhaite pas stocker la diagonal.

Merci d'avance, Léo
21/08/2008, 15h27
pseudocode

Je pense que je l'ai déja ecrit sur ce forum...

Matrice triangulaire de taille N avec sa diagonale (row<=col)

index = row*(N-1) - TriangularNumber(row-1) + col

Matrice triangulaire de taille N sans sa diagonale (row<col)

index = row*(N-2) - TriangularNumber(row-1) + col-1

avec TriangularNumber(a) = a*(a+1)/2

Salut!
Tu peux t'inspirer de ceci (extrait du sous-programme DSPCO, de la bibliothèque LinPack disponible gratuitement sur le site www.netlib.org

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
C     PACKED STORAGE
C
C          THE FOLLOWING PROGRAM SEGMENT WILL PACK THE UPPER
C          TRIANGLE OF A SYMMETRIC MATRIX.
C
C                K = 0
C                DO 20 J = 1, N
C                   DO 10 I = 1, J
C                      K = K + 1
C                      AP(K) = A(I,J)
C             10    CONTINUE
C             20 CONTINUE
C
C     LINPACK. THIS VERSION DATED 08/14/78 .
C     CLEVE MOLER, UNIVERSITY OF NEW MEXICO, ARGONNE NATIONAL LAB.

Jean-Marc Blanc