Télécharger la plaquette Cast3M


fr en

Télécharger Fonction_Parallele.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : Fonction_Parallele.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5.  
  6. OPTI DIME 2 ELEM SEG2;
  7.  
  8. P1 = 0. 0. ;
  9. P2 = 1. 0. ;
  10. COMPLET = FAUX; 'COMM' 'Cas trop gros pour passer avec le "ulimit -m"';
  11. BCHPO   = VRAI; 'COMM' 'CHPOINT'            ;
  12. BLREE   = VRAI; 'COMM' 'LISTREEL'           ;
  13. BMCH1   = VRAI; 'COMM' 'MCHAML de REAL*8'   ;
  14. BMCH2   = VRAI; 'COMM' 'MCHAML de POINTEURS';
  15. BEVO1   = VRAI; 'COMM' 'EVOLUTIO'           ;
  16.  
  17.  
  18. SI BCHPO ;
  19.  
  20.   M1 = 'MANU' 'POI1' P1;
  21.  
  22.   CHP4 = 'MANU' 'CHPO' M1 4
  23.         'TREC' 0.5  'TRIC' 1. 'TRAC' 9.5  'TROC' 2.1  'NATURE' 'DIFFUS';
  24.   CHP6 = 'MANU' 'CHPO' M1 4
  25.         'TREC' 0.5  'TRIC' 1. 'TRAC' .21  'TROC' .67  'NATURE' 'DIFFUS';
  26.   CHP7 = 'MANU' 'CHPO' M1 4
  27.       'TREC' 1.0  'TRIC' 1.21 'TRAC' 2.67 'TROC' 4.85 'NATURE' 'DIFFUS';
  28.   CHP8 = 'MANU' 'CHPO' M1 4
  29.       'TREC' -.98  'TRIC' -.761 'TRAC' 0. 'TROC' .47  'NATURE' 'DIFFUS';
  30.  
  31. ****************
  32. * OPERATION '**'
  33. ****************
  34.   CPTEST =CHP4 ** 0    ;
  35.   CPTEST =CHP4 ** 1    ;
  36.   CPTEST =CHP4 ** 2    ;
  37.   CPTEST =CHP4 ** 2.   ;
  38.   CPTEST =CHP4 ** 0.5  ;
  39.   CPTEST =CHP4 ** 100.5;
  40.   CPTEST =(-1*CHP4) ** (1.D0/9.D0);
  41.   CPTEST =CHP4 ** 2.1  ;
  42.  
  43.  
  44.   CPTEST =2   ** CHP4;
  45.   CPTEST =2.  ** CHP4;
  46.   CPTEST =2.1 ** CHP4;
  47.   CPTEST =0.  ** CHP4;
  48.  
  49. ***************
  50. * OPERATION '*'
  51. ***************
  52.   CPTEST =CHP4 * 0    ;
  53.   CPTEST =CHP4 * 2    ;
  54.   CPTEST =CHP4 * 0.   ;
  55.   CPTEST =CHP4 * 2.1  ;
  56.  
  57. ***************
  58. * OPERATION '+'
  59. ***************
  60.   CPTEST =CHP4 + 0    ;
  61.   CPTEST =CHP4 + 1    ;
  62.   CPTEST =CHP4 + 0.   ;
  63.   CPTEST =CHP4 + 1.   ;
  64.  
  65. ***************
  66. * OPERATION '-'
  67. ***************
  68.   CPTEST =CHP4 - 0    ;
  69.   CPTEST =CHP4 - 1    ;
  70.   CPTEST =CHP4 - 0.   ;
  71.   CPTEST =CHP4 - 1.   ;
  72.  
  73.   CPTEST =0 - CHP4    ;
  74.   CPTEST =1 - CHP4    ;
  75.   CPTEST =0.- CHP4    ;
  76.   CPTEST =1.- CHP4    ;
  77.  
  78. ***************
  79. * OPERATION '/'
  80. ***************
  81.   CPTEST =CHP4 / 2    ;
  82.   CPTEST =CHP4 / 0.5  ;
  83.  
  84.   CPTEST =0 / CHP4    ;
  85.   CPTEST =0./ CHP4    ;
  86.   CPTEST =1 / CHP4    ;
  87.   CPTEST =1./ CHP4    ;
  88.  
  89. ***********
  90. * FONCTIONS
  91. ***********
  92.   CPTEST = COS  CHP4    ;
  93.   CPTEST = SIN  CHP4    ;
  94.   CPTEST = TAN  CHP4    ;
  95.   CPTEST = ACOS CHP6    ;
  96.   CPTEST = ASIN CHP6    ;
  97.   CPTEST = ATG  CHP4    ;
  98.   CPTEST = EXP  CHP4    ;
  99.   CPTEST = LOG  CHP4    ;
  100.   CPTEST = ABS  CHP4    ;
  101.   CPTEST = COSH CHP4    ;
  102.   CPTEST = SINH CHP4    ;
  103.   CPTEST = TANH CHP4    ;
  104.   CPTEST = ERF  CHP4    ;
  105.   CPTEST = ERFC CHP4    ;
  106.   CPTEST = ACOH CHP7    ;
  107.   CPTEST = ASIH CHP4    ;
  108.   CPTEST = ATAH CHP8    ;
  109.   CPTEST = SIGN CHP8    ;
  110.  
  111.  
  112.  
  113. * Test parallelisme
  114.   SI COMPLET;
  115.     NBSEG=100000000;
  116.     L2   = 'DROI' (NBSEG - 1) P1 P2 ;
  117.     CHP2 = 'MANU' 'CHPO' L2 2 'TROK' .25 'TRUC' .125 'NATURE' 'DIFFUS' ;
  119.     REPE SURI (10);
  120.       TOTO = CHP2 * 2.1354 ;
  121.       'DETR' TOTO;
  122.     FIN SURI;
  123.     TEMP 'IMPR' 'MAXI' 'HORL';
  124.   FINS;
  125.   FIN;
  126. FINS;
  127.  
  128. SI BLREE ;
  129.   LREE4 = 'PROG'  0.5    1.     9.5    2.1   ;
  130.   LREE6 = 'PROG'  0.5    1.     0.21   0.67  ;
  131.   LREE7 = 'PROG'  1.     1.21   2.67   4.85  ;
  132.   LREE8 = 'PROG'  -0.98 -0.761  0.     0.47  ;
  133.  
  134. ****************
  135. * OPERATION '**'
  136. ****************
  137.   LRTEST =LREE4 ** 0    ;
  138.   LRTEST =LREE4 ** 1    ;
  139.   LRTEST =LREE4 ** 2    ;
  140.   LRTEST =LREE4 ** 2.   ;
  141.   LRTEST =LREE4 ** 0.5  ;
  142.   LRTEST =LREE4 ** 100.5;
  143.   LRTEST =(-1*LREE4) ** (1.D0/9.D0);
  144.   LRTEST =LREE4 ** 2.1  ;
  145.  
  146.  
  147.   LRTEST =2   ** LREE4;
  148.   LRTEST =2.  ** LREE4;
  149.   LRTEST =2.1 ** LREE4;
  150.   LRTEST =0.  ** LREE4;
  151.  
  152. ***************
  153. * OPERATION '*'
  154. ***************
  155.   LRTEST =LREE4 * 0    ;
  156.   LRTEST =LREE4 * 2    ;
  157.   LRTEST =LREE4 * 0.   ;
  158.   LRTEST =LREE4 * 2.1  ;
  159.  
  160. ***************
  161. * OPERATION '+'
  162. ***************
  163.   LRTEST =LREE4 + 0    ;
  164.   LRTEST =LREE4 + 1    ;
  165.   LRTEST =LREE4 + 0.   ;
  166.   LRTEST =LREE4 + 1.   ;
  167.  
  168. ***************
  169. * OPERATION '-'
  170. ***************
  171.   LRTEST =LREE4 - 0    ;
  172.   LRTEST =LREE4 - 1    ;
  173.   LRTEST =LREE4 - 0.   ;
  174.   LRTEST =LREE4 - 1.   ;
  175.  
  176.   LRTEST =0 - LREE4    ;
  177.   LRTEST =1 - LREE4    ;
  178.   LRTEST =0.- LREE4    ;
  179.   LRTEST =1.- LREE4    ;
  180.  
  181. ***************
  182. * OPERATION '/'
  183. ***************
  184.   LRTEST =LREE4 / 2    ;
  185.   LRTEST =LREE4 / 0.5  ;
  186.  
  187.   LRTEST =0 / LREE4    ;
  188.   LRTEST =0./ LREE4    ;
  189.   LRTEST =1 / LREE4    ;
  190.   LRTEST =1./ LREE4    ;
  191.  
  192. ***********
  193. * FONCTIONS
  194. ***********
  195.   LRTEST = COS  LREE4    ;
  196.   LRTEST = SIN  LREE4    ;
  197.   LRTEST = TAN  LREE4    ;
  198.   LRTEST = ACOS LREE6    ;
  199.   LRTEST = ASIN LREE6    ;
  200.   LRTEST = ATG  LREE4    ;
  201.   LRTEST = EXP  LREE4    ;
  202.   LRTEST = LOG  LREE4    ;
  203.   LRTEST = ABS  LREE4    ;
  204.   LRTEST = COSH LREE4    ;
  205.   LRTEST = SINH LREE4    ;
  206.   LRTEST = TANH LREE4    ;
  207.   LRTEST = ERF  LREE4    ;
  208.   LRTEST = ERFC LREE4    ;
  209.   LRTEST = ACOH LREE7    ;
  210.   LRTEST = ASIH LREE4    ;
  211.   LRTEST = ATAH LREE8    ;  
  212.   LRTEST = SIGN LREE8    ;
  213.  
  214. * Test parallelisme
  215.   SI COMPLET;
  216.     NBVAL = 125000;
  217.     LREE1 = 'PROG' 1.  'PAS' (1. / (NBVAL-1)) 2.  ;
  219.     REPE SURI (30*(32**2));
  220.       TOTO = LREE1 * 2.1354 ;
  221.      'DETR' TOTO;
  222.     FIN SURI;
  223.     TEMP 'IMPR' 'MAXI' 'HORL';
  224.   FINS;
  225. FINS;
  226.  
  227. SI BMCH1 ;
  228.   NBSEG=2000000;
  229.  
  230.   L1 = 'DROI'    NBSEG  P1 P2 ;
  231.   L2 = 'DROI' (2*NBSEG) P1 P2 ;
  232.   MOD1  = 'MODE' L1 'MECANIQUE' 'ELASTIQUE';
  233.   MOD2  = 'MODE' L2 'MECANIQUE' 'ELASTIQUE';
  234.  
  235.   COOR1 = 'COOR' 1 L1 ;
  236.   COOR2 = 'COOR' 1 L2 ;
  237.  
  238. *----------------------------------------------------------------------*
  239. * FABRICATION DE MCHAML CONSTANT, VARIABLE ET CONSTANT PAR SOUS-ZONE
  240. *----------------------------------------------------------------------*
  241. *  MCHAML CONSTANT PAR SOUS-ZONE
  242.   CHAM1  = 'MANU' 'CHML' MOD1 'TREC' 1.   'STRESSES';
  243.   CHAM2  = 'MANU' 'CHML' MOD2 'TROK' .25  'STRESSES';
  244.  
  245. *  MCHAML VARIABLE
  246.   CHAM3 = 'CHAN' 'COMP' ('CHAN' 'CHAM' COOR1 MOD1 'STRESSES') 'TREC';
  247.   CHAM4 = 'CHAN' 'COMP' ('CHAN' 'CHAM' COOR2 MOD2 'STRESSES') 'TROK';
  248.  
  249. * MCHAML CONSTANT PAR ELEMENT
  250.   CHAM5 = 'CHAN' 'COMP' ('CHAN' 'GRAVITE' CHAM3 MOD1) 'TREC';
  251.   CHAM6 = 'CHAN' 'COMP' ('CHAN' 'GRAVITE' CHAM4 MOD2) 'TROK';
  252.   CHAM5 = 'CHAN' 'STRESSES' CHAM5 MOD1;
  253.   CHAM6 = 'CHAN' 'STRESSES' CHAM6 MOD2;
  254.  
  255. *----------------------------------------------------------------------*
  256. * OPERATION A 2 ARGUMENTS
  257. *----------------------------------------------------------------------*
  258. *  'LIST' 'RESU' (ATG CHAM1 (CHAM1 * 2.)) ; COMM 'CST    CST'   ;
  259. *  'LIST' 'RESU' (ATG CHAM1 CHAM5       ) ; COMM 'CST    CST_EL';
  260. *  'LIST' 'RESU' (ATG CHAM5 CHAM1       ) ; COMM 'CST_EL CST'   ;
  261. *  'LIST' 'RESU' (ATG CHAM1 CHAM3       ) ; COMM 'CST    VARIAB';
  262. *  'LIST' 'RESU' (ATG CHAM3 CHAM1       ) ; COMM 'VARIAB CST'   ;
  263. *  'LIST' 'RESU' (ATG CHAM5 CHAM3       ) ; COMM 'CST_EL VARIAB';
  264. *  'LIST' 'RESU' (ATG CHAM3 CHAM5       ) ; COMM 'VARIAB CST_EL';
  265.  
  267.   REPE SURI 30;
  268.     TOTO = CHAM4 ** 2.5 ;
  269.     'DETR' TOTO;
  270.   FIN SURI;
  271.   TEMP 'IMPR' 'MAXI' 'HORL';
  272.  
  273. * Cohabitation pthread / ASSISTANTS
  274.   MAI1 = EXTR CHAM4 'MAIL';
  275.   OPTI 'PARA' VRAI;
  276.   MOD3 = PART 4 MOD2;
  277.   CHAM7 = REDU CHAM4 MOD3;
  278.  
  279. * Test parallelisme
  280.   SI COMPLET;
  282.   REPE SURI 30;
  283.     TOTO = CHAM7 ** 2.5 ;
  284.     'DETR' TOTO;
  285.   FIN SURI;
  286.   TEMP 'IMPR' 'MAXI' 'HORL';
  287.   FINS;  
  288. FINS;
  289.  
  290. SI BMCH2 ;
  291.   NBSEG=1;
  292.  
  293.   L1 = 'DROI'    NBSEG  P1 P2 ;
  294.   L2 = 'DROI' (2*NBSEG) P1 P2 ;
  295.   MOD1  = 'MODE' L1 'MECANIQUE' 'ELASTIQUE';
  296.   MOD2  = 'MODE' L2 'MECANIQUE' 'ELASTIQUE';
  297.  
  298.   COOR1 = 'COOR' 1 L1 ;
  299.   COOR2 = 'COOR' 1 L2 ;
  300.  
  301.   NBVAL = 10000000;
  302.   LREE1 = 'PROG' 1. 'PAS'  (1. / NBVAL) 2.  ;
  303.   EVOL1 = 'EVOL' 'MANU' LREE1 LREE1;
  304.   CHAM1LR= 'MANU' 'CHML' MOD1 'TROK' LREE1  'STRESSES';
  305.   CHAM1EV= 'MANU' 'CHML' MOD1 'TROK' EVOL1  'STRESSES';
  306.  
  307.   LENT1 = 'LECT' 1 'PAS'  1 NBVAL  ;
  308.   EVOL2 = 'EVOL' 'MANU' LENT1 LENT1;
  309.   CHAM1LE= 'MANU' 'CHML' MOD1 'TROK' LENT1  'STRESSES';
  310.   CHAM2EV= 'MANU' 'CHML' MOD1 'TROK' EVOL2  'STRESSES';
  311. *  LIST RESU CHAM1LR;
  312. *  LIST RESU (CHAM1LR ** 2);
  313.  
  314. * Test parallelisme
  315.   SI COMPLET;
  317.   REPE SURI 50;
  318.     TOTO = 'ASIH' CHAM1EV ;
  319.     'DETR' TOTO;
  320.   FIN SURI;
  321.   TEMP 'IMPR' 'MAXI' 'HORL';
  322.   FINS;
  323. FINS;
  324.  
  325. SI BEVO1;
  326.   LMOT1 = 'MOTS' 'TOTO' 'TITI' 'TATA' 'TRUC' ;
  327.   LENT2 = 'LECT'  1      2      3      4     ;
  328.   LREE3 = 'PROG'  0.5    1.     9.5    2.1   ;
  329.   LREE4 = 'PROG'  0.5    1.     0.21   0.67  ;
  330.   LREE5 = 'PROG'  1.     1.21   2.67   4.85  ;
  331.   LREE6 = 'PROG'  -0.98 -0.761  0.     0.47  ;
  332.   EVO3  = 'EVOL' 'MANU' 'XX' LMOT1 'IY' LENT2;
  333.   EVO4  = 'EVOL' 'MANU' 'XX' LMOT1 'IY' LREE3;
  334.   EVO5  = 'EVOL' 'MANU' 'XX' LENT2 'IY' LMOT1;
  335.   EVO6  = 'EVOL' 'MANU' 'XX' LMOT1 'IY' LREE4;
  336.   EVO7  = 'EVOL' 'MANU' 'XX' LMOT1 'IY' LREE5;
  337.   EVO8  = 'EVOL' 'MANU' 'XX' LMOT1 'IY' LREE6;
  338.  
  339. ****************
  340. * OPERATION '**'
  341. ****************
  342.   EVTEST =EVO4 ** 0    ;
  343.   EVTEST =EVO4 ** 1    ;
  344.   EVTEST =EVO4 ** 2    ;
  345.   EVTEST =EVO4 ** 2.   ;
  346.   EVTEST =EVO4 ** 0.5  ;
  347.   EVTEST =EVO4 ** 100.5;
  348.   EVTEST =(-1*EVO4) ** (1.D0/9.D0);
  349.   EVTEST =EVO4 ** 2.1  ;
  350.  
  351.   EVTEST =2   ** EVO4;
  352.   EVTEST =2.  ** EVO4;
  353.   EVTEST =2.1 ** EVO4;
  354.   EVTEST =0.  ** EVO4;
  355.  
  356. ***************
  357. * OPERATION '*'
  358. ***************
  359.   EVTEST =EVO4 * 0    ;
  360.   EVTEST =EVO4 * 2    ;
  361.   EVTEST =EVO4 * 0.   ;
  362.   EVTEST =EVO4 * 2.1  ;
  363.  
  364. ***************
  365. * OPERATION '+'
  366. ***************
  367.   EVTEST =EVO4 + 0    ;
  368.   EVTEST =EVO4 + 1    ;
  369.   EVTEST =EVO4 + 0.   ;
  370.   EVTEST =EVO4 + 1.   ;
  371.  
  372. ***************
  373. * OPERATION '-'
  374. ***************
  375.   EVTEST =EVO4 - 0    ;
  376.   EVTEST =EVO4 - 1    ;
  377.   EVTEST =EVO4 - 0.   ;
  378.   EVTEST =EVO4 - 1.   ;
  379.  
  380.   EVTEST =0 - EVO4    ;
  381.   EVTEST =1 - EVO4    ;
  382.   EVTEST =0.- EVO4    ;
  383.   EVTEST =1.- EVO4    ;
  384.  
  385. ***************
  386. * OPERATION '/'
  387. ***************
  388.   EVTEST =EVO4 / 2    ;
  389.   EVTEST =EVO4 / 0.5  ;
  390.  
  391.   EVTEST =0 / EVO4    ;
  392.   EVTEST =0./ EVO4    ;
  393.   EVTEST =1 / EVO4    ;
  394.   EVTEST =1./ EVO4    ;
  395.  
  396. ***********
  397. * FONCTIONS
  398. ***********
  399.   EVTEST = COS  EVO4    ;
  400.   EVTEST = SIN  EVO4    ;
  401.   EVTEST = TAN  EVO4    ;
  402.   EVTEST = ACOS EVO6    ;
  403.   EVTEST = ASIN EVO6    ;
  404.   EVTEST = ATG  EVO4    ;
  405.   EVTEST = EXP  EVO4    ;
  406.   EVTEST = LOG  EVO4    ;
  407.   EVTEST = ABS  EVO4    ;
  408.   EVTEST = COSH EVO4    ;
  409.   EVTEST = SINH EVO4    ;
  410.   EVTEST = TANH EVO4    ;
  411.   EVTEST = ERF  EVO4    ;
  412.   EVTEST = ERFC EVO4    ;
  413.   EVTEST = ACOH EVO7    ;
  414.   EVTEST = ASIH EVO4    ;
  415.   EVTEST = ATAH EVO8    ; 
  416.   EVTEST = SIGN EVO8    ;
  417.  
  418. * Test parallelisme
  419.   SI COMPLET;
  420.     NBVAL = 101280;
  421.     LREE1 = 'PROG' 1. 'PAS'  (1. / (NBVAL-1)) 2.  ;
  422.     LREE2 = LREE1 + 1 ;
  423.     EVO1  = 'EVOL' 'MANU' 'XX' LREE1 'YY' LREE2;
  424.  
  425.     LENT1 = 'ENTI' (100 * LREE1);
  426.     EVO2  = 'EVOL' 'MANU' 'XX' LREE1 'IY' LENT1;
  427.  
  429.     REPE SURI (10*(32**2)*10);
  430.       TOTO = EVO1  * 2.156 ;
  431.      'DETR' TOTO;
  432.     FIN SURI;
  433.     TEMP 'IMPR' 'MAXI' 'HORL';
  434.   FINS;
  435. FINS;
  436.  
  437.  
  438. FIN;
  439.  
  440.  
  441.  
  442.  
  443.  

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