Télécharger la plaquette Cast3M


fr en

Télécharger deadresi.procedur

Retour à la liste

Numérotation des lignes : 

  1. * DEADRESI  PROCEDUR  GOUNAND   26/01/09    21:15:16     12443          
  2. ************************************************************************
  3. * NOM         : DEADRESI
  4. * DESCRIPTION : Le résidu à annuler (Huang)
  5. *
  6. *
  7. *
  8. * LANGAGE     : GIBIANE-CAST3M
  9. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  10. *               mél : gounand@semt2.smts.cea.fr
  11. **********************************************************************
  12. * VERSION    : v1, 30/04/2007, version initiale
  13. * HISTORIQUE : v1, 30/04/2007, création
  14. * HISTORIQUE : 2018/01/22 : chgt nom composante hors diago G21 au lieu
  15. *                           de G12
  16. * HISTORIQUE : 2018/10/10 : NLINP au lieu de NLIN
  17. * HISTORIQUE :
  18. * HISTORIQUE :
  19. ************************************************************************
  20. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  21. * en cas de modification de ce sous-programme afin de faciliter
  22. * la maintenance !
  23. ************************************************************************
  24. *
  25. *
  26. 'DEBPROC' DEADRESI ;
  27. 'ARGUMENT' _mt*'MAILLAGE' ;
  28. 'ARGUMENT' theta*'FLOTTANT' ;
  29. 'ARGUMENT' gamma*'FLOTTANT' ;
  30. 'ARGUMENT' dpdua*'LISTMOTS' ;
  31. 'ARGUMENT' met/'CHPOINT' ;
  32. lmet = 'EXISTE' met ;
  33. idim  = 'VALEUR' 'DIME' ;
  34. 'ARGUMENT' idir/'ENTIER' ;
  35. 'SI' ('EXISTE' idir) ;
  36.    ldir = 'LECT' idir ;
  37. 'SINON' ;
  38.    'ARGUMENT' ldir/'LISTENTI' ;
  39.    'SI' ('NON' ('EXISTE' ldir)) ;
  40.       ldir = 'LECT' 1 'PAS' 1 idim ;
  41.    'FINSI' ;
  42. 'FINSI' ;
  43. *
  44. gdisc = DEADUTIL 'TYPM' _mt ;
  45. 'SI' ('EGA' gdisc 'LINE') ;
  46.    methgau = 'GAM1' ;
  47. 'SINON' ;
  48.    methgau = 'GAM2' ;
  49. 'FINSI' ;
  50. vdim  = DEADUTIL 'DIMM' _mt ;
  51. *
  52. debloi = 'CHAINE' 'AHU' ;
  53. ncmet  = '/' ('*' idim ('+' idim 1)) 2 ;
  54. loi = 'CHAINE' debloi 'R' ;
  55. *
  56. *
  57. * Calcul du résidu
  58. *
  59. nldir = 'DIME' ldir ;
  60. *
  61. numop  = '*' nldir vdim ;
  62. numder = vdim ;
  63. numvar = 1 ;
  64. numdat = 0 ;
  65. numcof = 0 ;
  66. *
  67.    A = ININLIN numop numvar numdat numcof numder ;
  68.    A . 'VAR' . 1 . 'VALEUR' = 1.D0 ;
  69. *
  70.    numvar = idim ;
  71.    numdat = '+' ncmet 2 ;
  72.    numcof = '*' nldir vdim ;
  73.    B = ININLIN numop numvar numdat numcof numder ;
  74.    'REPETER' ivar numvar ;
  75.       B . 'VAR' . &ivar . 'NOMDDL' = 'EXTRAIRE' dpdua &ivar ;
  76.       B . 'VAR' . &ivar . 'DISC'   = gdisc ;
  77.    'FIN' ivar ;
  78. *
  79.   idat = 0 ;
  80.    'REPETER' idi idim ;
  81.       nomdat = 'CHAINE' 'G' &idi &idi ;
  82. *         'MESSAGE' ('CHAINE' 'nomdat=' nomdat) ;
  83.       idat = '+' idat 1 ;
  84.       B . 'DAT' . idat . 'NOMDDL' = nomdat ;
  85.       B . 'DAT' . idat . 'DISC'   = gdisc ;
  86.       'SI' lmet ;
  87. *         B . 'DAT' . idat . 'VALEUR' = 'EXCO' nomdat met nomdat ;
  88.          B . 'DAT' . idat . 'VALEUR' = met ;
  89.       'SINON' ;
  90.          B . 'DAT' . idat . 'VALEUR' = 1.D0 ;
  91.       'FINSI' ;
  92.    'FIN' idi ;
  93.    'REPETER' idi idim ;
  94.       nj = '-' idim &idi ;
  95. *      'MESSAGE' ('CHAINE' 'nj=' nj) ;
  96.       'SI' ('>EG' nj 1) ;
  97.          'REPETER' jdi nj ;
  98. * Mise en cohérence avec Castem
  99. *            nomdat = 'CHAINE' 'G' &idi ('+' &idi &jdi) ;
  100.             nomdat = 'CHAINE' 'G' ('+' &idi &jdi) &idi ;
  101. *               'MESSAGE' ('CHAINE' 'nomdat=' nomdat) ;
  102.             idat = '+' idat 1 ;
  103.             B . 'DAT' . idat . 'NOMDDL' = nomdat ;
  104.             B . 'DAT' . idat . 'DISC'   = gdisc ;
  105.             'SI' lmet ;
  106. *               B . 'DAT' . idat . 'VALEUR' = 'EXCO' nomdat met nomdat ;
  107.                B . 'DAT' . idat . 'VALEUR' = met ;
  108.             'SINON' ;
  109.                B . 'DAT' . idat . 'VALEUR' = 0.D0 ;
  110.             'FINSI' ;
  111.          'FIN' jdi ;
  112.       'FINSI' ;
  113.    'FIN' idi ;
  114.    lisdat = 'LECT' 1 'PAS' 1 'NPAS' ('-' numdat 1) ;
  115.    idat = '+' idat 1 ;
  116.    B . 'DAT' . idat . 'VALEUR' = theta ;
  117.    idat = '+' idat 1 ;
  118.    B . 'DAT' . idat . 'VALEUR' = gamma ;
  119. *
  120.    iop = 0 ;
  121.    'REPETER' k nldir ;
  122.       idir = 'EXTRAIRE' ldir &k ;
  123.       'REPETER' l vdim ;
  124.          iop = '+' iop 1 ;
  125.          A . iop . 1 . 0  = 0 ;
  126.          nomloi = 'CHAINE' loi idir &l ;
  127. *         'MESSAGE' ('CHAINE' 'Nomloi=' nomloi) ;
  128.          B . 'COF' . iop . 'COMPOR' = nomloi ;
  129.          B . 'COF' . iop . 'LDAT'   = lisdat ;
  130.          B . iop . idir . &l  = iop ;
  131.       'FIN' l ;
  132.    'FIN' k ;
  133. *
  134. res = NLINP gdisc _mt A B 'EREF' methgau ;
  135. *
  136. 'RESPRO' res ;
  137. *
  138. * End of procedure file DEADRESI
  139. *
  140. 'FINPROC' ;
  141.  
  142.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales