1 #@ MODIF fiabilite_mefisto Macro DATE 14/09/2004 AUTEUR MCOURTOI M.COURTOIS
2 # -*- coding: iso-8859-1 -*-
3 # CONFIGURATION MANAGEMENT OF EDF VERSION
4 # ======================================================================
5 # COPYRIGHT (C) 1991 - 2004 EDF R&D WWW.CODE-ASTER.ORG
6 # THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
7 # IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
8 # THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
9 # (AT YOUR OPTION) ANY LATER VERSION.
11 # THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
12 # WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
13 # MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
14 # GENERAL PUBLIC LICENSE FOR MORE DETAILS.
16 # YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
17 # ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
18 # 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
19 # ======================================================================
22 # RESPONSABLE GNICOLAS G.NICOLAS
24 def fiabilite_mefisto ( self, Rep_Calc_LOGICIEL_global,
31 # valeurs_lois est un dictionnaire indexé sur les variables.
32 # Chaque case, valeurs_lois[m], est un dictionnaire contenant :
33 # d["v_moy_physique"] = valeur moyenne physique
34 # d["v_moy_loi"] = valeur moyenne de la loi
35 # d["v_min_loi"] = valeur minimale de la loi
36 # d["v_max_loi"] = valeur maximale de la loi
37 # d["sigma_loi"] = ecart type de la loi
39 # args est le dictionnaire des arguments optionnels
40 # args.keys() est la liste des mots-clés
41 # args.keys()[0] est la premiere valeur de cette liste
42 # args.keys()[1:] est la liste des valeurs suivantes dans cette liste
43 # args.keys(mot_cle) représente le contenu de la variable mot_cle dans la macro appelante.
45 """ Ecriture des données spécifiques à MEFISTO. """
47 from Macro import fiabilite_fichier
51 #____________________________________________________________________
54 #____________________________________________________________________
57 messages_erreur = { 0 : "Tout va bien",
58 1 : "Fichier inconnu.",
59 2 : "Problème d'ouverture de fichier.",
60 10 : "Problème d'ouverture de fichier.",
61 11 : "Problème de fermeture de fichier.",
62 20 : "Problème d'impression de fichier.",
63 50 : "Donnée inacceptable.",
66 trad_oui_non = { "OUI" : 1,
73 #____________________________________________________________________
75 # 2. Les fichiers pour le logiciel de fiabilité
76 # Ils sont créés dans le répertoire d'exécution du logiciel de fiabilité, avec leurs noms officiels
77 #____________________________________________________________________
80 fic_dataMenu = "dataMenu"
81 fic_dataStoch = "dataStoch"
82 fic_dataNum = "dataNum"
83 fic_dataGrad = "dataGrad"
85 #____________________________________________________________________
87 # 3. Construction du fichier 'dataMenu'
88 #____________________________________________________________________
90 # 3.1 ==> Ouverture du fichier
92 f_menu = fiabilite_fichier.fiabilite_fichier ( self, Rep_Calc_LOGICIEL_global, fic_dataMenu, INFO )
93 erreur = f_menu.Ouvre_Fichier ( "w" )
97 # 3.2 ==> Ecriture des données nécessaires
99 f_menu.Ecrit_Titre ("MENU DU PROGRAMME MEFISTO")
100 f_menu.Ecrit_Titre ("1 <=> OUI et 0 <=> NON (entiers)")
102 f_menu.Ecrit_Titre ("Recherche du point de conception")
103 aux = trad_oui_non[args["RECH_PT_CONCEPT"]]
104 f_menu.Ecrit_Valeurs (aux)
106 f_menu.Ecrit_Titre ("First Order Reliability Analyses")
107 aux = trad_oui_non[args["METHODE_FORM"]]
108 f_menu.Ecrit_Valeurs (aux)
110 f_menu.Ecrit_Titre ("Second Order Reliability Analyses")
111 aux = trad_oui_non[args["METHODE_SORM"]]
112 f_menu.Ecrit_Valeurs (aux)
114 f_menu.Ecrit_Titre ("Importance Sampling Analyses")
115 aux = trad_oui_non[args["TIRAGE_IMPORTANCE"]]
116 f_menu.Ecrit_Valeurs (aux)
118 f_menu.Ecrit_Titre ("Optimality Test (1) : Hessian Test")
119 aux = trad_oui_non[args["T_HESSIEN"]]
120 f_menu.Ecrit_Valeurs (aux)
122 f_menu.Ecrit_Titre ("Optimality Test (2) : Sphere Test")
123 aux = trad_oui_non[args["T_SPHERE"]]
124 f_menu.Ecrit_Valeurs (aux)
126 f_menu.Ecrit_Titre ("Optimality Test (3) : Strong Max Test")
127 aux = trad_oui_non[args["T_MAXIMUM_FORT"]]
128 f_menu.Ecrit_Valeurs (aux)
130 f_menu.Ecrit_Titre ("Experiment Plan")
131 aux = trad_oui_non[args["PLAN_EXPERIENCE"]]
132 f_menu.Ecrit_Valeurs (aux)
134 f_menu.Ecrit_Titre ("Polynomial Taylor Approximation (order 2)")
135 aux = trad_oui_non[args["POLYNOME_TAYLOR"]]
136 f_menu.Ecrit_Valeurs (aux)
138 # 3.3 ==> Fermeture du fichier
140 erreur = f_menu.Ferme_Fichier ( )
145 erreur = f_menu.Imprime_Fichier ( )
149 #____________________________________________________________________
151 # 4. Construction du fichier "dataStoch"
152 #____________________________________________________________________
154 # 4.1 ==> Ouverture du fichier
156 f_stoch = fiabilite_fichier.fiabilite_fichier ( self, Rep_Calc_LOGICIEL_global, fic_dataStoch, INFO )
157 erreur = f_stoch.Ouvre_Fichier ( "w" )
161 # 4.2 ==> Nombre d'occurence de VARIABLE
163 nb_occu_variable = len(VARIABLE)
165 # 4.3 ==> Ecriture des données nécessaires
167 f_stoch.Ecrit_Titre ("Code name")
168 aux=string.replace(VERSION,"_",".")
169 aux=string.replace(aux,"N","n")
170 aux=string.replace(aux,"V","v")
171 f_stoch.Ecrit_Valeurs ( "aster_" + aux )
173 aux = [ "Gradients evaluated by the code" ]
174 aux.append("1 : Au moins 1 ; 0 : aucun")
175 f_stoch.Ecrit_Titre (aux)
178 if m["GRADIENT"] == "OUI" : gradient = 1
179 f_stoch.Ecrit_Valeurs (gradient)
181 f_stoch.Ecrit_Titre ("Variates number")
182 f_stoch.Ecrit_Valeurs ( nb_occu_variable )
184 aux = [ "Stochastic Variates" ]
185 aux.append("1: Uniforme (min, max)")
186 aux.append("2: Normal (mean, std dev)")
187 aux.append("3: LogNormal (mean, std dev, min)")
188 aux.append("4: Normal Truncated (mean, std dev, min, max)")
189 f_stoch.Ecrit_Titre (aux)
194 if m["LOI"] == "UNIFORME" :
195 f_stoch.Ecrit_Valeurs ( [ m["NOM"], 1 , d["v_min_loi"] , d["v_max_loi"] ] )
196 elif m["LOI"] == "NORMALE" :
197 f_stoch.Ecrit_Valeurs ( [ m["NOM"], 2 , d["v_moy_loi"] , d["sigma_loi"] ] )
198 elif m["LOI"] == "LOGNORMALE" :
199 f_stoch.Ecrit_Valeurs ( [ m["NOM"], 3 , d["v_moy_loi"] , d["sigma_loi"] , d["v_min_loi"] ] )
200 elif m["LOI"] == "NORMALE_TRONQUEE" :
201 f_stoch.Ecrit_Valeurs ( [ m["NOM"], 4 , d["v_moy_loi"] , d["sigma_loi"] , d["v_min_loi"] , d["v_max_loi"] ] )
208 f_stoch.Ecrit_Titre ("Initial Points")
210 if m["POINT_INI"] is None :
211 aux = valeurs_lois[m]["v_moy_physique"]
214 f_stoch.Ecrit_Valeurs ( aux )
216 f_stoch.Ecrit_Titre ("Reference Points")
218 if m["POINT_REF"] is None :
219 aux = valeurs_lois[m]["v_moy_physique"]
222 f_stoch.Ecrit_Valeurs ( aux )
224 f_stoch.Ecrit_Titre ("Design Points")
226 if args["RECH_PT_CONCEPT"] == "OUI" :
228 elif m["POINT_CONCEPT"] is None :
229 aux = valeurs_lois[m]["v_moy_physique"]
231 aux = m["POINT_CONCEPT"]
232 f_stoch.Ecrit_Valeurs ( aux )
234 f_stoch.Ecrit_Titre ("Correlation matrix fictive")
235 for m in range(nb_occu_variable) :
237 for n in range(nb_occu_variable) :
238 aux.append(args["MATRICE"][n + m*nb_occu_variable])
239 f_stoch.Ecrit_Valeurs ( aux )
241 f_stoch.Ecrit_Titre ("Parameter threshold value")
242 if SEUIL_TYPE == "MAXIMUM" :
246 f_stoch.Ecrit_Valeurs ( aux )
248 # 4.4 ==> Fermeture du fichier
250 erreur = f_stoch.Ferme_Fichier ( )
255 erreur = f_stoch.Imprime_Fichier ( )
259 #____________________________________________________________________
261 # 5. Construction du fichier 'dataNum'
262 #____________________________________________________________________
264 # 5.1 ==> Ouverture du fichier
266 f_num = fiabilite_fichier.fiabilite_fichier ( self, Rep_Calc_LOGICIEL_global, fic_dataNum, INFO )
267 erreur = f_num.Ouvre_Fichier ( "w" )
271 # 5.2 ==> Ecriture des données nécessaires
273 f_num.Ecrit_Titre ("Parameters : EpsU, EpsG, Tau, Omega, iterMax")
274 if args["RECH_PT_CONCEPT"] == "OUI" :
275 f_num.Ecrit_Valeurs (args["EPSILON_U"])
276 f_num.Ecrit_Valeurs (args["EPSILON_G"])
277 f_num.Ecrit_Valeurs (args["TAU"])
278 f_num.Ecrit_Valeurs (args["OMEGA"])
279 f_num.Ecrit_Valeurs (args["ITER_MAX"])
283 f_num.Ecrit_Valeurs (aux)
285 f_num.Ecrit_Titre ("Parameters : hgrad, hhess")
286 f_num.Ecrit_Valeurs (args["HGRAD"])
287 f_num.Ecrit_Valeurs (args["HHESS"])
289 aux = [ "Parameter Optimality Test(sphere)" ]
290 aux.append("1: Parametric Method (Point Number in each direction)")
291 aux.append("2: Gaussian Method (Total Point Number)")
292 aux.append("3: Rejection Method (Total Point Number)")
293 f_num.Ecrit_Titre (aux)
295 if args["T_SPHERE"] == "OUI" :
297 if args["METHODE_TEST"] == "PARAMETRIQUE" :
299 elif args["METHODE_TEST"] == "GAUSSIENNE" :
301 elif args["METHODE_TEST"] == "REJECTION" :
304 self.cr.warn("METHODE DE TEST : "+args["METHODE_TEST"])
308 aux2 = args["NB_POINT"]
312 # remarque : il faut mettre une valeur plausible en aux1, sinon plantage violent ...
316 f_num.Ecrit_Valeurs ( [ aux1 , aux2 ] )
318 aux = [ "Parameters : alpha, beta" ]
319 aux.append("alpha: common net")
320 aux.append("beta: extreme net")
321 f_num.Ecrit_Titre (aux)
322 if args["PLAN_EXPERIENCE"] == "OUI" :
328 f_num.Ecrit_Valeurs ( aux1 )
329 f_num.Ecrit_Valeurs ( aux2 )
331 f_num.Ecrit_Titre ("Parameters Strong Max Test : cosLim, dProb")
332 if args["T_MAXIMUM_FORT"] == "OUI" :
333 aux1 = args["COS_LIM"]
338 f_num.Ecrit_Valeurs ( aux1 )
339 f_num.Ecrit_Valeurs ( aux2 )
341 f_num.Ecrit_Titre ("Parameter Importance Samplings : Simulation Number")
342 if args["TIRAGE_IMPORTANCE"] == "OUI" :
343 aux = args["NB_SIMULATION"]
346 f_num.Ecrit_Valeurs ( aux )
348 # 5.3 ==> Fermeture du fichier
350 erreur = f_num.Ferme_Fichier ( )
355 erreur = f_num.Imprime_Fichier ( )
359 #____________________________________________________________________
361 # 6. Construction du fichier 'dataGrad'
362 #____________________________________________________________________
364 # 6.1 ==> Création du fichier
366 f_grad = fiabilite_fichier.fiabilite_fichier ( self, Rep_Calc_LOGICIEL_global, fic_dataGrad, INFO )
367 erreur = f_grad.Ouvre_Fichier ( "w" )
371 # 6.2 ==> Ecriture des données nécessaires
373 f_grad.Ecrit_Titre ("Commentaires")
376 f_grad.Ecrit_Commentaires (m["NOM"])
377 if m["GRADIENT"] == "OUI" :
382 increment = m["INCREMENT"]
383 aux = [gradient,increment]
384 f_grad.Ecrit_Valeurs (aux)
386 # 6.3 ==> Fermeture du fichier
388 erreur = f_grad.Ferme_Fichier ( )
393 erreur = f_grad.Imprime_Fichier ( )
397 #____________________________________________________________________
400 #____________________________________________________________________
405 if not messages_erreur.has_key(erreur) :
407 self.cr.warn(messages_erreur[erreur])
412 ########################## Fin de la fonction##################################
414 ########################## Auto-test##################################
416 if __name__ == "__main__" :
422 Rep_Calc_LOGICIEL_global = tempfile.mktemp()
423 os.mkdir(Rep_Calc_LOGICIEL_global)
428 SEUIL_TYPE = "MAXIMUM"
433 erreur = fiabilite_mefisto ( None, Rep_Calc_LOGICIEL_global,
439 ### print "Erreur = ", erreur
440 Liste = os.listdir(Rep_Calc_LOGICIEL_global)
442 for nomfic in Liste :
443 fic_total = os.path.join(Rep_Calc_LOGICIEL_global,nomfic)
444 os.chmod (fic_total,0755)
445 os.remove (fic_total)
446 os.rmdir (Rep_Calc_LOGICIEL_global)