![[BIG]](../icons/zoom18.gif){kind=icon}
![[Normal]](../icons/zoom14.gif){kind=icon}
![[small]](../icons/zoom10.gif){kind=icon}
Next: 1.3 The non-symmetric model problem Up: 1 Examples Prev: 1.1 Introduction Contents
Next: 1.3 The non-symmetric model problem
Up: 1 Examples
Prev: 1.1 Introduction
Contents
The example chosen is the solution of a realistic linear elasticity problem by the finite element method: a hook, greased on part of its boundary, is submitted to a vertical traction force on the utility part. Given the envisaged boundary conditions, and the material assumed homogeneous isotropic, the problem exhibits a symmetry plane, so that the domain is reduced to a half-hook (see Figure 1).
For more details about the variational formulation of this problem, as well as its finite element approximation and the different aspects of the computer implementation, consult also [1]. Other numerical experiences concerning the solution of this problem are also found in this reference.
During numerical implementation, several distinctive steps appear before the actual solution:
(data structure NOPO)
by module APNOPO,
Consult the corresponding guides for the utilization of these modules.
Figure 1.1: Mesh of the domain
The geometry of the model problem is described by a NOPO data structure, and a view of the object is presented in Figure 1.1. The mesh contains 1006 elements and 937 nodes, the finite elements used being of degree 1.
At each of the nodes, we need to compute the three displacement components, so that the resulting linear system has 2811 unknowns. The structure of the matrix of the linear system is represented in figure 1.2.
To solve the linear system corresponding to the model problem, we need to:
We have already mentioned that each of these steps are performed by a distinct module; to simplify the utilization of the MODULEF code, sequences of these modules are programmed for each option available (the preprocessors of library PPAL). The user can also create his/her own sequences from the MODULEF programs.
We will now look at the different possibilities:
Main program CHOLXX successively calls modules PREPAC, ASSMUA, ASEMBV, CHOLPC and, lastly, DRCHPC. The listing of the data of CHOLXX is shown below:
mail $ FILE NAME
1 $ AND LEVEL OF THE D.S. MAIL
1 5 3 $ NDSM NTYP ND
tae $ FILE NAME
1 $ AND LEVEL OF THE D.S. TAE
1 $ 1 IF BDCL IS USED , 0 IF NOT
bdcl $ FILE NAME
1 $ AND LEVEL OF THE D.S. BDCL
0 $ 1 IF BOUNDARY CONDITIONS I.T.O.
$ LINEAR RELATIONS EXIST
B $ FILE NAME
1 $ AND LEVEL OF THE D.S. B
3 $ IMPREB
The results corresponding to an execution on a Hewlett-Packard Apollo series 400 workstation is given below:
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 93
DATE : 14/06/93
AUTEUR : dutoit
**************************************
EXECUTING THE CHOLESKY SOLUTION METHOD
**************************************
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
E
-- NAME OF DATA FILE ?
dchol.data_gb
-- PRINT PARAMETER FOR THE EXECUTION ?
3
++ 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='crochet.bchol',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAC :
NUMBER OF D.O.F. PER NODE (ND) : 3
NUMBER OF LOAD CASES (NDSM) : 1
PRESENCE OF LINEAR RELATIONS (NCLRL) : 0
MATRIX STORAGE CODE (NCODSA) : 1
TYPE OF TABLES (NTY) : 5
INPUT D.S. MAIL (NFMAIL ET NIMAIL) : 11 1
INPUT D.S. NDL1 (NFNDL1 ET NINDL1) : 0 0
INPUT D.S. BDCL (NFBDCL ET NIBDCL) : 13 1
OUTPUT D.S. MUA (NFMUA ET NIMUA ) : 0 0
NUMBER OF ASSOCIATED TABLES (NTMUA) : 0
NUMBER OF WORDS OF THE MATRIX (LMUA5) : 737202
MAX DIFFERENCE BETWEEN 2 NODES OF AN ELEMENT (LBDP) : 70
MAX DIFFERENCE BETWEEN 2 D.O.F. OF AN ELEMENT (LBDPDL) : 210
MAX NUMBER OF NODES PER ELEMENT (NNOMAX) : 8
NUMBER OF RHS VECTORS (NDSM) : 1
NUMBER OF WORDS IN C.M. FOR D.S. NDL1 (MCNDL1) : 0
NUMBER OF WORDS IN C.M. FOR D.S. MUA (MCMUA) : 740060
NUMBER OF WORDS IN C.M. FOR D.S. B (MCB) : 5667
END OF MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSMUA :
ARRAY OF TAE TO ASSEMBLE (NOT) : 1
INPUT D.S. TAE (NFTAE AND NITAE) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE AND NIMUAE) : 0 0
OUTPUT D.S. MUA (NFMUAS AND NIMUAS) : 0 0
ASSEMBLY OF ELEMENT ARRAY 1 HAS 1 INDEX(CES), WITH STORAGE CODE 1
STORAGE CODE OF THE GLOBAL MATRIX 1
END OF MODULE ASSMUA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
ARRAY TO ASSEMBLE (NOT) : 2
NUMBER OF LOAD CASES (NDSM) : 1
INPUT D.S. TAE (NFTAE AND NITAE) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. B (NFB AND NIB) : 0 1
NUMBER OF ASSOCIATED TABLES (NTB) : 0
NUMBER OF WORDS IN CENTRAL MEMORY FOR THE D.S. B : 5667
END OF MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PROCESSING OPTION (NIVO) : 3
PENALTY VALUE (VTG) : 1.000000
INPUT D.S. NDL1 (NFNDL1 ET NINDL1 ) : 0 0
INPUT D.S. MUA (NFMUAE ET NIMUAE ) : 0 0
INPUT D.S. B (NFBE ET NIBE ) : 0 1
INPUT D.S. BDCL (NFBDCL ET NIBDCL ) : 13 1
OUTPUT D.S. MUA (NFMUAS ET NIMUAS ) : 0 0
OUTPUT D.S. B (NFBS ET NIBS ) : 0 1
END OF MODULE CLIMPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CHOLPC :
PRECISION (EPS) : 0.1000000E-14
RETURN OPTION (NENTRE) : 0
INPUT D.S. NDL1 (NFNDL1 ET NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE ET NIMUAE) : 0 0
OUTPUT D.S. MUA (NFMUAS ET NIMUAS) : 0 0
** FACTORISATION EXECUTION TIME 0.504265E+02 S **
END OF MODULE CHOLPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE DRCHPC :
PROCESSING OPTION (NIVEAU) : 3
PRESENCE OF BOUNDARY CONDITIONS IN L.R. (NCLRL) : 0
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE AND NIMUAE) : 0 0
INPUT D.S. B (NFBE AND NIBE) : 0 1
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
OUTPUT D.S. B (NFBS AND NIBS) : 14 1
** RESOLUTION EXECUTION TIME 0.382422E+01 S **
END OF MODULE DRCHPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PRINT OUT OF D.S. B OF LEVEL 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITLE :
DATE AND USER'S NAME : 06/05/93 pate
DATA STRUCTURE TYPE : B
LEVEL AND STATE NUMBER : 1 2
NUMBER OF ASSOCIATED TABLES : 0
TABLE B 2
----------
TYPE OF THE TABLE (NTYT) : 5
NUMBER OF INDICES AND THEIR MAX VALUES (NIND..) : 2 1 2811
OPTION (1: PAGES OF SAME SIZE , 0: IF NOT) (NCOD) : 1
NUMBER OF PAGES OF THE TABLE B4 (NBLOC) : 1
NUMBER OF NODES (NOE) : 937
NUMBER OF TABLES B4 IN THIS D.S. (NBBLOC) : 1
CONSTANT NUMBER OF D.O.F. AT EACH NODE OR 0 (ND) : 3
NUMBER OF D.O.F. OR LENGHT OF A B4 PAGE (NTDL) : 2811
STORAGE OPTION FOR B4 (NCODSB) : -1
TABLE B 3
----------
POINTER ON THE LAST ROW OF EACH PAGE OF B4
1 0 2 2811
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
F
Main program CROUXX successively calls modules PREPAC, ASSMUA, ASEMBV, CROUPC and, lastly, DRCRPC. The listing of the data of CROUXX is shown below:
crochet.mail $ FILE NAME
1 $ AND LEVEL OF THE D.S.
1 5 3 $ NDSM NTYP ND
crochet.tae $ FILE NAME
1 $ AND LEVEL OF THE D.S.
1 $ 1 IF BDCL EIS USED , 0 IF NOT
crochet.bdcl $ FILE NAME
1 $ AND LEVEL OF THE D.S. BDCL
0 $ 1 IF B.C. EXIST AS L.R
crochet.b $ FILE NAME
1 $ AND LEVEL OF THE D.S.
3 $ IMPREB
The corresponding results are given below:
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 93
DATE : 14/06/93
AUTEUR : dutoit
************************************************************
EXECUTING THE CROUT SOLUTION METHOD WITH SEQUENTIAL FILES
************************************************************
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
E
-- NAME OF DATA FILE ?
dcrou.data_gb
-- PRINT PARAMETER FOR THE EXECUTION ?
3
++ 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='crochet.bcrou',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAC :
NUMBER OF D.O.F. PER NODE (ND) : 3
NUMBER OF LOAD CASES (NDSM) : 1
PRESENCE OF LINEAR RELATIONS (NCLRL) : 0
MATRIX STORAGE CODE (NCODSA) : 1
TYPE OF TABLES (NTY) : 5
INPUT D.S. MAIL (NFMAIL ET NIMAIL) : 11 1
INPUT D.S. NDL1 (NFNDL1 ET NINDL1) : 0 0
INPUT D.S. BDCL (NFBDCL ET NIBDCL) : 13 1
OUTPUT D.S. MUA (NFMUA ET NIMUA ) : 0 0
NUMBER OF ASSOCIATED TABLES (NTMUA) : 0
NUMBER OF WORDS OF THE MATRIX (LMUA5) : 737202
MAX DIFFERENCE BETWEEN 2 NODES OF AN ELEMENT (LBDP) : 70
MAX DIFFERENCE BETWEEN 2 D.O.F. OF AN ELEMENT (LBDPDL) : 210
MAX NUMBER OF NODES PER ELEMENT (NNOMAX) : 8
NUMBER OF RHS VECTORS (NDSM) : 1
NUMBER OF WORDS IN C.M. FOR D.S. NDL1 (MCNDL1) : 0
NUMBER OF WORDS IN C.M. FOR D.S. MUA (MCMUA) : 740060
NUMBER OF WORDS IN C.M. FOR D.S. B (MCB) : 5667
END OF MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSMUA :
ARRAY OF TAE TO ASSEMBLE (NOT) : 1
INPUT D.S. TAE (NFTAE AND NITAE) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE AND NIMUAE) : 0 0
OUTPUT D.S. MUA (NFMUAS AND NIMUAS) : 0 0
ASSEMBLY OF ELEMENT ARRAY 1 HAS 1 INDEX(CES), WITH STORAGE CODE 1
STORAGE CODE OF THE GLOBAL MATRIX 1
END OF MODULE ASSMUA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
ARRAY TO ASSEMBLE (NOT) : 2
NUMBER OF LOAD CASES (NDSM) : 1
INPUT D.S. TAE (NFTAE AND NITAE) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. B (NFB AND NIB) : 0 1
NUMBER OF ASSOCIATED TABLES (NTB) : 0
NUMBER OF WORDS IN CENTRAL MEMORY FOR THE D.S. B : 5667
END OF MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PROCESSING OPTION (NIVO) : 3
PENALTY VALUE (VTG) : 1.000000
INPUT D.S. NDL1 (NFNDL1 ET NINDL1 ) : 0 0
INPUT D.S. MUA (NFMUAE ET NIMUAE ) : 0 0
INPUT D.S. B (NFBE ET NIBE ) : 0 1
INPUT D.S. BDCL (NFBDCL ET NIBDCL ) : 13 1
OUTPUT D.S. MUA (NFMUAS ET NIMUAS ) : 0 0
OUTPUT D.S. B (NFBS ET NIBS ) : 0 1
END OF MODULE CLIMPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CROUPC :
PRECISION OF MEMORY WORD (EPS) : 0.1000000E-14
RETURN OPTION (NENTRE) : 0
IDS NDL1 (NFNDL1 AND NINDL1) : 0 0
IDS MUA (NFMUAE AND NIMUAE) : 0 0
IDS MUA (NFMUAS AND NIMUAS) : 0 0
** FACTORISATION EXECUTION TIME 0.884851E+02 S **
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE CROUPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE DRCRPC :
PROCESSING OPTION (NIVEAU) : 3
PRESENCE OF BOUNDARY CONDITIONS IN L.R. (NCLRL) : 0
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE AND NIMUAE) : 0 0
INPUT D.S. B (NFBE AND NIBE) : 0 1
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
OUTPUT D.S. B (NFBS AND NIBS) : 14 1
** RESOLUTION EXECUTION TIME 0.224702E+01 S **
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE DRCRPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PRINT OUT OF D.S. B OF LEVEL 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITLE :
DATE AND USER'S NAME : 06/05/93 pate
DATA STRUCTURE TYPE : B
LEVEL AND STATE NUMBER : 1 2
NUMBER OF ASSOCIATED TABLES : 0
TABLE B 2
----------
TYPE OF THE TABLE (NTYT) : 5
NUMBER OF INDICES AND THEIR MAX VALUES (NIND..) : 2 1 2811
OPTION (1: PAGES OF SAME SIZE , 0: IF NOT) (NCOD) : 1
NUMBER OF PAGES OF THE TABLE B4 (NBLOC) : 1
NUMBER OF NODES (NOE) : 937
NUMBER OF TABLES B4 IN THIS D.S. (NBBLOC) : 1
CONSTANT NUMBER OF D.O.F. AT EACH NODE OR 0 (ND) : 3
NUMBER OF D.O.F. OR LENGHT OF A B4 PAGE (NTDL) : 2811
STORAGE OPTION FOR B4 (NCODSB) : -1
TABLE B 3
----------
POINTER ON THE LAST ROW OF EACH PAGE OF B4
1 0 2 2811
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
F
Since the memory space required to store this matrix is quite large, it is advantageous to use module CHOLPS where the matrix is stored in secondary memory. Main program CHOSXX successively calls modules PREPAC, ASMAPS, ASMBMS CHOLPS and, lastly, DRCHPS. The listing of the data of CHOSXX is given below:
crochet.mail $ FILE NAME
1 $ AND LEVEL OF THE D.S. MAIL
1 5 3 $ NDSM NTYP ND
crochet.tae $ FILE NAME
1 $ AND LEVEL OF THE D.S. TAE
1 $ 1 IF BDCL IS USED , 0 IF NOT
crochet.bdcl $ FILE NAME
1 $ AND LEVEL OF THE D.S. BDCL
0 $ 1 IF BOUNDARY CONDITIONS I.T.O.
$ LINEAR RELATIONS EXIST
crochet.mua $ FILE NAME
1 $ AND LEVEL OF THE D.S. MUA
crochet.b $ FILE NAME
1 $ AND LEVEL OF THE D.S. B
4 $ IMPREB
4096 $ LPAGE
The corresponding results are given below:
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 93
DATE : 15/06/93
AUTEUR : dutoit
**************************************
EXECUTING THE CHOLESKY SOLUTION METHOD
**************************************
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
E
-- NAME OF DATA FILE ?
dchos.data-gb
-- PRINT PARAMETER FOR THE EXECUTION ?
3
++ 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 :
NUMBER OF D.O.F. PER NODE (ND) : 3
NUMBER OF LOAD CASES (NDSM) : 1
PRESENCE OF LINEAR RELATIONS (NCLRL) : 0
MATRIX STORAGE CODE (NCODSA) : 1
TYPE OF TABLES (NTY) : 5
INPUT D.S. MAIL (NFMAIL ET NIMAIL) : 11 1
INPUT D.S. NDL1 (NFNDL1 ET NINDL1) : 0 0
INPUT D.S. BDCL (NFBDCL ET NIBDCL) : 13 1
OUTPUT D.S. MUA (NFMUA ET NIMUA ) : 0 0
NUMBER OF ASSOCIATED TABLES (NTMUA) : 0
NUMBER OF WORDS OF THE MATRIX (LMUA5) : 737202
MAX DIFFERENCE BETWEEN 2 NODES OF AN ELEMENT (LBDP) : 70
MAX DIFFERENCE BETWEEN 2 D.O.F. OF AN ELEMENT (LBDPDL) : 210
MAX NUMBER OF NODES PER ELEMENT (NNOMAX) : 8
NUMBER OF RHS VECTORS (NDSM) : 1
NUMBER OF WORDS IN C.M. FOR D.S. NDL1 (MCNDL1) : 0
NUMBER OF WORDS IN C.M. FOR D.S. MUA (MCMUA) : 740060
NUMBER OF WORDS IN C.M. FOR D.S. B (MCB) : 5667
END OF MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
lpage = 4096
nbpag = 6
++ OPEN(14,FILE='mua',SPEC='DIRECT',RECL=16384)
++ OPEN(15,FILE='crochet.bchos',SPEC='DIRECT',RECL=16384)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAC :
NUMBER OF D.O.F. PER NODE (ND) : 3
NUMBER OF LOAD CASES (NDSM) : 1
PRESENCE OF LINEAR RELATIONS (NCLRL) : 0
MATRIX STORAGE CODE (NCODSA) : 1
TYPE OF TABLES (NTY) : 5
INPUT D.S. MAIL (NFMAIL ET NIMAIL) : 11 1
INPUT D.S. NDL1 (NFNDL1 ET NINDL1) : 0 0
INPUT D.S. BDCL (NFBDCL ET NIBDCL) : 13 1
OUTPUT D.S. MUA (NFMUA ET NIMUA ) : -14 1
NUMBER OF ASSOCIATED TABLES (NTMUA) : 0
NUMBER OF WORDS OF THE MATRIX (LMUA5) : 737202
MAX DIFFERENCE BETWEEN 2 NODES OF AN ELEMENT (LBDP) : 70
MAX DIFFERENCE BETWEEN 2 D.O.F. OF AN ELEMENT (LBDPDL) : 210
MAX NUMBER OF NODES PER ELEMENT (NNOMAX) : 8
NUMBER OF RHS VECTORS (NDSM) : 1
NUMBER OF WORDS IN C.M. FOR D.S. NDL1 (MCNDL1) : 0
NUMBER OF WORDS IN C.M. FOR D.S. MUA (MCMUA) : 740247
NUMBER OF WORDS IN C.M. FOR D.S. B (MCB) : 5667
END OF MODULE PREPAC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASMAPS :
ARRAY TAE TO ASSEMBLE (NOT) : 1
I.D.S. TAE (NFTAE AND NITAE ) : 12 1
O.D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
O.D.S. MUA (NFMUAE AND NIMUAE) : -14 0
O.D.S. MUA (NFMUAS AND NIMUAS) : -14 0
ASSEMBLY OF ELEMENT ARRAY 1 WITH 1 INDEX(S),
WITH STORAGE CODE 1
STORAGE CODE OF GLOBAL MATRIX 1
SP INITMP : NUMBER OF PAGES IN MAIN MEMORY : 35
NUMBER OF WORDS IN MAIN MEMORY OF MUA 7141
END OF MODULE ASMAPS
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASMBMS :
ARRAY TAE TO ASSEMBLE (NOT) : 2
NUMBER OF LOAD CASES (NDSM) : 1
I.D.S. TAE (NFTAE AND NITAE ) : 12 1
I.D.S. NDL1 (NFNDL1 ET NINDL1) : 0 0
O.D.S. B (NFB ET NIB ) : -15 1
NUMBER OF ASSOCIATED TABLES (NTB) : 0
SP INITMP : NUMBER OF PAGES IN MAIN MEMORY : 2
ASSEMBLY OF ELEMENT ARRAY 1 WITH 1 INDEX(S),
WITH STORAGE CODE 1
STORAGE CODE OF RIGHT HAND SIDE -1
NUMBER OF WORDS IN CENTRAL MEMORY FOR THE D.S. B : 4142
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
FIN DU MODULE ASMBMS
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMPS :
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PROCESSING OPTION (NIVO) : 3
PENALTY VALUE (VTG) : 1.000000
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE AND NIMUAE) : -14 0
INPUT D.S. B (NFBE AND NIBE) : -15 1
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
INPUT D.S. MUA (NFMUAS AND NIMUAS) : -14 0
OUTPUT D.S. B (NFBS AND NIBS) : -15 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE CLIMPS
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CHOLPS :
PRECISION (EPS) : 0.1000000E-14
RETURN OPTION (NENTRE) : 0
INPUT D.S. NDL1 (NFNDL1 ET NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE ET NIMUAE) : -14 0
OUTPUT D.S. MUA (NFMUAS ET NIMUAS) : -14 0
SP INITMP : NUMBER OF PAGES IN MAIN MEMORY : 35
** FACTORISATION EXECUTION TIME 0.996348E+02 S **
END OF MODULE CHOLPC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE DRCHPS :
PROCESSING OPTION (NIVEAU) : 3
PRESENCE OF BOUNDARY CONDITIONS IN L.R. (NCLRL) : 0
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. MUA (NFMUAE AND NIMUAE) : -14 0
INPUT D.S. B (NFBE AND NIBE) : -15 1
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
OUTPUT D.S. B (NFBS AND NIBS) : -15 1
SP INITMP : NUMBER OF PAGES IN MAIN MEMORY : 2
** SOLUTION EXECUTION TIME 0.267583E+02 S **
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PRINT OUT OF D.S. B OF LEVEL 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITLE :
DATE AND USER'S NAME : 09/07/93 pate
DATA STRUCTURE TYPE : B
LEVEL AND STATE NUMBER : 1 2
NUMBER OF ASSOCIATED TABLES : 0
TABLE B 2
----------
TYPE OF THE TABLE (NTYT) : 5
NUMBER OF INDICES AND THEIR MAX VALUES (NIND..) : 2 1 2048
OPTION (1: PAGES OF SAME SIZE , 0: IF NOT) (NCOD) : 1
NUMBER OF PAGES OF THE TABLE B4 (NBLOC) : 2
NUMBER OF NODES (NOE) : 937
NUMBER OF TABLES B4 IN THIS D.S. (NBBLOC) : 1
CONSTANT NUMBER OF D.O.F. AT EACH NODE OR 0 (ND) : 3
NUMBER OF D.O.F. OR LENGHT OF A B4 PAGE (NTDL) : 2811
STORAGE OPTION FOR B4 (NCODSB) : -1
TABLE B 3
----------
POINTER ON THE LAST ROW OF EACH PAGE OF B4
1 0 2 2048 3 2811
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
F
The optimization of memory space is accomplished by the frontal method, which manages this resource to the best, by storing only the indispensable rows of the linear system.
Main program FRONXX successively calls modules PREPAF and FRONT. The listing of the data of FRONXX is shown below:
mail $ FILE NAME
1 $ AND LEVEL OF THE D.S.
1 5 3 $ NDSM NTYP ND
tae $ FILE NAME
1 $ AND LEVEL OF THE D.S.
1 $ 1 IF BOUNDARY CONDITIONS I.T.O.
$ LINEAR RELATIONS EXIST
bdcl $ FILE NAME
1 $ AND LEVEL OF THE D.S.
crochet.bfron $ FILE NAME
1 $ AND LEVEL OF THE D.S.
3 $ IMPREB
2 $ INCREASE OF THE NUMBER OF FRONTS
The corresponding results are given below. Files fron and gaus are temporary files used by modules PREPAF and 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 93
DATE : 15/06/93
AUTEUR : dutoit
*****************************************************************
EXECUTION OF THE FRONTAL SOLUTION METHOD WITH SEQUENTIAL FILES
*****************************************************************
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
E
-- NAME OF DATA FILE ?
dfron.data_gb
-- PRINT PARAMETER FOR THE EXECUTION ?
3
++ 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='crochet.bfron',SPEC='UNFORMATTED',RECL=0)
++ OPEN(15,FILE='',SPEC='SCRATCH,UNFORMATTED',RECL=0)
++ OPEN(16,FILE='',SPEC='SCRATCH,UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPAF :
MAXIMUM NUMBER OF NODES IN AN ELEMENT (NNOMAX) : 8
CONSTANT NUMBER OF D.O.F. PER NODES (OR 0) (ND) : 3
MAXIMUM NUMBER OF D.O.F. AT A NODE (NDLMAX) : 3
TOTAL NUMBER OF FRONTS (NTYPFR) : 1
NUMBER OF D.O.F. PER FRONT :
FRONT NUMBER 1 : 3
NUMBER OF ELEMENTS : 1006
NUMBER OF NODES : 937
NUMBER OF FRONTS : 1
LENGTH OF FRONT EXPECTED ( FOR ANY TYPE OF FRONT ) : 937
FRONT NUMBER : 1 (NODES SUPPORTING 3 D.O.F.)
LENGTH OF FRONT : 242
NUMBER OF DEGREES OF FREEDOM OF SYSTEM : 2811
NUMBER OF DEGREES OF FREEDOM OF FRONTAL SYSTEM : 726
SIZE OF NON-SYMMETRIC FRONTAL MATRIX : 527076
SIZE OF SYMMETRIC FRONTAL MATRIX : 263901
SIZE OF FRONTAL R.H.S. : 726
** PREPARATION EXECUTION TIME 0.166020E+01 S **
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE PREPAF
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE FRONT :
SOLUTION OF A SYMMETRIC SYSTEM
TOTAL NO OF ELEMENTSB 1006
TOTAL NUMBER OF NODES 937
RANK OF SYSTEM 726
NUMBER OF R.H.S. 1
** RESOLUTION EXECUTION TIME 0.780065E+03 S **
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE FRONT
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PRINT OUT OF D.S. B OF LEVEL 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITLE :
DATE AND USER'S NAME : 06/05/93 pate
DATA STRUCTURE TYPE : B
LEVEL AND STATE NUMBER : 1 0
NUMBER OF ASSOCIATED TABLES : 0
TABLE B 2
----------
TYPE OF THE TABLE (NTYT) : 5
NUMBER OF INDICES AND THEIR MAX VALUES (NIND..) : 2 1 2811
OPTION (1: PAGES OF SAME SIZE , 0: IF NOT) (NCOD) : 1
NUMBER OF PAGES OF THE TABLE B4 (NBLOC) : 1
NUMBER OF NODES (NOE) : 937
NUMBER OF TABLES B4 IN THIS D.S. (NBBLOC) : 1
CONSTANT NUMBER OF D.O.F. AT EACH NODE OR 0 (ND) : 3
NUMBER OF D.O.F. OR LENGHT OF A B4 PAGE (NTDL) : 2811
STORAGE OPTION FOR B4 (NCODSB) : -1
TABLE B 3
----------
POINTER ON THE LAST ROW OF EACH PAGE OF B4
1 0 2 2811
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
F
If we want to solve the model problem by the conjugate gradient method, we need to assemble the matrix in a compact form, an AMAT data structure. Module GRADXX successively calls modules PREPGC, ASSAMA, ASEMBV, then depending on the type of preconditioning chosen, one of the following modules:
The model problem is firstly solved by conjugate gradient iterations without preconditioning (module SIMPGC), using the following data:
crochet.nopo $ FILE NAME
1 $ AND LEVEL OF THE D.S. NOPO
1 5 3 $ NDSM NTYP ND
crochet.tae $ FILE NAME
1 $ AND LEVEL OF THE D.S. TAE
1 $ 1 IF BDCL EXISTS
crochet.bdcl $ FILE NAME
1 $ AND LEVEL OF THE D.S. BDCL
0 $ 1 IF BOUNDARY CONDITIONS
$ I.T.O. LINEAR RELATIONS EXIST
crochet.b $ FILE NAME
1 $ AND LEVEL OF THE D.S. B
3 $ PRINT PARAMETER
0 $ METHOD
0.1000000E-03 $ EPS
The corresponding results are given below:
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 93
DATE : 18/06/93
AUTEUR : dutoit
**************************************************************
EXECUTING OF CONJUGATE GRADIENT SOLUTION WITH SEQUENTIAL FILES
**************************************************************
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
E
-- NAME OF DATA FILE ?
dgrad.data_gb
--PRINT PARAMETER FOR EXECUTION ?
3
++ 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='crochet.bgrad',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPGC :
NUMBER OF D.O.F. PER NODE (ND) : 3
MATRIX STRORAGE CODE (NCODSA) : 1
ARRAY TYPE (NTY) : 5
PRESENCE OF LINEAR RELATIONS (NCLRL) : 0
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. NOPO (NFNOPO AND NINOPO) : 11 1
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
NUMBER OF ASSOCIATED TABLES (NTAMAT) : 0
NUMBER OF WORDS IN C.M. OF THE ASSEMBLED D.S. AMAT : 198030
END OF MODULE PREPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSAMA :
ARRAY TAE TO ASSEMBLE (NOT) : 1
INPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
INPUT D.S. TAE (NFTAE AND NITAE ) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
ASSEMBLY OF ELEMENT ARRAY 1 HAS 1 INDEX(ES), OF STORAGE CODE 1
STORAGE CODE FOR THE GLOBAL MATRIX 1
END OF MODULE ASSAMA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
ARRAY TO ASSEMBLE (NOT) : 2
NUMBER OF LOAD CASES (NDSM) : 1
INPUT D.S. TAE (NFTAE AND NITAE) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. B (NFB AND NIB) : 0 1
NUMBER OF ASSOCIATED TABLES (NTB) : 0
NUMBER OF WORDS IN CENTRAL MEMORY FOR THE D.S. B : 5667
END OF MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMGC :
PROCESSING OPTION (NIVO) : 3
PENALTY VALUE (VTG) : 1.000000
INPUT D.S. NDL1 (NFNDL1 AND NINDL1 ) : 0 0
INPUT D.S. AMAT (NFAMAT AND NIAMAT ) : 0 0
INPUT D.S. B (NFBE AND NIBE ) : 0 1
INPUT D.S. BDCL (NFBDCL AND NIBDCL ) : 13 1
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT ) : 0 0
OUTPUT D.S. B (NFBS AND NIBS ) : 0 1
END OF MODULE CLIMGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE SIMPGC :
STOP PARAMETER (EPS) : 0.100000E-03
NAME OF INITIALISATION ARRAY : NULL
NUMBER OF LINEAR COMBINATIONS (NCLRL) : 0
I.D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
I.D.S. B (NFBE AND NIBE ) : 0 1
I.D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
O.D.S. B (NFBS AND NIBS ) : 14 1
INITIALISATION OF C.G. TO ZERO
SYSTEM NUMBER : 1
--------------
NORM OF B : 0.229213E+09 MAXIMUM NUMBER OF ITERATIONS : 2811
STOP TEST (EPS1) : 2.29213
NUMBER OF ITERATIONS OF SYSTEM NUMBER 1 : 1211
** RESOLUTION EXECUTION TIME 0.432811E+03 S **
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE SIMPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PRINT OUT OF D.S. B OF LEVEL 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITLE :
DATE AND USER'S NAME : 06/05/93 pate
DATA STRUCTURE TYPE : B
LEVEL AND STATE NUMBER : 1 2
NUMBER OF ASSOCIATED TABLES : 0
TABLE B 2
----------
TYPE OF THE TABLE (NTYT) : 5
NUMBER OF INDICES AND THEIR MAX VALUES (NIND..) : 2 1 2811
OPTION (1: PAGES OF SAME SIZE , 0: IF NOT) (NCOD) : 1
NUMBER OF PAGES OF THE TABLE B4 (NBLOC) : 1
NUMBER OF NODES (NOE) : 937
NUMBER OF TABLES B4 IN THIS D.S. (NBBLOC) : 1
CONSTANT NUMBER OF D.O.F. AT EACH NODE OR 0 (ND) : 3
NUMBER OF D.O.F. OR LENGHT OF A B4 PAGE (NTDL) : 2811
STORAGE OPTION FOR B4 (NCODSB) : -1
TABLE B 3
----------
POINTER ON THE LAST ROW OF EACH PAGE OF B4
1 0 2 2811
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
f
The model problem is then solved by conjugate gradient iterations with SSOR preconditioning (module SSORGC), using the following data:
crochet.nopo $ FILE NAME
1 $ AND LEVEL OF THE D.S. NOPO
1 5 3 $ NDSM NTYP ND
crochet.tae $ FILE NAME
1 $ AND LEVEL OF THE D.S. TAE
1 $ 1 IF BDCL EXISTS
crochet.bdcl $ FILE NAME
1 $ AND LEVEL OF THE D.S. BDCL
0 $ 1 IF BOUNDARY CONDITIONS
$ I.T.O. LINEAR RELATIONS EXIST
crochet.b $ FILE NAME
1 $ AND LEVEL OF THE D.S. B
3 $ PRINT PARAMETER
1 $ METHOD
0.1000000E-03 1.000000 $ EPS OMEGA
The corresponding results are given below:
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 93
DATE : 18/06/93
AUTEUR : dutoit
**************************************************************
EXECUTING OF CONJUGATE GRADIENT SOLUTION WITH SEQUENTIAL FILES
**************************************************************
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
E
-- NAME OF DATA FILE ?
dgrad1.data_gb
--PRINT PARAMETER FOR EXECUTION ?
3
++ 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='crochet.bsor',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPGC :
NUMBER OF D.O.F. PER NODE (ND) : 3
MATRIX STRORAGE CODE (NCODSA) : 1
ARRAY TYPE (NTY) : 5
PRESENCE OF LINEAR RELATIONS (NCLRL) : 0
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. NOPO (NFNOPO AND NINOPO) : 11 1
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
NUMBER OF ASSOCIATED TABLES (NTAMAT) : 0
NUMBER OF WORDS IN C.M. OF THE ASSEMBLED D.S. AMAT : 198030
END OF MODULE PREPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSAMA :
ARRAY TAE TO ASSEMBLE (NOT) : 1
INPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
INPUT D.S. TAE (NFTAE AND NITAE ) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
ASSEMBLY OF ELEMENT ARRAY 1 HAS 1 INDEX(ES), OF STORAGE CODE 1
STORAGE CODE FOR THE GLOBAL MATRIX 1
END OF MODULE ASSAMA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
ARRAY TO ASSEMBLE (NOT) : 2
NUMBER OF LOAD CASES (NDSM) : 1
INPUT D.S. TAE (NFTAE AND NITAE) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. B (NFB AND NIB) : 0 1
NUMBER OF ASSOCIATED TABLES (NTB) : 0
NUMBER OF WORDS IN CENTRAL MEMORY FOR THE D.S. B : 5667
END OF MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMGC :
PROCESSING OPTION (NIVO) : 3
PENALTY VALUE (VTG) : 1.000000
INPUT D.S. NDL1 (NFNDL1 AND NINDL1 ) : 0 0
INPUT D.S. AMAT (NFAMAT AND NIAMAT ) : 0 0
INPUT D.S. B (NFBE AND NIBE ) : 0 1
INPUT D.S. BDCL (NFBDCL AND NIBDCL ) : 13 1
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT ) : 0 0
OUTPUT D.S. B (NFBS AND NIBS ) : 0 1
END OF MODULE CLIMGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE SSORGC :
RELAXATION PARAMETER (OMEGA) : 1.00000
STOP PARAMETER (EPS) : 0.100000E-03
NAME OF INITIALISATION ARRAY : NULL
NUMBER OF LINEAR COMBINATIONS (NCLRL) : 0
I.D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
I.D.S. B (NFBE AND NIBE ) : 0 1
I.D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
O.D.S. B (NFBS AND NIBS ) : 14 1
INITIALISATION OF C.G. TO ZERO
SYSTEM NUMBER : 1
--------------
OMEGA : 1.00000 NORM OF B : 0.229213E+09 MAXIMUM NUMBER OF ITERATIONS : 1405
STOP TEST (EPS1) : 2.29213
NUMBER OF ITERATIONS OF SYSTEM NUMBER 1 : 231
** RESOLUTION EXECUTION TIME 0.163186E+03 S **
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE SSORGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PRINT OUT OF D.S. B OF LEVEL 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITLE :
DATE AND USER'S NAME : 06/05/93 pate
DATA STRUCTURE TYPE : B
LEVEL AND STATE NUMBER : 1 2
NUMBER OF ASSOCIATED TABLES : 0
TABLE B 2
----------
TYPE OF THE TABLE (NTYT) : 5
NUMBER OF INDICES AND THEIR MAX VALUES (NIND..) : 2 1 2811
OPTION (1: PAGES OF SAME SIZE , 0: IF NOT) (NCOD) : 1
NUMBER OF PAGES OF THE TABLE B4 (NBLOC) : 1
NUMBER OF NODES (NOE) : 937
NUMBER OF TABLES B4 IN THIS D.S. (NBBLOC) : 1
CONSTANT NUMBER OF D.O.F. AT EACH NODE OR 0 (ND) : 3
NUMBER OF D.O.F. OR LENGHT OF A B4 PAGE (NTDL) : 2811
STORAGE OPTION FOR B4 (NCODSB) : -1
TABLE B 3
----------
POINTER ON THE LAST ROW OF EACH PAGE OF B4
1 0 2 2811
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
f
Lastly, the model problem is solved by conjugate gradient iterations with
preconditioning by incomplete Cholesky
factorisation by level (module ICHRGC). We use an incomplete factorisation
of level 0, and the preconditioning matrix is compressed using
module CCAMAT, with a sort coefficient equal to 
.
crochet.nopo $ FILE NAME
1 $ AND LEVEL OF THE D.S. NOPO
1 5 3 $ NDSM NTYP ND
crochet.tae $ FILE NAME
1 $ AND LEVEL OF THE D.S. TAE
1 $ 1 IF BDCL EXISTS
crochet.bdcl $ FILE NAME
1 $ AND LEVEL OF THE D.S. BDCL
0 $ 1 IF BOUNDARY CONDITIONS
$ I.T.O. LINEAR RELATIONS EXIST
crochet.b $ FILE NAME
1 $ AND LEVEL OF THE D.S. B
3 $ PRINT PARAMETER
2 $ METHOD
1 $ FACTORISATION
1 0.1000000E-04 $ 1 IF COMPRESSION RATIO
0.1000000E-03 $ EPS
The results corresponding to the last method is given below:
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 93
DATE : 18/06/93
AUTEUR : dutoit
**************************************************************
EXECUTING OF CONJUGATE GRADIENT SOLUTION WITH SEQUENTIAL FILES
**************************************************************
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
e
-- NAME OF DATA FILE ?
dgrad2.data-gb
--PRINT PARAMETER FOR EXECUTION ?
3
++ 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='crochet.bich',SPEC='UNFORMATTED',RECL=0)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE PREPGC :
NUMBER OF D.O.F. PER NODE (ND) : 3
MATRIX STRORAGE CODE (NCODSA) : 1
ARRAY TYPE (NTY) : 5
PRESENCE OF LINEAR RELATIONS (NCLRL) : 0
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
INPUT D.S. NOPO (NFNOPO AND NINOPO) : 11 1
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
NUMBER OF ASSOCIATED TABLES (NTAMAT) : 0
NUMBER OF WORDS IN C.M. OF THE ASSEMBLED D.S. AMAT : 198030
END OF MODULE PREPGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASSAMA :
ARRAY TAE TO ASSEMBLE (NOT) : 1
INPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
INPUT D.S. TAE (NFTAE AND NITAE ) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
ASSEMBLY OF ELEMENT ARRAY 1 HAS 1 INDEX(ES), OF STORAGE CODE 1
STORAGE CODE FOR THE GLOBAL MATRIX 1
END OF MODULE ASSAMA
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ASEMBV :
ARRAY TO ASSEMBLE (NOT) : 2
NUMBER OF LOAD CASES (NDSM) : 1
INPUT D.S. TAE (NFTAE AND NITAE) : 12 1
INPUT D.S. NDL1 (NFNDL1 AND NINDL1) : 0 0
OUTPUT D.S. B (NFB AND NIB) : 0 1
NUMBER OF ASSOCIATED TABLES (NTB) : 0
NUMBER OF WORDS IN CENTRAL MEMORY FOR THE D.S. B : 5667
END OF MODULE ASEMBV
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CLIMGC :
PROCESSING OPTION (NIVO) : 3
PENALTY VALUE (VTG) : 1.000000
INPUT D.S. NDL1 (NFNDL1 AND NINDL1 ) : 0 0
INPUT D.S. AMAT (NFAMAT AND NIAMAT ) : 0 0
INPUT D.S. B (NFBE AND NIBE ) : 0 1
INPUT D.S. BDCL (NFBDCL AND NIBDCL ) : 13 1
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT ) : 0 0
OUTPUT D.S. B (NFBS AND NIBS ) : 0 1
END OF MODULE CLIMGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE FANIGC :
CHOLESKY INCOMPLETE FACTORISATION : 1
INPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
OUTPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 1
LEVEL OF FILL-IN : 0
NUMBER OF UNKNOWNS 2811
NUMBER OF NONZEROS OF A 65058
INCOMPLETE FACTORISATION AT LEVEL 0
NUMBER OF NONZEROS OF CA 65058
** FACTORISATION CPU TIME 0.400442E+01 SECONDS **
END OF MODULE FANIGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE CCAMAT :
COMPRESSION RATIO (RAPORT) : 0.100000E-04
I.D.S AMAT (NFAMAT AND NIAMAT ) : 0 1
O.D.S. AMAT (NFAMAT AND NIAMAT ) : 0 1
NUMBER OF COEFFICIENTS ELIMINATED : 9724
END OF MODULE CCAMAT
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE ICHRGC :
TERMINATION PARAMETER (EPS) : 0.100000E-03
NAME OF INITIALISATION ARRAY : NULL
NUMBER OF LINEAR COMBINATIONS (NCLRL) : 0
INPUT D.S. BDCL (NFBDCL AND NIBDCL) : 13 1
INPUT D.S. B (NFBE AND NIBE ) : 0 1
INPUT D.S. AMAT (NFAMAC AND NIAMAC) : 0 1
INPUT D.S. AMAT (NFAMAT AND NIAMAT) : 0 0
OUTPUT D.S. B (NFBS AND NIBS ) : 14 1
FACTORISATION : 1 CHOLESKY, 2 CROUT (NFACTO) : 1
INITIALISATION OF CONJUGATE GRADIENT TO ZERO
SYSTEM NUMBER : 1
--------------
NORM OF B : 0.229213E+09 MAXIMUM NUMBER OF ITERATIONS : 2811
STOP TEST (EPS1) : 2.29213
NUMBER OF ITERATIONS OF SYSTEM NUMBER 1 : 80
** RESOLUTION EXECUTION TIME 0.523324E+02 S **
END OF MODULE ICHRGC
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
PRINT OUT OF D.S. B OF LEVEL 1
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
TITLE :
DATE AND USER'S NAME : 06/05/93 pate
DATA STRUCTURE TYPE : B
LEVEL AND STATE NUMBER : 1 1
NUMBER OF ASSOCIATED TABLES : 0
TABLE B 2
----------
TYPE OF THE TABLE (NTYT) : 5
NUMBER OF INDICES AND THEIR MAX VALUES (NIND..) : 2 1 2811
OPTION (1: PAGES OF SAME SIZE , 0: IF NOT) (NCOD) : 1
NUMBER OF PAGES OF THE TABLE B4 (NBLOC) : 1
NUMBER OF NODES (NOE) : 937
NUMBER OF TABLES B4 IN THIS D.S. (NBBLOC) : 1
CONSTANT NUMBER OF D.O.F. AT EACH NODE OR 0 (ND) : 3
NUMBER OF D.O.F. OR LENGHT OF A B4 PAGE (NTDL) : 2811
STORAGE OPTION FOR B4 (NCODSB) : -1
TABLE B 3
----------
POINTER ON THE LAST ROW OF EACH PAGE OF B4
1 0 2 2811
-- CREATE THE DATA =-= EXECUTE MODULE (DATA EXISTING) =-= END =-= ?
- - ---
f