Télécharger la plaquette Cast3M


fr en

Télécharger echi_som.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  echi_som.dgibi
  2. **** test de non-regression concernant l'operateur ECHI
  3.  
  4.  
  5.  
  6. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  7.  
  8. EPSI = 1.E-8 ;
  9.  
  10. P0 = 0.0 0.0 ;
  11. P1 = 1.0 0.0 ;
  12. P2 = 1.0 1.0 ;
  13. P3 = 0.0 1.0 ;
  14.  
  15. S1 = 'MANU' 'QUA4' P0 P1 P2 P3 ;
  16.  
  17. Mparoif = 'CHAN' S1 'QUAF' ;
  18. $paroif = 'MODE' Mparoif 'NAVIER_STOKES' 'LINE' ;
  19. paroif  = 'DOMA' $paroif 'MAILLAGE' ;
  20. Diag    = 'DOMA' $paroif 'XXDIAGSI' ;
  21.  
  22.    rtf = 'EQEX' 
  23.    'OPTI' 'EFM1' 'CENTREE' 'IMPL'
  24.    'ZONE' $paroif 'OPER' 'ECHI' 'KHEW'  'TBPW' 'INCO' 'TF'
  25.    'OPTI' 'EFM1' 'CENTREE' 'IMPL'
  26.    'ZONE' $paroif 'OPER' 'ECHI' 'KHW1'  'TBPW' 'INCO' 'TF'
  27.    'OPTI' 'EF' 'CENTREE' 'IMPL'
  28.    'ZONE' $paroif 'OPER' 'ECHI' 'KHEW'  'TBPW' 'INCO' 'TF'
  29.    'OPTI' 'EF' 'CENTREE' 'IMPL'
  30.    'ZONE' $paroif 'OPER' 'ECHI' 'KHW1'  'TBPW' 'INCO' 'TF' ;
  31.  
  32. rtf.'INCO' = 'TABLE' 'INCO' ;
  33. rtf.'INCO'.'TF'   = 'KCHT' $paroif 'SCAL' 'SOMMET' 10.0 ;
  34. rtf.'INCO'.'TBPW' = 'KCHT' $paroif 'SCAL' 'SOMMET'  5.0 ;
  35. XX1 = 'COORD' 1 paroif ;
  36. rtf.'INCO'.'KHEW' = 'KCHT' $paroif 'SCAL' 'CENTRE' 1.5 ;
  37. rtf.'INCO'.'KHW1' = 'KCHT' $paroif 'SCAL' 'SOMMET' (XX1 + 1.0 ) ;
  38.  
  39. B1 A1 = 'ECHI' rtf.'1ECHI' ;
  40. B2 A2 = 'ECHI' rtf.'2ECHI' ;
  41. B3 A3 = 'ECHI' rtf.'3ECHI' ;
  42. B4 A4 = 'ECHI' rtf.'4ECHI' ;
  43.  
  44. TF1 = 'COPIER' rtf.'INCO'.'TBPW' ;
  45. TF1 = 'NOMC' TF1 'TF' ;
  46.  
  47. BB2 = 'KMF' A2 TF1 ;
  48. BB4 = 'KMF' A4 TF1 ;
  49.  
  50. BB1 = 'KMF' A1 TF1 ;
  51. BB3 = 'KMF' A3 TF1 ;
  52.  
  53. ERROR = 0 ;
  54.  
  55. BBB2 = Diag * (rtf.'INCO'.'KHW1') * TF1 ;
  56.  
  57. 'SI' (('MAXI' (BB2 - B2) 'ABS') '>' EPSI) ;
  58. ERROR = ERROR + 1 ;
  59. 'MESS' ('MAXI' (BB2 - B2) 'ABS') ;
  60. 'MESS' 'Probleme echi (EFM1) coeff mult aux sommets !' ;
  61. 'FINSI' ;
  62.  
  63. 'SI' (('MAXI' (BBB2 - BB2) 'ABS') '>' EPSI) ;
  64. ERROR = ERROR + 1 ;
  65. 'MESS' ('MAXI' (BBB2 - BB2) 'ABS') ;
  66. 'MESS' 'Probleme echi (EFM1) coeff mult aux sommets !' ;
  67. 'FINSI' ;
  68.  
  69. 'SI' (('MAXI' (BB4 - B4) 'ABS') '>' EPSI) ;
  70. ERROR = ERROR + 1 ;
  71. 'MESS' ('MAXI' (BB4 - B4) 'ABS') ;
  72. 'MESS' 'Probleme echi (EF)   coeff mult aux sommets !' ;
  73. 'FINSI' ;
  74.  
  75. 'SI' (('MAXI' (BB1 - B1) 'ABS') '>' EPSI) ;
  76. ERROR = ERROR + 1 ;
  77. 'MESS' ('MAXI' (BB1 - B1) 'ABS') ;
  78. 'MESS' 'Probleme echi (EFM1) coeff mult aux centres !' ;
  79. 'FINSI' ;
  80.  
  81. 'SI' (('MAXI' (BB3 - B3) 'ABS') '>' EPSI) ;
  82. ERROR = ERROR + 1 ;
  83. 'MESS' ('MAXI' (BB3 - B3) 'ABS') ;
  84. 'MESS' 'Probleme echi (EF)   coeff mult aux centres !' ;
  85. 'FINSI' ;
  86.  
  87.  
  88. 'SI' (ERROR '>' 0) ;
  89. 'ERRE' 5 ;
  90. 'FINSI' ;
  91.  
  92.  
  93. 'FIN' ;
  94.  
  95.  
  96.  

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