Modulef
Suiv.: 1.12 La SD NDL1 Sup.: 1 Description des SD par type Préc.: 1.10 La SD MPLR Index Table des matières

1.11 La SD MUA

1.11.1 Contenu

Cette SD est l'organisation permettant de stocker, pour chaque ligne, les coefficients compris entre la première colonne de coefficient a priori non nul et la diagonale d'une matrice creuse de type "éléments finis". Ce type de stockage est dit profil ou skyline.

La SD MUA est composée de six tableaux dont l'organisation est prédéfinie.

Tableau MUA0 : Informations générales.

De type entier, ce tableau possède 32 variables. Il contient une description générale du travail (titre, date, nom), de la SD MUA (type, niveau, ...) et indique la présence ou non du tableau MUA1.

1:20 TITRE

le titre du travail en 20 mots de 4 caractères,

21:22 DATE

la date de création en 2 mots de 4 caractères,

23:28 NOMCRE

le nom du créateur en 6 mots de 4 caractères,

29 'MUA'

le type de la SD,

30 NIVEAU

le paramètre de niveau de la SD,

31 ETAT

un paramètre réservé,

32 NTACM

le nombre de tableaux supplémentaires associés à la SD
(ils seront décrits dans le tableau MUA1).

Tableau MUA1 : Descripteur des éventuels tableaux supplémentaires.

Ce tableau est analogue au tableau B1 de la SD B (voir cette SD).

Tableau MUA2 : Description générale de la matrice.

Ce tableau de type entier contient 12 valeurs.

1 NTYP

le type des coefficients de la matrice,

2 NBLOC

le nombre de blocs ou pages de la matrice,

3 NTCOL

le nombre de colonnes du plus grand bloc,

4 NMOPB

le nombre de mots nécessaires en m.c. pour stocker le plus grand bloc,

5 NMATA

le nombre de matrices sur le fichier,

6 NCODSA

le type de stockage de la matrice,

1 : matrice symétrique, seule la partie triangulaire inférieure est stockée, ligne par ligne,
0 : matrice diagonale,
-1 : matrice non symétrique, rangée ligne par ligne, le coefficient diagonal étant rangé en dernier,

7 LBDP

la plus grande différence + 1 entre les numéros de 2 noeuds d'un même élément si NCODSA est non nul, 1 si NCODSA=0, c'est à dire la largeur de la demi-bande en terme de noeuds.

8 LBDPDL

la plus grande différence + 1 entre les numéros des degrés de liberté de 2 noeuds d'un même élément si NCODSA est non nul, 1 si NCODSA=0, c'est à dire et la demi-largeur de bande en terme de degrés de liberté.

9 ND

le nombre de degrés de liberté par noeud s'il est constant, 0 sinon,

10 NOE

le nombre de noeuds,

11 NTDL

l'ordre de la matrice,

12 MUA5AR

l'assemblage est réalisé si 1, n'a pas déjà été fait si 0.

Tableau MUA3 : Pointeur sur le dernier degré de chaque bloc.

Ce tableau est de type entier, de longueur NBLOC+1, il contient :

MUA3(1) = 0,
MUA3(i+1) = le numéro du dernier degré de liberté du bloc i de la matrice,

Tableau MUA4 : Pointeur sur le coefficient diagonal de chaque ligne.

Ce tableau est de type entier et de longueur NTDL+1 ou 2 selon NCODSA, il contient :

NCODSA = 0
- MUA4(1) = 0
- MUA4(2) = NTDL,
NCODSA 0
- MUA4(1) = 0
- MUA4(i+1) = l'adresse du coefficient diagonal i.

Tableau MUA5 : Les coefficients de la matrice.

De type NTYP, ce tableau contient les coefficients de la matrice :

MUA5(k) = valeur du k-ème coefficient de la matrice (si elle n'est pas paginée),
ou MUA5(k-MUA3(i)) le k-ème coefficient du bloc i (si elle est paginée).

Outils associés

Ecriture et lecture d'une SD MUA

Cette SD est de catégorie 2. Une SD, résidant en mémoire centrale, est écrite partiellement (seuls ses 5 tableaux génériques ainsi que les éventuels tableaux associés décrits dans son tableau MUA1) sur fichier séquentiel ou en accès direct via le module SDECRI.

De même, une SD, résidant sur un fichier, est lue et reportée partiellement (ses 5 tableaux génériques ainsi que les éventuels tableaux associés décrits dans son tableau MUA1) en mémoire centrale via le module SDLECT, le tableau MUA5 est adressé seulement :

Pour une lecture, le tableau MUA5 est adressé seulement et la lecture est à effectuer par le programme via une instruction du type READ :
```
  
      READ(NFMUA) LE,(M(IAMUA5-1+I),I=1,LE)
```
ou si le fichier est en accès direct :
```
  
      READ(NFMUA,NEMUA) LE,(M(IAMUA5-1+I),I=1,LE)
```
Pour une écriture, le tableau MUA5 est à reporter via une instruction du type WRITE.