1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
|
# CONFIGURATION MANAGEMENT OF EDF VERSION
# ======================================================================
# COPYRIGHT (C) 1991 - 2012 EDF R&D WWW.CODE-ASTER.ORG
# THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
# IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
# THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
# (AT YOUR OPTION) ANY LATER VERSION.
#
# THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
# WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
# GENERAL PUBLIC LICENSE FOR MORE DETAILS.
#
# YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
# ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
# 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
# ======================================================================
# person_in_charge: jacques.pellet at edf.fr
def factoriser_prod(MATR_ASSE,**args):
if AsType(MATR_ASSE) == matr_asse_depl_r : return matr_asse_depl_r
if AsType(MATR_ASSE) == matr_asse_depl_c : return matr_asse_depl_c
if AsType(MATR_ASSE) == matr_asse_temp_r : return matr_asse_temp_r
if AsType(MATR_ASSE) == matr_asse_temp_c : return matr_asse_temp_c
if AsType(MATR_ASSE) == matr_asse_pres_r : return matr_asse_pres_r
if AsType(MATR_ASSE) == matr_asse_pres_c : return matr_asse_pres_c
raise AsException("type de concept resultat non prevu")
FACTORISER=OPER(nom="FACTORISER",op=14,sd_prod=factoriser_prod,
fr="Factoriser une matrice assemblée en un produit de deux matrices triangulaires"+
"ou construire une matrice de préconditionnement pour une résolution par gradient conjugué",
reentrant='f', UIinfo={"groupes":("Résolution",)},
regles=(EXCLUS('BLOC_DEBUT','DDL_DEBUT'),
EXCLUS('BLOC_FIN','DDL_FIN'),),
MATR_ASSE =SIMP(statut='o',typ=(matr_asse_depl_r,matr_asse_depl_c,matr_asse_temp_r,
matr_asse_temp_c,matr_asse_pres_r,matr_asse_pres_c) ),
# mots clés pour solveur LDLT et MULT_FRONT et MUMPS:
STOP_SINGULIER =SIMP(statut='f',typ='TXM',defaut="OUI",into=("OUI","NON") ),
NPREC =SIMP(statut='f',typ='I',defaut=8),
# mots clés pour solveur LDLT :
BLOC_DEBUT =SIMP(statut='f',typ='I',val_min=1,),
DDL_DEBUT =SIMP(statut='f',typ='I',val_min=1,),
BLOC_FIN =SIMP(statut='f',typ='I',val_min=1,),
DDL_FIN =SIMP(statut='f',typ='I',val_min=1,),
# mots clés pour solveur MUMPS :
TYPE_RESOL =SIMP(statut='f',typ='TXM',defaut='AUTO',into=('NONSYM','SYMGEN','SYMDEF','AUTO')),
PRETRAITEMENTS =SIMP(statut='f',typ='TXM',defaut="AUTO",into=("SANS","AUTO")),
PCENT_PIVOT =SIMP(statut='f',typ='I',defaut=20,),
ELIM_LAGR2 =SIMP(statut='f',typ='TXM',defaut="OUI",into=("OUI","NON")),
GESTION_MEMOIRE =SIMP(statut='f',typ='TXM',defaut="IN_CORE",into=("IN_CORE","OUT_OF_CORE","EVAL")),
# mots clés pour solveur GCPC et PETSc :
PRE_COND =SIMP(statut='f',typ='TXM',into=("LDLT_INC","JACOBI","SOR","LDLT_SP"),defaut="LDLT_INC" ),
b_ldlt_inc =BLOC(condition = "PRE_COND == 'LDLT_INC' ",
NIVE_REMPLISSAGE = SIMP(statut='f',typ='I',defaut= 0 ),
REMPLISSAGE = SIMP(statut='f',typ='R',defaut= 1.0),
),
b_ldlt_sp =BLOC(condition = "PRE_COND == 'LDLT_SP' ", fr="Paramètres de la factorisation simple précision",
REAC_PRECOND =SIMP(statut='f',typ='I',defaut=30, ),
),
b_jacobi_sor =BLOC(condition = "PRE_COND == 'JACOBI' or PRE_COND == 'SOR'",
fr="Paramètres des préconditionneurs JACOBI et SOR",
RENUM =SIMP(statut='f',typ='TXM',defaut="SANS",into=("SANS","RCMK") ),
),
TITRE =SIMP(statut='f',typ='TXM',max='**'),
INFO =SIMP(statut='f',typ='I',into=(1,2) ),
) ;
|
