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

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