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. *
  10. *
  11. * support géométrique
  12. p1 = 0. 0. ;
  13. p2 = 0. 1. ;
  14. p3 = 1. 1. ;
  15. p4 = 1. 0. ;
  16. qua1 = manu qua4 p1 p2 p3 p4 ;
  17. *
  18. * champs argumets
  19. chp1 = manu chpo p1 2 'X' 0.26 'Y' 0.26 ;
  20. mo1 = modele qua1 mecanique elastique ;
  21. chml1 = manu 'CHML' mo1 'X' 0.25 'Y' 0.25 ;
  22. *
  23. * le nuage stocke la fonction f(x,y)=z=x+y+1
  24. *
  25. nua1 = nuage 'X'*'FLOTTANT' 'Y'*'FLOTTANT' 'Z'*'FLOTTANT'
  26. 0. 0. 1. 1. 0. 2. 1. 1. 3. 0. 1. 2.
  27. 2. 3. 6. 0.5 0.5 2. 1.8 -0.2 2.6 1000. 0. 1001. ;
  28. fmax = 1000. ;
  29. tol1 = ('VALE' 'PREC') '*' 4. '*' fmax ;
  30. 'MESS' 'tol1 =' tol1 ;
  31. *
  32. * avec un chpoint
  33. chpi1 = ipol nua1 chp1 ;
  34. * on doit obtenir un chpo constant valant z = 1.52 ;
  35. un52 = extr chpi1 p1 'Z' ;
  36. si ( (abs (un52 - 1.52)) > tol1 ) ;
  37. erreur 'Mauvaise interpolation';
  38. finsi ;
  39. *
  40. * avec un chamelem
  41. chmli1 = ipol nua1 chml1 ;
  42. * on doit obtenir un chml constant valant z = 1.5' ;
  43. un5 = extr chmli1 'Z' 1 1 1 ;
  44. si ( (abs (un5 - 1.5)) > tol1 ) ;
  45. erreur 'Mauvaise interpolation';
  46. finsi ;
  47. *
  48. fin;
  49.  
  50.  
  51.  

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