Télécharger la plaquette Cast3M


fr en

Télécharger projgril_1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : projgril_1.dgibi
  2. OPTI 'DIME' 2 'ECHO' 1 ;
  3.  
  4. ************************************************************************
  5. *                    Test de la procedure PROJGRIL                     *
  6. *    Projection dans 2 dimensions d'un nuage representant une grille   *
  7. *                           de n dimensions                            *
  8. *        utilise dans le cas de l'operateur IPOL option 'GRILL'        *
  9. *                                                                      *
  10. * - test avec fonction de 2, 3, 4 et meme 5 variables                  *
  11. ************************************************************************
  12. MESS ; MESS ; MESS ;
  13.  
  14.  
  15.  
  16.  
  17. ** Tolerance pour les tests et indicateur de trace
  18. ** -----------------------------------------------
  19. tol1 = 1.E-14 ;
  20. itrac = FAUX ;
  21.  
  22.  
  23.  
  24.  
  25. ** TEST sur une grille de dimension 2
  26. ** ----------------------------------
  27. MESS 'Grille de dimension 2' ;
  29. l1  = PROG 1. 5. 12. ;
  30. l2  = PROG 0. 2.5 ;
  31. lf  = PROG 4. 8. 15.   16. 23. 42. ;
  32. nu1 = NUAG 'COMP' 'JIN' l1
  33.            'COMP' 'SUN' l2
  34.            'COMP' 'JACK' lf ;
  35. mail1 chp1 = PROJGRIL nu1 ;
  36. tp = TABL ;
  37. tf = TABL ;
  38. tp . 1 = 1. 0. ;
  39. tf . 1 = 4. ;
  40. tp . 2 = 5. 0. ;
  41. tf . 2 = 8. ;
  42. tp . 3 = 12. 0. ;
  43. tf . 3 = 15. ;
  44. tp . 4 = 1. 2.5 ;
  45. tf . 4 = 16. ;
  46. tp . 5 = 5. 2.5 ;
  47. tf . 5 = 23. ;
  48. tp . 6 = 12. 2.5 ;
  49. tf . 6 = 42. ;
  50. MESS ' Point  |   Valeur       |   Valeur       | Erreur' ;
  51. MESS '        |   theorique    |   interpolee   |       ' ;
  52. MESS '--------|----------------|----------------|-------' ;
  53. REPE b1 (DIME tp) ;
  54.   ft = tf . &b1 ;
  55.   pt = mail1 POIN 'PROC' (tp . &b1) ;
  56.   fc = EXTR chp1 'JACK' pt ;
  57.   err1 = ABS ((fc - ft) / ft) ;
  58.   MESS &b1 '|' ft '|' fc '|' err1 ;
  59.   SI (err1 > tol1) ;
  60.     MESS ; MESS ;
  61.     MESS 'ECHEC DU CAS TEST !' ;
  62.     ERREUR 4 ;
  63.   FINSI ;
  64. FIN b1 ;
  65. SI itrac ;
  66.   TRAC chp1 mail1 'TITR' 'Test sur une grille en dimension 2' ;
  67. FINSI ;
  68. MESS ; MESS ; MESS ;
  69.  
  70.  
  71.  
  72.  
  73. ** TEST sur une grille de dimension 3
  74. ** ----------------------------------
  75. MESS 'Grille de dimension 3' ;
  77. l3 = PROG 0. 10. ;
  78. lf = lf ET (PROG 8. 42. 4.   16. 15. 23.) ;
  79. nu1 = NUAG 'COMP' 'EKO' l1
  80.            'COMP' 'BEN' l2
  81.            'COMP' 'BOON' l3
  82.            'COMP' 'KATE' lf ;
  83. mail1 chp1 = PROJGRIL nu1 (MOTS 'BEN') (PROG 1.25) ;
  84. tp = TABL ;
  85. tf = TABL ;
  86. tp . 1 = 1. 0. ;
  87. tf . 1 = 10. ;
  88. tp . 2 = 5. 0. ;
  89. tf . 2 = 15.5 ;
  90. tp . 3 = 12. 0. ;
  91. tf . 3 = 28.5 ;
  92. tp . 4 = 1. 10. ;
  93. tf . 4 = 12. ;
  94. tp . 5 = 5. 10. ;
  95. tf . 5 = 28.5 ;
  96. tp . 6 = 12. 10. ;
  97. tf . 6 = 13.5 ;
  98. MESS ' Point  |   Valeur       |   Valeur       | Erreur' ;
  99. MESS '        |   theorique    |   interpolee   |       ' ;
  100. MESS '--------|----------------|----------------|-------' ;
  101. REPE b1 (DIME tp) ;
  102.   ft = tf . &b1 ;
  103.   pt = mail1 POIN 'PROC' (tp . &b1) ;
  104.   fc = EXTR chp1 'KATE' pt ;
  105.   err1 = ABS ((fc - ft) / ft) ;
  106.   MESS &b1 '|' ft '|' fc '|' err1 ;
  107.   SI (err1 > tol1) ;
  108.     MESS ; MESS ;
  109.     MESS 'ECHEC DU CAS TEST !' ;
  110.     ERREUR 4 ;
  111.   FINSI ;
  112. FIN b1 ;
  113. SI itrac ;
  114.   TRAC chp1 mail1 'TITR' 'Test sur une grille en dimension 3' ;
  115. FINSI ;
  116. MESS ; MESS ; MESS ;
  117.  
  118.  
  119.  
  120. ** TEST sur une grille de dimension 4
  121. ** ----------------------------------
  122. MESS 'Grille de dimension 4' ;
  124. l4 = PROG 1977. 2004. ;
  125. lf = lf ET (PROG 42.  4. 23.   8. 15. 16.
  126.                   4. 23. 16.   8. 15. 42.) ;
  127. nu1 = NUAG 'COMP' 'WALT' l1
  128.            'COMP' 'ROSE' l2
  129.            'COMP' 'JACO' l3
  130.            'COMP' 'SAYI' l4
  131.            'COMP' 'HUGO' lf ;
  132. mail1 chp1 = PROJGRIL nu1 (MOTS 'ROSE' 'WALT') (PROG 1.25 8.5) ;
  133. tp = TABL ;
  134. tf = TABL ;
  135. tp . 1 = 0. 1977. ;
  136. tf . 1 = 22. ;
  137. tp . 2 = 10. 1977. ;
  138. tf . 2 = 21. ;
  139. tp . 3 = 0. 2004. ;
  140. tf . 3 = 14.5 ;
  141. tp . 4 = 10. 2004. ;
  142. tf . 4 = 24. ;
  143. MESS ' Point  |   Valeur       |   Valeur       | Erreur' ;
  144. MESS '        |   theorique    |   interpolee   |       ' ;
  145. MESS '--------|----------------|----------------|-------' ;
  146. REPE b1 (DIME tp) ;
  147.   ft = tf . &b1 ;
  148.   pt = mail1 POIN 'PROC' (tp . &b1) ;
  149.   fc = EXTR chp1 'HUGO' pt ;
  150.   err1 = ABS ((fc - ft) / ft) ;
  151.   MESS &b1 '|' ft '|' fc '|' err1 ;
  152.   SI (err1 > tol1) ;
  153.     MESS ; MESS ;
  154.     MESS 'ECHEC DU CAS TEST !' ;
  155.     ERREUR 4 ;
  156.   FINSI ;
  157. FIN b1 ;
  158. SI itrac ;
  159.   TRAC chp1 mail1 'TITR' 'Test sur une grille en dimension 4' ;
  160. FINSI ;
  161. MESS ; MESS ; MESS ;
  162.  
  163.  
  164.  
  165.  
  166. ** TEST sur une grille de dimension 5
  167. ** ----------------------------------
  168. MESS 'Grille de dimension 5' ;
  170. l5 = PROG 1. 4. 6. ;
  171. lf = lf ET (PROG  4. 23. 15.    4.  8. 15.   4.  8. 15.   42. 23. 16.
  172.                   8. 42. 16.   42. 23. 16.  16. 23. 16.   15.  8. 15.
  173.                  15.  4. 23.    4.  8. 15.  15. 42. 23.   16.  4.  8.
  174.                  16.  8. 42.   42. 23. 16.   8.  4. 42.   23. 42.  4.) ;
  175. nu1 = NUAG 'COMP' 'ALEX' l1
  176.            'COMP' 'ANA' l2
  177.            'COMP' 'SAWY' l3
  178.            'COMP' 'CHAN' l4
  179.            'COMP' 'DESM' l5
  180.            'COMP' 'JOHN' lf ;
  181. mail1 chp1 = PROJGRIL nu1 (MOTS 'ALEX' 'CHAN' 'SAWY' )
  182.                           (PROG   5.    1984.     10.   ) ;
  183. xb = (1984. - 1977.) / (2004. - 1977.) ;
  184. xbb = 1. - xb ;
  185. tp = TABL ;
  186. tf = TABL ;
  187. tp . 1 = 0. 1. ;
  188. tf . 1 = (42. * xbb) + (23. * xb) ;
  189. tp . 2 = 2.5 1. ;
  190. tf . 2 = 15. ;
  191. tp . 3 = 0. 4. ;
  192. tf . 3 = (8. * xbb) + (23. * xb) ;
  193. tp . 4 = 2.5 4. ;
  194. tf . 4 = (23. * xbb) + (8. * xb) ;
  195. tp . 5 = 0. 6. ;
  196. tf . 5 = (42. * xbb) + (4. * xb) ;
  197. tp . 6 = 2.5 6. ;
  198. tf . 6 = (4. * xbb) + (42. * xb) ;
  199. MESS ' Point  |   Valeur       |   Valeur       | Erreur' ;
  200. MESS '        |   theorique    |   interpolee   |       ' ;
  201. MESS '--------|----------------|----------------|-------' ;
  202. REPE b1 (DIME tp) ;
  203.   ft = tf . &b1 ;
  204.   pt = mail1 POIN 'PROC' (tp . &b1) ;
  205.   fc = EXTR chp1 'JOHN' pt ;
  206.   err1 = ABS ((fc - ft) / ft) ;
  207.   MESS &b1 '|' ft '|' fc '|' err1 ;
  208.   SI (err1 > tol1) ;
  209.     MESS ; MESS ;
  210.     MESS 'ECHEC DU CAS TEST !' ;
  211.     ERREUR 4 ;
  212.   FINSI ;
  213. FIN b1 ;
  214. SI itrac ;
  215.   TRAC chp1 mail1 'TITR' 'Test sur une grille en dimension 5' ;
  216. FINSI ;
  217. MESS ; MESS ; MESS ;
  218.  
  219.  
  220.  
  221.  
  223. MESS 'SUCCES DU CAS TEST' ;
  224.  
  225.  
  226.  
  227.  
  228. FIN ;
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  

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