Aster/Cata/cataSTA74/Macro/impr_table_ops.py

   1 #@ MODIF impr_table_ops Macro  DATE 30/11/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.
  10 #
  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.
  15 #
  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 # ======================================================================
  20
  21 # RESPONSABLE MCOURTOI M.COURTOIS
  22
  23 import os.path
  24 import re
  25
  26 from types import ListType, TupleType, StringTypes
  27 EnumTypes=(ListType, TupleType)
  28
  29
  30 # ------------------------------------------------------------------------------
  31 def impr_table_ops(self, FORMAT, TABLE, INFO, **args):
  32    """
  33    Macro IMPR_TABLE permettant d'imprimer une table dans un fichier.
  34    Erreurs<S> dans IMPR_TABLE pour ne pas perdre la base.
  35    """
  36    macro='IMPR_TABLE'
  37    import aster
  38    from Accas import _F
  39    from Cata.cata import table_jeveux
  40    from Utilitai.Utmess import UTMESS
  41    ier=0
  42    # La macro compte pour 1 dans la numerotation des commandes
  43    self.set_icmd(1)
  44
  45    # On importe les definitions des commandes a utiliser dans la macro
  46    # Le nom de la variable doit etre obligatoirement le nom de la commande
  47    DETRUIRE         = self.get_cmd('DETRUIRE')
  48    DEFI_FICHIER     = self.get_cmd('DEFI_FICHIER')
  49    RECU_FONCTION    = self.get_cmd('RECU_FONCTION')
  50
  51    #----------------------------------------------
  52    # 0. Traitement des arguments, initialisations
  53    # unité logique des fichiers réservés
  54    ul_reserve=(8,)
  55
  56    # 0.1. Fichier
  57    nomfich=None
  58    if args['UNITE'] and args['UNITE']<>6:
  59       nomfich='fort.'+str(args['UNITE'])
  60    if nomfich and os.path.exists(nomfich):
  61       if FORMAT=='XMGRACE':
  62          UTMESS('A',macro,'Le fichier '+nomfich+' existe déjà, on écrit ' \
  63                 'à la suite.')
  64
  65    # 0.2. Création des dictionnaires des FILTRES
  66    Filtre=[]
  67    if args['FILTRE']:
  68       for Fi in args['FILTRE']:
  69          dF = Fi.cree_dict_valeurs(Fi.mc_liste)
  70          for mc in dF.keys():
  71             if dF[mc]==None: del dF[mc]
  72          Filtre.append(dF)
  73    # format pour l'impression des filtres
  74    form_filtre='\nFILTRE -> NOM_PARA: %-16s CRIT_COMP: %-4s VALE: %s'
  75
  76    # 0.3. Création de la liste des tables (une seule sans SENSIBILITE)
  77    form_sens='\n... SENSIBILITE AU PARAMETRE %s'
  78    ltab=[]
  79    if args['SENSIBILITE']:
  80       nmemo='&NOSENSI.MEMO.CORR'.ljust(24)
  81       vect=aster.getvectjev(nmemo)
  82       if vect:
  83          lps=args['SENSIBILITE']
  84          if not type(lps) in EnumTypes:
  85             lps=[lps,]
  86          for ps in [ps.get_name() for ps in lps]:
  87             trouv=False
  88             for ch in vect[0:len(vect):2]:
  89                if ch[0:8].strip()==TABLE.get_name() and ch[8:16].strip()==ps:
  90                   trouv=True
  91                   ncomp=ch[16:24].strip()
  92                   sdtab=table_jeveux(ncomp)
  93                   tabs=sdtab.EXTR_TABLE()
  94                   tabs.titr+=form_sens % ps
  95                   ltab.append([tabs, sdtab])
  96             if not trouv:
  97                UTMESS('A',macro,'Dérivée de %s par rapport à %s non disponible'\
  98                      % (TABLE.get_name(), ps))
  99       else:
 100          UTMESS('A',macro,'Pas de calcul de sensibilité accessible.')
 101    else:
 102       ltab.append([TABLE.EXTR_TABLE(), TABLE])
 103
 104    if len(ltab)<1:
 105       return ier
 106
 107    # 0.4.1. liste des paramètres à conserver
 108    nom_para=ltab[0][0].para
 109    if args['NOM_PARA']:
 110       nom_para=args['NOM_PARA']
 111
 112    # 0.4.2. Traiter le cas des UL réservées
 113    if args['UNITE'] and args['UNITE'] in ul_reserve:
 114       DEFI_FICHIER( ACTION='LIBERER', UNITE=args['UNITE'], )
 115
 116    #----------------------------------------------
 117    # Boucle sur les tables
 118    for tab, sdtab in ltab:
 119
 120       # ----- 1. Infos de base
 121       if INFO==2:
 122          print 'IMPRESSION DE LA TABLE : %s' % sdtab.get_name()
 123
 124       if args['TITRE']:
 125          tab.titr=args['TITRE'] + '\n' + tab.titr
 126
 127       # ----- 2. Filtres
 128       for Fi in Filtre:
 129          col = getattr(tab, Fi['NOM_PARA'])
 130          # peu importe le type
 131          opts=[Fi[k] for k in ('VALE','VALE_I','VALE_C','VALE_K') if Fi.has_key(k)]
 132          kargs={}
 133          for k in ('CRITERE','PRECISION'):
 134             if Fi.has_key(k):
 135                kargs[k]=Fi[k]
 136          tab = tab & ( getattr(col, Fi['CRIT_COMP'])(*opts,**kargs) )
 137          # trace l'operation dans le titre
 138          #if FORMAT in ('TABLEAU','ASTER'):
 139          tab.titr+=form_filtre % (Fi['NOM_PARA'], Fi['CRIT_COMP'], \
 140                ' '.join([str(v) for v in opts]))
 141
 142       # ----- 3. Tris
 143       if args['TRI']:
 144          # une seule occurence de TRI
 145          T0=args['TRI'][0]
 146          dT=T0.cree_dict_valeurs(T0.mc_liste)
 147          tab.sort(CLES=dT['NOM_PARA'], ORDRE=dT['ORDRE'])
 148
 149       # ----- 4. Impression
 150       timp=tab[nom_para]
 151       # passage des mots-clés de mise en forme à la méthode Impr
 152       kargs=args.copy()
 153       kargs.update({
 154          'FORMAT'    : FORMAT,
 155          'FICHIER'   : nomfich,
 156          'dform'     : {},
 157       })
 158       # pour l'impression des fonctions
 159       kfonc={
 160          'FORMAT'    : FORMAT,
 161          'FICHIER'   : nomfich,
 162       }
 163
 164       # 4.1. au format AGRAF
 165       if FORMAT=='AGRAF':
 166          kargs['dform']={ 'formR' : '%12.5E' }
 167          kfonc['FORMAT']='TABLEAU'
 168
 169       # 4.2. au format XMGRACE et dérivés
 170       elif FORMAT=='XMGRACE':
 171          kargs['dform']={ 'formR' : '%.8g' }
 172          kargs['PILOTE']=args['PILOTE']
 173          kfonc['PILOTE']=args['PILOTE']
 174
 175       # 4.3. format spécifié dans les arguments
 176       if args['FORMAT_R']:
 177          kargs['dform'].update({ 'formR' : fmtF2PY(args['FORMAT_R']) })
 178
 179       # 4.4. regroupement par paramètre : PAGINATION
 180       if args['PAGINATION']:
 181          kargs['PAGINATION']=args['PAGINATION']
 182
 183       timp.Impr(**kargs)
 184
 185       # ----- 5. IMPR_FONCTION='OUI'
 186       if args['IMPR_FONCTION'] and args['IMPR_FONCTION']=='OUI':
 187          # cherche parmi les cellules celles qui contiennent un nom de fonction
 188          dfon={}
 189          for row in timp:
 190             for par,cell in row.items():
 191                if type(cell) in StringTypes:
 192                 if aster.getvectjev(cell.strip().ljust(19)+'.PROL')<>None:
 193                   dfon[cell.strip().ljust(19)]=par
 194          # impression des fonctions trouvées
 195          for f,par in dfon.items():
 196             __fonc=RECU_FONCTION(
 197                TABLE=sdtab,
 198                FILTRE=_F(
 199                   NOM_PARA=par,
 200                   VALE_K=f,
 201                ),
 202                NOM_PARA_TABL=par,
 203             )
 204             __fonc.Trace(**kfonc)
 205             DETRUIRE(CONCEPT=_F(NOM=('__fonc',),), ALARME='NON', INFO=1,)
 206
 207    # 99. Traiter le cas des UL réservées
 208    if args['UNITE'] and args['UNITE'] in ul_reserve:
 209       DEFI_FICHIER( ACTION='ASSOCIER', UNITE=args['UNITE'],
 210             TYPE='ASCII', ACCES='APPEND' )
 211
 212    return ier
 213
 214 # ------------------------------------------------------------------------------
 215 def fmtF2PY(fformat):
 216    """Convertit un format Fortran en format Python (printf style).
 217    Gère uniquement les fortrans réels, par exemple : E12.5, 1PE13.6, D12.5...
 218    """
 219    fmt=''
 220    matP=re.search('([0-9]+)P',fformat)
 221    if matP:
 222       fmt+=' '*int(matP.group(1))
 223    matR=re.search('([eEdDfFgG]{1})([\.0-9]+)',fformat)
 224    if matR:
 225       fmt+='%'+matR.group(2)+re.sub('[dD]+','E',matR.group(1))
 226    try:
 227       s=fmt % -0.123
 228    except (ValueError, TypeError), msg:
 229       fmt='%12.5E'
 230       print 'Error :',msg
 231       print 'Format par défaut utilisé :',fmt
 232    return fmt