Télécharger maxi.dgibi

Retour à la liste

Numérotation des lignes :

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

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