Télécharger la plaquette Cast3M


fr en

Télécharger maxi.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : maxi.dgibi
  2.  
  3. ************************************************************************
  4. ************************************************************************
  5. ************************************************************************
  6. * 
  7. * Mots-clés : MAXI, MANU 'CHPO'
  8. * 
  9. * Test de verification de les operateurs MAXI et MANU 'CHPO'
  10. * 
  11. * creation : bp, 2015-06-09
  12. * modifs   : ... a completer ...
  13. * 
  14. ************************************************************************
  15.  
  16. OPTI DIME 2 ELEM QUA4;
  17. NBERR = 0; lerr = lect;
  18. XPREC = 1.E-10;
  19.  
  20. ************************************************************************
  21. * syntaxe de MAXI : 
  22. * Xmax = MAXI X1 X2 ... Xn
  23. ************************************************************************
  24. X1 = MAXI -10. 1. 0. 2.5 7. -5.5; 
  25. X2 = MAXI 'ABS' -10. 1. 0. 2.5 7. -5.5; 
  26. X3 = MINI  -10. 1. 0. 2.5 7. -5.5; 
  27. X4 = MAXI  -10 1 0 2 7 -5; 
  28. X5 = MAXI 'ABS' -10 1 0 2 7 -5; 
  29. X6 = MINI  -10 1 0 2 7 -5; 
  30. X1a6 = PROG X1  X2  X3  X4  X5 X6;
  31. V1a6 = PROG 7. 10. -10.  7  10 -10;
  32. Xerr = MAXI 'ABS' (X1a6 - V1a6);
  33. SI (Xerr > XPREC);
  34.   NBERR=NBERR + 1; lerr= lerr et 0;
  35. FINSI;
  36.  
  37. ************************************************************************
  38. * syntaxe de MAXI : 
  39. * OBJET3   =  MAXI  OBJET1  OBJET2 (objeti ..)
  40. ************************************************************************
  41.  
  42. *** LISTREEL ***
  43.  
  44. pipi = prog  1 2 3 4;                                                   
  45. pipi2 = 2.* pipi;                                                       
  46. pipi3 = 3.* pipi;  
  47.  
  48. popo = maxi pipi pipi2;                                                 
  49. popo3 = maxi pipi pipi2 pipi3;
  50.  
  51. * list popo;
  52. xpopo = SOMM popo;
  53. SI ((ABS (xpopo - 20.)) > XPREC); 
  54.   NBERR=NBERR + 1; lerr= lerr et 1;
  55. FINSI;
  56. * list popo3;
  57. xpopo3 = SOMM popo3;
  58. SI ((ABS (xpopo3 - 30.)) > XPREC); 
  59.   NBERR=NBERR + 1; lerr= lerr et 2; 
  60. FINSI;
  61.  
  62.  
  63. *** CHPOINT ***
  64.  
  65. *************** créé via MANU CHPO + maillage de POI1 :
  66.  
  67. p0 = 0. 0.;
  68. mesh1 = MANU 'POI1' p0 ;
  69. kr = MANU 'CHPO' mesh1 1 scal 0.19;
  70. va = MANU 'CHPO' mesh1 1 scal 0.10;
  71. * list kr ;
  72. * list va ;
  73. * on peut aussi verifier ce qu il y a comme maillage
  74. * kr1 = extr  kr 'MAIL'; list kr1;
  75. * va1 = extr  va 'MAIL'; list va1;
  76. obj = MAXI va kr ;
  77. * list obj;
  78. xobj = EXTR obj 'SCAL' p0;
  79. SI ((ABS (xobj - 0.19)) > XPREC); 
  80.   NBERR=NBERR + 1; lerr= lerr et 3; 
  81. FINSI;
  82.  
  83. *************** créé via COOR + maillage d elements
  84.  
  85. p1 = 0 0; 
  86. p2 = 2 1;
  87. l12 = p1 droi 10 p2;
  88. x = coor l12 1;
  89. list x;
  90. * meshx = extr x MAIL; list meshx;
  91. x2 = x**2;
  92. x3 = x**3;
  93. toto =  maxi x x2;
  94. list toto;
  95.  
  96. toto3 =  maxi x x2 x3;
  97. list toto3;
  98.  
  99. * test
  100. xtoto = MAXI  toto  'ABS';
  101. xtoto3 = MAXI toto3 'ABS';
  102. SI ((ABS (xtoto  - 4.)) > XPREC); 
  103.   NBERR=NBERR + 1; lerr= lerr et 4; 
  104. FINSI;
  105. SI ((ABS (xtoto3 - 8.)) > XPREC); 
  106.   NBERR=NBERR + 1; lerr= lerr et 5; 
  107. FINSI;
  108. xtoto = MAXI (RESU  toto)  ;
  109. xtoto3 = MAXI (RESU toto3) ;
  110. SI ((ABS (xtoto  - 16.2)) > XPREC); 
  111.   NBERR=NBERR + 1; lerr= lerr et 6; 
  112. FINSI;
  113. SI ((ABS (xtoto3 - 25.4)) > XPREC); 
  114.   NBERR=NBERR + 1; lerr= lerr et 7; 
  115. FINSI;
  116.  
  117. * apres un REDU (le meleme change)
  118. pmil = l12 poin 6;
  119. d1mil = l12 elem compris p1 pmil;
  120. titi = maxi (redu x d1mil) (redu x2 d1mil);
  121. xtiti = MAXI (RESU  titi)  ;
  122. SI ((ABS (xtiti  - 3.0)) > XPREC); 
  123.   NBERR=NBERR + 1; lerr= lerr et 8;
  124. FINSI;
  125.  
  126.  
  127.  
  128. *************** créé via MANU CHPO + maillage d elements :
  129.  
  130. kr2 = MANU 'CHPO' l12 2 'UX' 0.19 'UY'  0.29;
  131. va2 = MANU 'CHPO' l12 2 'UX' 0.10 'UY' -0.45;
  132. obj2 = MAXI va2 kr2 ;
  133. xobj2 = EXTR obj2 'UY' p2;
  134. SI ((ABS (xobj2 - 0.29)) > XPREC); 
  135.   NBERR=NBERR + 1; lerr= lerr et 9; 
  136. FINSI;
  137.  
  138.  
  139. *************** créé via MANU CHPO + maillage de POI1 avec des doublons:
  140.  
  141. mesh11 = (MANU 'POI1' p0) et p1 'TELQUEL';
  142. mesh11 = mesh11 et p2 'TELQUEL';
  143. mesh11 = mesh11 et p0 'TELQUEL';
  144. kr1 = MANU 'CHPO' mesh11 1 scal 0.19 'NATURE' 'DIFFUS';
  145. va1 = MANU 'CHPO' mesh11 1 scal 0.10 'NATURE' 'DIFFUS';
  146. max11 = MAXI kr1 va1;
  147. xmax11 = MAXI max11;
  148. SI ((ABS (xmax11 - 0.19)) > XPREC); 
  149.   NBERR=NBERR + 1; lerr= lerr et 10; 
  150. FINSI;
  151.  
  152.  
  153. opti echo 0;
  155. mess '>>>>> ' NBERR ' ERREURS DETECTEES  <<<<<' ;
  156. mess 'pour les cas tests :'; list lerr;
  158. SI (NBERR EGA 0); ERRE 0;
  159. SINON;            ERRE 5;
  160. FINSI;
  161.  
  162. FIN ;
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  

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