Télécharger la plaquette Cast3M


fr en

Télécharger ipol1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ipol1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *   Fichier DGIBI pour tester l'operateur 'IPOL' avec un nuage
  6. *   la methode utilisée est celle des éléments finis diffus avec
  7. *   des polynomes d'ordre 1
  8. *   l'interpolation de toutes les fonctions lineaire est donc exacte
  9. *
  11. *
  12. * support géométrique
  13. p1 = 0. 0. ;
  14. p2 = 0. 1. ;
  15. p3 = 1. 1. ;
  16. p4 = 1. 0. ;
  17. qua1 = manu qua4 p1 p2 p3 p4 ;
  18. *
  19. * champs argumets
  20. chp1 = manu chpo p1 2 'X' 0.26 'Y' 0.26 ;
  21. mo1 = modele qua1 mecanique elastique ;
  22. chml1 = manu 'CHML' mo1 'X' 0.25 'Y' 0.25 ;
  23. *
  24. *  le nuage stocke la fonction f(x,y)=z=x+y+1
  25. *
  26. nua1 = nuage 'X'*'FLOTTANT' 'Y'*'FLOTTANT' 'Z'*'FLOTTANT'
  27.               0. 0. 1.   1. 0. 2.  1. 1. 3.   0. 1. 2.
  28.               2. 3. 6.  0.5 0.5 2.  1.8 -0.2 2.6  1000. 0. 1001. ;
  29. fmax = 1000.  ;
  30. tol1 = ('VALE' 'PREC') '*' 4. '*' fmax ;
  31. 'MESS' 'tol1 =' tol1 ;
  32. *
  33. *     avec un chpoint
  34. chpi1 = ipol nua1 chp1 ;
  35. * on doit obtenir un chpo constant valant z = 1.52 ;
  36. un52 = extr chpi1 p1 'Z' ;
  37. si ( (abs (un52 - 1.52))  > tol1 ) ;
  38.   erreur 'Mauvaise interpolation';
  39. finsi ;
  40. *
  41. *   avec un chamelem
  42. chmli1 = ipol nua1 chml1 ;
  43. * on doit obtenir un chml constant valant z = 1.5' ;
  44. un5 = extr chmli1 'Z' 1 1 1  ;
  45. si ( (abs (un5 - 1.5))  > tol1 ) ;
  46.   erreur 'Mauvaise interpolation';
  47. finsi ;
  48. *
  49. fin;
  50.  
  51.  
  52.  

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