![[BIG]](../icons/zoom18.gif){kind=icon}
![[Normal]](../icons/zoom14.gif){kind=icon}
![[small]](../icons/zoom10.gif){kind=icon}
Suiv.: Le problème modèle non symétrique Sup.: 1 Exemples d'apprentissage Préc.: 1.1 Introduction Index Table des matières
Suiv.: Le problème modèle non symétrique
Sup.: 1 Exemples d'apprentissage
Préc.: 1.1 Introduction
Index
Table des matières
L'exemple choisi est la résolution par la méthode des éléments finis d'un problème d'élasticité linéaire assez réaliste : il s'agit d'un crochet encastré sur une partie de sa frontière, et soumis à une force de traction verticale sur la partie utilitaire. Etant données les conditions aux limites envisagées, et le matériau étant supposé homogène isotrope, le problème admet un plan de symétrie, et le domaine est réduit à un demi-crochet (voir la Figure 1).
Pour plus de détails sur la formulation variationnelle de ce problème, ainsi que sur son approximation par la méthode des éléments finis et différents aspects de la mise en oeuvre sur ordinateur, voir aussi [1]. On trouvera aussi dans cette référence d'autres expériences numériques concernant la résolution de ce problème.
Au cours de la mise en oeuvre, plusieurs étapes distinctes apparaissent avant la résolution proprement dite :
(structure de données NOPO)
par le module APNOPO,
Se reporter aux guides correspondants pour l'utilisation de ces modules.
Figure 1.1: Maillage du domaine
La géométrie du problème modèle est décrite par une structure de données NOPO, et une vue de l'objet est présenté sur la Figure 1.1. Le maillage comprend 1006 éléments et 937 noeuds, les éléments finis utilisés étant de degré 1.
En chacun des noeuds, il faut calculer les trois composantes du déplacement, le système linéaire résultant a donc 2811 inconnues. La structure de la matrice du système linéaire est représentée sur la Figure 1.2.
Figure 1.2: Structure de la matrice
Pour résoudre le système linéaire associé au problème modèle, il faut maintenant
On a déjà signalé que chacune de ces étapes était réalisée par un module distinct; pour simplifier l'utilisation du code Modulef, des enchaînements de ces modules sont programmés pour chaque option disponible (ce sont les préprocesseurs de la bibliothèque PPAL). L'utilisateur peut également créer ses propres enchaînements à partir des programmes Modulef.
Nous allons maintenant passer en revue les différentes possibilités :
Le module CHOLXX appelle successivement les modules PREPAC, ASSMUA, ASEMBV, CHOLPC et enfin DRCHPC. Voici la liste des données de CHOLXX :
mail $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD MAIL
1 5 3 $ NDSM NTYP ND
tae $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD TAE
1 $ 1 SI BDCL EST UTILISE , 0 SINON
bdcl $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD BDCL
0 $ 1 SI CL. EN RL. EXISTE
b $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD B
3 $ IMPREB
et les résultats correspondant à une exécution sur une station de travail Hewlett-Packard série 400 :
M M OOO DDDD U U L EEEEE FFFFF
MM MM O O D D U U L E F
M M M O O D D U U L EEEE FFFF
M M O O D D U U L E F
M M OOO DDDD UUU LLLLL EEEEE F VERSION 92
DATE : 01/06/92
AUTEUR : joly
**********************************************
APPEL DE LA RESOLUTION PAR METHODE DE CHOLESKY
**********************************************
-- CREATION DU DATA =-= EXECUTION MODULE (DATA EXISTANT) =-= FIN
- - -
E
-- NOM DU FICHIER CONTENANT LES DONNEES ?
data/chol
++ OPEN(10,FILE='data/chol',SPEC='OLD',RECL=0)
-- PARAMETRE D'IMPRESSION POUR L'EXECUTION ?
1
++ OPEN(11,FILE='mail',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(12,FILE='tae',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(13,FILE='bdcl',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(14,FILE='b',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAC :
NOMBRE DE MOTS DE LA MATRICE (LMUA5) : 737202
DIFFERENCE MAX ENTRE 2 NOEUDS D'UN ELEMENT (LBDP) : 70
DIFFERENCE MAX ENTRE 2 D.L. D'UN ELEMENT (LBDPDL) : 210
NOMBRE MAX DE NOEUDS D'UN ELEMENT (NNOMAX) : 8
NOMBRE DE SECONDS MEMBRES (NDSM) : 1
NOMBRE DE MOTS EN M.C. POUR LA S.D. NDL1 (MCNDL1) : 0
NOMBRE DE MOTS EN M.C. POUR LA S.D. MUA (MCMUA) : 740060
NOMBRE DE MOTS EN M.C. POUR LA S.D. B (MCB) : 5667
FIN DU MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSMUA :
FIN DU MODULE ASSMUA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
FIN DU MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMPC :
FIN DU MODULE CLIMPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CHOLPC :
** NOMBRE D'INCONNUES 2811
** LONGUEUR DU TABLEAU MUA5 737202
** TEMPS CALCUL FACTORISATION 0.565175E+02 S.
FIN DU MODULE CHOLPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE DRCHPC :
** TEMPS CALCUL RESOLUTION 0.316320E+01 S.
FIN DU MODULE DRCHPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
IMPRESSION DE LA S.D. B DE NIVEAU 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITRE :
DATE ET NOM UTILISATEUR : 01/06/92 joly
TYPE DE LA STRUCTURE DE DONNEES : B
NIVEAU ET NUMERO D'ETAT : 1 2
NOMBRE DE TABLEAUX ASSOCIES : 0
TABLEAU B 2
------------
TYPE DU TABLEAU (NTYT) : 5
NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2
TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1
NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 1
NOMBRE DE NOEUDS (NOE) : 937
NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1
NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3
NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 2811
CODE DE STOCKAGE DE B (NCODSB) : -1
TABLEAU B 3
------------
POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B
1 0 2 2811
Le module CROUXX appelle successivement les modules PREPAC, ASSMUA, ASEMBV, CROUPC et enfin DRCRPC. Voici la liste des données de CROUXX :
mail $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD MAIL
1 5 3 $ NDSM NTYP ND
tae $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD TAE
1 $ 1 SI BDCL EST UTILISE , 0 SINON
bdcl $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD BDCL
0 $ 1 SI CL. EN RL. EXISTE
b $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD B
3 $ IMPREB
et les résultats correspondants :
M M OOO DDDD U U L EEEEE FFFFF
MM MM O O D D U U L E F
M M M O O D D U U L EEEE FFFF
M M O O D D U U L E F
M M OOO DDDD UUU LLLLL EEEEE F VERSION 92
DATE : 01/06/92
AUTEUR : joly
*******************************************
APPEL DE LA RESOLUTION PAR METHODE DE CROUT
*******************************************
-- CREATION DU DATA =-= EXECUTION MODULE (DATA EXISTANT) =-= FIN
- - -
E
-- NOM DU FICHIER CONTENANT LES DONNEES ?
data/crout
++ OPEN(10,FILE='data/crout',SPEC='OLD',RECL=0)
-- PARAMETRE D'IMPRESSION POUR L'EXECUTION ?
1
++ OPEN(11,FILE='mail',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(12,FILE='tae',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(13,FILE='bdcl',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(14,FILE='b',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAC :
NOMBRE DE MOTS DE LA MATRICE (LMUA5) : 737202
DIFFERENCE MAX ENTRE 2 NOEUDS D'UN ELEMENT (LBDP) : 70
DIFFERENCE MAX ENTRE 2 D.L. D'UN ELEMENT (LBDPDL) : 210
NOMBRE MAX DE NOEUDS D'UN ELEMENT (NNOMAX) : 8
NOMBRE DE SECONDS MEMBRES (NDSM) : 1
NOMBRE DE MOTS EN M.C. POUR LA S.D. NDL1 (MCNDL1) : 0
NOMBRE DE MOTS EN M.C. POUR LA S.D. MUA (MCMUA) : 740060
NOMBRE DE MOTS EN M.C. POUR LA S.D. B (MCB) : 5667
FIN DU MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSMUA :
FIN DU MODULE ASSMUA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
FIN DU MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMPC :
FIN DU MODULE CLIMPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CROUPC :
** NOMBRE D'INCONNUES 2811
** LONGUEUR DU TABLEAU MUA5 737202
** TEMPS CALCUL FACTORISATION 0.917402E+02 S.
FIN DU MODULE CROUPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE DRCRPC :
** TEMPS CALCUL RESOLUTION 0.241842E+01 S.
FIN DU MODULE DRCRPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
IMPRESSION DE LA S.D. B DE NIVEAU 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITRE :
DATE ET NOM UTILISATEUR : 01/06/92 joly
TYPE DE LA STRUCTURE DE DONNEES : B
NIVEAU ET NUMERO D'ETAT : 1 2
NOMBRE DE TABLEAUX ASSOCIES : 0
TABLEAU B 2
------------
TYPE DU TABLEAU (NTYT) : 5
NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2
TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1
NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 1
NOMBRE DE NOEUDS (NOE) : 937
NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1
NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3
NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 2811
CODE DE STOCKAGE DE B (NCODSB) : -1
TABLEAU B 3
------------
POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B
1 0 2 2811
La place mémoire nécessaire au stockage de cette matrice étant assez importante, il est intéressant d'utiliser le module CHOLPS avec stockage de la matrice sur une mémoire secondaire. Le module CHOSXX appelle successivement les modules PREPAC, ASMAPS, ASMBMS CHOLPS et enfin DRCHPS. Voici la liste des données de CHOSXX :
mail $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD MAIL
1 5 3 $ NDSM NTYP ND
tae $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD TAE
1 $ 1 SI BDCL EST UTILISE , 0 SINON
bdcl $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD BDCL
0 $ 1 SI CL. EN RL. EXISTE
b $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD B
3 $ IMPREB
4096 $ LPAGE
et les résultats correspondants :
M M OOO DDDD U U L EEEEE FFFFF
MM MM O O D D U U L E F
M M M O O D D U U L EEEE FFFF
M M O O D D U U L E F
M M OOO DDDD UUU LLLLL EEEEE F VERSION 92
DATE : 01/06/92
AUTEUR : joly
**********************************************
APPEL DE LA RESOLUTION PAR METHODE DE CHOLESKY
**********************************************
-- CREATION DU DATA =-= EXECUTION MODULE (DATA EXISTANT) =-= FIN
- - -
E
-- NOM DU FICHIER CONTENANT LES DONNEES ?
data/chos
++ OPEN(10,FILE='data/chos',SPEC='OLD',RECL=0)
-- PARAMETRE D'IMPRESSION POUR L'EXECUTION ?
1
++ OPEN(11,FILE='mail',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(12,FILE='tae',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(13,FILE='bdcl',SPEC='OLD,UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAC :
NOMBRE DE MOTS DE LA MATRICE (LMUA5) : 737202
DIFFERENCE MAX ENTRE 2 NOEUDS D'UN ELEMENT (LBDP) : 70
DIFFERENCE MAX ENTRE 2 D.L. D'UN ELEMENT (LBDPDL) : 210
NOMBRE MAX DE NOEUDS D'UN ELEMENT (NNOMAX) : 8
NOMBRE DE SECONDS MEMBRES (NDSM) : 1
NOMBRE DE MOTS EN M.C. POUR LA S.D. NDL1 (MCNDL1) : 0
NOMBRE DE MOTS EN M.C. POUR LA S.D. MUA (MCMUA) : 740060
NOMBRE DE MOTS EN M.C. POUR LA S.D. B (MCB) : 5667
FIN DU MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
++ OPEN(14,FILE='b',SPEC='DIRECT',RECL=16384)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAC :
NOMBRE DE MOTS DE LA MATRICE (LMUA5) : 737202
DIFFERENCE MAX ENTRE 2 NOEUDS D'UN ELEMENT (LBDP) : 70
DIFFERENCE MAX ENTRE 2 D.L. D'UN ELEMENT (LBDPDL) : 210
NOMBRE MAX DE NOEUDS D'UN ELEMENT (NNOMAX) : 8
NOMBRE DE SECONDS MEMBRES (NDSM) : 1
NOMBRE DE MOTS EN M.C. POUR LA S.D. NDL1 (MCNDL1) : 0
NOMBRE DE MOTS EN M.C. POUR LA S.D. MUA (MCMUA) : 740247
NOMBRE DE MOTS EN M.C. POUR LA S.D. B (MCB) : 5667
FIN DU MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASMAPS :
SP INITMP : NOMBRE DE PAGES EN MEMOIRE CENTRALE : 17
FIN DU MODULE ASMAPS
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASMBMS :
SP INITMP : NOMBRE DE PAGES EN MEMOIRE CENTRALE : 2
FIN DU MODULE ASMBMS
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMPS :
FIN DU MODULE CLIMPS
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CHOLPS :
SP INITMP : NOMBRE DE PAGES EN MEMOIRE CENTRALE : 16
** TEMPS CALCUL FACTORISATION 0.131826E+03 S.
FIN DU MODULE CHOLPS
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE DRCHPS :
SP INITMP : NOMBRE DE PAGES EN MEMOIRE CENTRALE : 2
** TEMPS CALCUL RESOLUTION 0.266093E+02 S.
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
IMPRESSION DE LA S.D. B DE NIVEAU 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITRE :
DATE ET NOM UTILISATEUR : 01/06/92 joly
TYPE DE LA STRUCTURE DE DONNEES : B
NIVEAU ET NUMERO D'ETAT : 1 2
NOMBRE DE TABLEAUX ASSOCIES : 0
TABLEAU B 2
------------
TYPE DU TABLEAU (NTYT) : 5
NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2
TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1
NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 2
NOMBRE DE NOEUDS (NOE) : 937
NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1
NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3
NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 2811
CODE DE STOCKAGE DE B (NCODSB) : -1
TABLEAU B 3
------------
POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B
1 0 2 2048 3 2811
L'optimisation de la place mémoire est réalisée par la méthode frontale qui gère au mieux cette ressource, en ne stockant que les lignes indispensables du système linéaire. Le module FRONXX appelle successivement les modules PREPAF et FRONT. Voici la liste des données de FRONXX :
mail $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD MAIL
1 5 3 $ NDSM NTYP ND
tae $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD TAE
1 $ 1 SI BDCL EST UTILISE , 0 SINON
bdcl $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD BDCL
0 $ 1 SI CL. EN RL. EXISTE
b $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD B
3 $ IMPREB
2 $ MAJORANT DU NOMBRE DE FRONTS
et les résultats correspondants. Les fichiers fron et gaus sont des fichiers temporaires gérés par les modules PREPAF et FRONT.
M M OOO DDDD U U L EEEEE FFFFF
MM MM O O D D U U L E F
M M M O O D D U U L EEEE FFFF
M M O O D D U U L E F
M M OOO DDDD UUU LLLLL EEEEE F VERSION 92
DATE : 01/06/92
AUTEUR : joly
*******************************************
APPEL DE LA RESOLUTION PAR METHODE FRONTALE
*******************************************
-- CREATION DU DATA =-= EXECUTION MODULE (DATA EXISTANT) =-= FIN
- - -
E
-- NOM DU FICHIER CONTENANT LES DONNEES ?
data/fron
++ OPEN(10,FILE='data/fron',SPEC='OLD',RECL=0)
-- PARAMETRE D'IMPRESSION POUR L'EXECUTION ?
1
++ OPEN(11,FILE='mail',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(12,FILE='tae',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(13,FILE='bdcl',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(14,FILE='b',SPEC='UNFORMATTED',RECL=0)
++ OPEN(50,FILE='fron',SPEC='UNFORMATTED',RECL=0)
++ OPEN(51,FILE='gaus',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAF :
** NOMBRE D'INCONNUES 2811
** LONGUEUR DU FRONT 68
** TEMPS CALCUL PREPARATION 0.981024E+00 S.
FIN DU MODULE PREPAF
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE FRONT :
** TEMPS CALCUL RESOLUTION 0.127254E+03 S.
FIN DU MODULE FRONT
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
IMPRESSION DE LA S.D. B DE NIVEAU 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITRE :
DATE ET NOM UTILISATEUR : 01/06/92 joly
TYPE DE LA STRUCTURE DE DONNEES : B
NIVEAU ET NUMERO D'ETAT : 1 0
NOMBRE DE TABLEAUX ASSOCIES : 0
TABLEAU B 2
------------
TYPE DU TABLEAU (NTYT) : 5
NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2
TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1
NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 1
NOMBRE DE NOEUDS (NOE) : 937
NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1
NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3
NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 2811
CODE DE STOCKAGE DE B (NCODSB) : -1
TABLEAU B 3
------------
POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B
1 0 2 2811
Si on veut maintenant résoudre le problème modèle par la méthode du gradient conjugué, il faut assembler la matrice sous la forme d'une structure de données AMAT. Le module GRADXX appelle successivement les modules PREPGC, ASSAMA, ASEMBV, puis suivant le type de préconditionnement choisi, l'un des modules suivants :
Le problème modèle est d'abord résolu par itérations de gradient conjugué sans préconditionnement (module SIMPGC), avec les données suivantes :
nopo $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD NOPO
1 5 3 $ NDSM NTYP ND
tae $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD TAE
1 $ 1 SI BDCL EXISTE
bdcl $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD BDCL
0 $ 1 SI CL. EN RL. EXISTE
b $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD B
3 $ IMPRESSION
0 $ METHODE
0.1000000E-03 $ EPS
et les résultats correspondants :
M M OOO DDDD U U L EEEEE FFFFF
MM MM O O D D U U L E F
M M M O O D D U U L EEEE FFFF
M M O O D D U U L E F
M M OOO DDDD UUU LLLLL EEEEE F VERSION 92
DATE : 01/06/92
AUTEUR : joly
*******************************************
APPEL DE LA RESOLUTION PAR GRADIENT CONJUGUE
*******************************************
-- CREATION DU DATA =-= EXECUTION MODULE (DATA EXISTANT) =-= FIN
- - -
E
-- NOM DU FICHIER CONTENANT LES DONNEES ?
data/simp
++ OPEN(10,FILE='data/simp',SPEC='OLD',RECL=0)
-- PARAMETRE D'IMPRESSION POUR L'EXECUTION ?
1
++ OPEN(11,FILE='nopo',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(12,FILE='tae',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(13,FILE='bdcl',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(14,FILE='b',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPGC :
FIN DU MODULE PREPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSAMA :
FIN DU MODULE ASSAMA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
FIN DU MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMGC :
FIN DU MODULE CLIMGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE SIMPGC :
INITIALISATION DU G.C. A ZERO
NOMBRE D'ITERATIONS DU SYSTEME NUMERO 1 : 1212
** TEMPS RESOLUTION 0.368821E+03 S.
FIN DU MODULE SIMPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
IMPRESSION DE LA S.D. B DE NIVEAU 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITRE :
DATE ET NOM UTILISATEUR : 01/06/92 joly
TYPE DE LA STRUCTURE DE DONNEES : B
NIVEAU ET NUMERO D'ETAT : 1 2
NOMBRE DE TABLEAUX ASSOCIES : 0
TABLEAU B 2
------------
TYPE DU TABLEAU (NTYT) : 5
NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2
TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1
NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 1
NOMBRE DE NOEUDS (NOE) : 937
NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1
NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3
NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 2811
CODE DE STOCKAGE DE B (NCODSB) : -1
TABLEAU B 3
------------
POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B
1 0 2 2811
Le problème modèle est ensuite résolu par itérations de gradient conjugué avec préconditionnement SSOR (module SSORGC), avec les données suivantes :
nopo $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD NOPO
1 5 3 $ NDSM NTYP ND
tae $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD TAE
1 $ 1 SI BDCL EXISTE
bdcl $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD BDCL
0 $ 1 SI CL. EN RL. EXISTE
b $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD B
3 $ IMPRESSION
1 $ METHODE
0.1000000E-03 1.000000 $ EPS OMEGA
et les résultats correspondants :
M M OOO DDDD U U L EEEEE FFFFF
MM MM O O D D U U L E F
M M M O O D D U U L EEEE FFFF
M M O O D D U U L E F
M M OOO DDDD UUU LLLLL EEEEE F VERSION 92
DATE : 01/06/92
AUTEUR : joly
*******************************************
APPEL DE LA RESOLUTION PAR GRADIENT CONJUGUE
*******************************************
-- CREATION DU DATA =-= EXECUTION MODULE (DATA EXISTANT) =-= FIN
- - -
E
-- NOM DU FICHIER CONTENANT LES DONNEES ?
data/ssor
++ OPEN(10,FILE='data/ssor',SPEC='OLD',RECL=0)
-- PARAMETRE D'IMPRESSION POUR L'EXECUTION ?
1
++ OPEN(11,FILE='nopo',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(12,FILE='tae',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(13,FILE='bdcl',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(14,FILE='b',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPGC :
FIN DU MODULE PREPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSAMA :
FIN DU MODULE ASSAMA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
FIN DU MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMGC :
FIN DU MODULE CLIMGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE SSORGC :
INITIALISATION DU G.C. A ZERO
NOMBRE D'ITERATIONS DU SYSTEME NUMERO 1 : 231
** TEMPS RESOLUTION 0.151587E+03 S.
FIN DU MODULE SSORGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
IMPRESSION DE LA S.D. B DE NIVEAU 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITRE :
DATE ET NOM UTILISATEUR : 01/06/92 joly
TYPE DE LA STRUCTURE DE DONNEES : B
NIVEAU ET NUMERO D'ETAT : 1 2
NOMBRE DE TABLEAUX ASSOCIES : 0
TABLEAU B 2
------------
TYPE DU TABLEAU (NTYT) : 5
NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2
TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1
NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 1
NOMBRE DE NOEUDS (NOE) : 937
NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1
NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3
NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 2811
CODE DE STOCKAGE DE B (NCODSB) : -1
TABLEAU B 3
------------
POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B
1 0 2 2811
Enfin le problème modèle est résolu par itérations de gradient
conjugué avec préconditionnement par factorisation incomplète de Cholesky
par niveau (module ICHRGC). On utilise une factorisation
incomplète de niveau 0, et la matrice de préconditionnement est compressée
à l'aide du module CCAMAT, avec un coefficient de tri égal à 
.
nopo $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD NOPO
1 5 3 $ NDSM NTYP ND
tae $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD TAE
1 $ 1 SI BDCL EXISTE
bdcl $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD BDCL
0 $ 1 SI CL. EN RL. EXISTE
b $ NOM DU FICHIER
1 $ ET NIVEAU DE LA SD B
3 $ IMPRESSION
2 $ METHODE
1 $ FACTORISATION
0 $ NIVEAU DE FACTORISATION
1.000000 $ FACTEUR DE STABILISATION
1 0.1000000E-04 $ 1 SI COMPRESSION RAPORT
0.1000000E-03 $ EPS
Les résultats correspondant à cette dernière méthode :
M M OOO DDDD U U L EEEEE FFFFF
MM MM O O D D U U L E F
M M M O O D D U U L EEEE FFFF
M M O O D D U U L E F
M M OOO DDDD UUU LLLLL EEEEE F VERSION 92
DATE : 01/06/92
AUTEUR : joly
********************************************
APPEL DE LA RESOLUTION PAR GRADIENT CONJUGUE
********************************************
-- CREATION DU DATA =-= EXECUTION MODULE (DATA EXISTANT) =-= FIN
- - -
E
-- NOM DU FICHIER CONTENANT LES DONNEES ?
data/ichr
++ OPEN(10,FILE='data/ichr',SPEC='OLD',RECL=0)
-- PARAMETRE D'IMPRESSION POUR L'EXECUTION ?
1
++ OPEN(11,FILE='nopo',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(12,FILE='tae',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(13,FILE='bdcl',SPEC='OLD,UNFORMATTED',RECL=0)
++ OPEN(14,FILE='b',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPGC :
FIN DU MODULE PREPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSAMA :
FIN DU MODULE ASSAMA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
FIN DU MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMGC :
FIN DU MODULE CLIMGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE FANIGC :
** NOMBRE D'INCONNUES 2811
** NIVEAU DE FACTORISATION 0
** LONGUEUR DU TABLEAU AMAT5 65058
** LONGUEUR DU TABLEAU AMAT6 130116
** LONGUEUR DU TABLEAU AMAC5 65058
** LONGUEUR DU TABLEAU AMAC6 130116
** TEMPS CALCUL FACTORISATION 0.567541E+01 S.
FIN DU MODULE FANIGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CCAMAT :
NOMBRE DE COEFFICIENTS ELIMINES : 9724
FIN DU MODULE CCAMAT
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ICHRGC :
FACTORISATION : 1 CHOLESKY, 2 CROUT (NFACTO) : 1
INITIALISATION DU GRADIENT CONJUGUE A ZERO
** TEMPS CALCUL RESOLUTION 0.522310E+02 S.
** NOMBRE D'ITERATIONS 80
FIN DU MODULE ICHRGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
IMPRESSION DE LA S.D. B DE NIVEAU 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITRE :
DATE ET NOM UTILISATEUR : 01/06/92 joly
TYPE DE LA STRUCTURE DE DONNEES : B
NIVEAU ET NUMERO D'ETAT : 1 2
NOMBRE DE TABLEAUX ASSOCIES : 0
TABLEAU B 2
------------
TYPE DU TABLEAU (NTYT) : 5
NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2
TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1
NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 1
NOMBRE DE NOEUDS (NOE) : 937
NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1
NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3
NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 2811
CODE DE STOCKAGE DE B (NCODSB) : -1
TABLEAU B 3
------------
POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B
1 0 2 2811