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. *
  29. * avec un chpoint
  30. chpi1 = ipol nua1 chp1 ;
  31. * on doit obtenir un chpo constant valant z = 1.52 ;
  32. un52 = extr chpi1 p1 'Z' ;
  33. si ( (abs (un52 - 1.52)) > 1.d-4 ) ;
  34. erreur 'Mauvaise interpolation';
  35. finsi ;
  36. *
  37. * avec un chamelem
  38. chmli1 = ipol nua1 chml1 ;
  39. * on doit obtenir un chml constant valant z = 1.5' ;
  40. un5 = extr chmli1 'Z' 1 1 1 ;
  41. si ( (abs (un5 - 1.5)) > 1.d-4 ) ;
  42. erreur 'Mauvaise interpolation';
  43. finsi ;
  44. *
  45. fin;
  46.  
  47.  
  48.  
  49.  

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