Télécharger la plaquette Cast3M


fr en

Télécharger proi-parallele.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : proi-parallele.dgibi
  2. 'OPTI' echo 0 ;
  3.  
  4. ************************************************************************
  5. * NOM         : PROI-PARALLELE
  6. * DESCRIPTION : On teste le parallélisme avec les assistants pour faire
  7. *               un PROI en parallèle
  8. *
  9. *
  10. *
  11. * LANGAGE     : GIBIANE-CAST3M
  12. * AUTEUR      : Stephane GOUNAND (CEA/DEN/DM2S/SEMT/LTA)
  13. *               mel : stephane.gounand@cea.fr
  14. **********************************************************************
  15. * VERSION    : v1, 31/07/2019, version initiale
  16. * HISTORIQUE : v1, 31/07/2019, création
  17. * HISTORIQUE :
  18. * HISTORIQUE :
  19. ************************************************************************
  20. *
  21. *
  22. interact= FAUX ;
  23. *
  24. 'OPTION' 'DIME' 2 'ELEM' 'QUA4' ;
  25. *
  26. p1 = 0. 0.   ; p2 = 1.5 0.1 ;
  27. p3 = 1.4 1.3 ; p4 = 0.1 1.3 ; 
  28. *
  29. n1 = 3 ; n2 = 4 ; n3 = 5 ; n4 = 6 ;
  30. cnt = 'DROI' n1 p1 p2 'DROI' n2 p3 'DROI' n3 p4 'DROI' n4 p1 ;
  31. sur = 'SURF' cnt ;
  32. *
  33. x y = 'COOR' sur ;
  34. r = '**' ('+' ('**' x 2) ('**' y 2)) 0.5 ;
  35. chamr = 'CHAN' 'CHAM' r sur ;
  36. *
  37. lok = vrai ;
  38. *
  39. * PROI sur les points de sur
  40. *
  41. mail = 'CHAN' 'POI1' sur ;
  42. *
  43. * Sequentiel
  44. *
  45. chpors   = 'PROI' chamr mail ;
  46. errs = 'MAXI' ('-' chpors r) 'ABS' ;
  47. * l'erreur vaut 5.d-9 sur ma machine au 31/07/2019
  48. * l'erreur vaut 5.44493E-08 sur semt2 au 01/08/2019
  49. crit1 = 1.D-7 ;
  50. 'MESS' 'errs=' errs ' crit1=' crit1 ;
  51. tst1 = errs '<' crit1 ;
  52. lok = lok 'ET' tst1 ;
  53. *
  54. * Parallele
  55. *
  56. nbpart = 'VALEUR' 'ASSI' ;
  57. mapart = 'PART' 'ARLE'  mail nbpart ;
  58. tchpo  = 'ASSI' 'TOUS' 'PROI' chamr mapart ;
  59. *sg 23/03/2016 recommande par Clement
  60. *      chpo   = 'ET' tchpo ;
  61. chporp   = 'ETG' tchpo ;
  62. * la difference entre PROI parallele et sequentiel doit etre nulle
  63. errps = 'MAXI' ('-' chporp chpors) 'ABS' ;
  64. crit2 = ('VALE' 'PETI') '*' 10. ;
  65. 'MESS' 'errps=' errps ' crit2=' crit2 ;
  66. tst2 = errps '<' crit2 ;
  67. *
  68. lok = lok 'ET' tst2 ;
  69. *
  70. 'SAUT' 1 'LIGNE' ;
  71. 'SI' lok ;
  72.    'MESSAGE' 'Tout sest bien passe' ;
  73. 'SINON' ;
  74.    'MESSAGE' '!!! Il y a eu des erreurs' ;
  75. 'FINSI' ;   
  76. 'SAUT' 1 'LIGNE' ;
  77. *
  78. 'SI' interact ;
  79.    'OPTION' 'DONN' 5 'ECHO' 1 ;
  80. 'FINSI' ;   
  81. 'SI' ('NON' lok) ;
  82.    'ERREUR' 5 ;
  83. 'FINSI' ;   
  84. *
  85. * End of dgibi file PROI-PARALLELE
  86. *
  87. 'FIN' ;
  88.  
  89.  
  90.  
  91.  

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