Télécharger la plaquette Cast3M


fr en

Télécharger konv_impl_murs.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  konv_impl_murs.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***********************************************************
  5. ***********************************************************
  6. **** APPROCHE VF "Cell-Centred Formulation" pour la    ****
  7. **** solution des                                      **** 
  8. **** Equations d'Euler pour un gaz parfait             ****
  9. **** OPERATEURS PRIM, PRET, KONV                       ****
  10. **** Implicit: calcul du jacobien du residu            ****
  11. ****                                                   ****
  12. **** Cas gaz monoespece, "calorically perfect"         ****
  13. ****                                                   ****
  14. **** Methodes: VLH                                     ****
  15. ****                                                   ****
  16. **** A. BECCANTINI DRN/DMT/SEMT/LTMF    AOUT  2000     ****
  17. ***********************************************************
  18. ***********************************************************
  19.  
  20.  'OPTION'  'DIME' 2 ;
  21.  'OPTION'  'ELEM' QUA4 ;
  22.  'OPTION'  'ECHO'  0 ;
  23.  'OPTION'  'TRAC' 'X' ;
  24.  
  25. *
  26. *** GRAPH
  27. *
  28.  
  29.  GRAPH = FAUX ;
  30. *  GRAPH = VRAI ;
  31.  
  32.  ERRTOL = 1.0D-3 ;
  33.  
  34. ***************************
  35. ***** DOMAINE SPATIAL  ****
  36. ***************************
  37.  
  38.  
  39.  A1 = 0.0D0 0.0D0;
  40.  A2 = 1.0D0 0.0D0;
  41.  A3 = 1.0D0 1.0D0;
  42.  A4 = 0.0D0 1.0D0;
  43.  
  44.  L12 = A1 'DROIT' 1 A2;
  45.  L23 = A2 'DROIT' 1 A3;
  46.  L34 = A3 'DROIT' 1 A4;
  47.  L41 = A4 'DROIT' 1 A1;
  48.  
  49.  DOM10 = 'DALL' L12 L23 L34 L41
  50.         'PLANE';
  51.  
  52. *
  53. *** Etat 
  54. *
  55.  
  56.  ro  = 1.11 ;
  57.  p   = 1234.7;
  58.  gam = 1.4D0;
  59.  cson = (gam '*' p '/' ro) '**' 0.5 ;
  60.  un  = 0.2 * cson ;
  61.  ut  = 0.3  * cson ;
  62.  
  63. *
  64. *** retgd
  65. *
  66.  
  67.  ecin = 0.5D0 '*' ro '*' ((un  '*' un ) '+' (ut  '*' ut ));
  68.  ret  = (p  '/' (gam  '-' 1.0)) '+' ecin ;
  69.  
  70.  
  71.  
  72. *
  73. ***  flux en (n,t)
  74. *
  75.  
  76. ****************************************************
  77. ****************************************************
  78. ******** Boucle sur les angles             *********
  79. ****************************************************
  80. ****************************************************
  81.  
  82.  DANGLE =  (360.0 '/' 7.15) ;
  83.  ANGLE = 0.0 ;
  84.  
  85.  'REPETER' BLOC  8 ;
  86.  
  87. *
  88. *** Rotation
  89. *
  90.  
  91.     ANGLE = ANGLE '+' DANGLE ;
  92.     ORIG = 0.0D0 0.0D0;
  93.  
  94.     'MESSAGE' ;
  95.     'MESSAGE' ('CHAINE'  'Angle de rotation= ' ANGLE);
  96.     'MESSAGE' ;
  97.  
  98.     DOMTOT = DOM10   'TOURNER' ANGLE ORIG;
  99.     $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  100.     TDOMTOT = 'DOMA'  $DOMTOT 'VF';
  101.     MDOMTOT = TDOMTOT . 'QUAF' ;
  102.     'ELIM' MDOMTOT  1.E-6 ;
  103.  
  104.     PCON = ('DOMA' $DOMTOT 'CENTRE') 'POIN' 1 ;
  105.  
  106.     'SI' GRAPH;
  107.       'TRACER' (('DOMA' $DOMTOT 'MAILLAGE' ) 'ET' 
  108.          ('DOMA' $DOMTOT 'CENTRE')) 'TITRE' 'Domaine et centre' ;
  109.     'FINSI' ;
  110.  
  111. ***********************
  112. **** Les CHPOINTs  ****
  113. ***********************
  114.  
  115.     ux  = (un  '*' ('COS' ANGLE)) '-' (ut  '*' ('SIN' ANGLE));
  116.     uy  = (un  '*' ('SIN' ANGLE)) '+' (ut  '*' ('COS' ANGLE));
  117.  
  118.  
  119.     ro0 = ro ;
  120.     gnx0 = ro '*' ux ;
  121.     gny0 = ro '*' uy ;
  122.     ret0 = ret ;
  123.  
  124.     GAMMAN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' gam ;
  125.  
  126.     RN  = 'KCHT' $DOMTOT   'SCAL' 'CENTRE' ro0 ;
  127.  
  128.     GN  = 'KCHT' $DOMTOT   'VECT' 'CENTRE' (gnx0 gny0) ;
  129.  
  130.     RETN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' ret0 ;
  131.  
  132. ***************************
  133. ****  L'operateur  PRIM****
  134. ***************************
  135.  
  136.     VITESSE PRES   = 'PRIM' 'PERFMONO'  RN GN RETN GAMMAN ;
  137.  
  138. ***************************
  139. ****  L'operateur  PRET****
  140. ***************************
  141.  
  142.     ORDESP = 1;
  143.     ORDTEM = 1;
  144.  
  145.     ROF VITF PF GAMF   =  'PRET' 'PERFMONO'  ORDESP ORDTEM
  146.                            $DOMTOT RN VITESSE PRES GAMMAN ;
  147.  
  148. ***************************
  149. ****  L'operateur  KONV****
  150. ***************************
  151.  
  152.     METO = 'VLH' ;
  153.  
  154.     LISTINCO = 'MOTS' 'RN' 'RUXN' 'RUYN' 'RETN' ;
  155.  
  156.     CHPRES DT = 'KONV' 'VF' 'PERFMONO' 'RESI' METO
  157.           $DOMTOT ROF VITF PF GAMF LISTINCO ;
  158.  
  159.     IJACO = 'KONV' 'VF' 'PERFMONO' 'JACOCONS'
  160.           $DOMTOT LISTINCO
  161.           METO RN VITESSE PRES GAMMAN ;
  162.  
  163. ***********************
  164. ***** Le residu *******
  165. ***********************
  166.  
  167.  
  168.    DEBRN0   = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 1) ;
  169.    DEBGNX0  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 2) ;
  170.    DEBGNY0  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 3) ;
  171.    DEBRETN0 = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 4) ;
  172.  
  173. *
  174. ***** Le jacobien du residu et le residu
  175. *
  176.  
  177.  
  178.    UNRO = 'MANUEL' 'CHPO' ('DOMA' $DOMTOT 'CENTRE') 4 'RN' 1.0
  179.       'RUXN' 0.0 'RUYN' 0.0 'RETN' 0.0
  180.         'NATURE' 'DISCRET' ;
  181.    UNGX = 'MANUEL' 'CHPO' ('DOMA' $DOMTOT 'CENTRE') 4 'RN' 0.0
  182.       'RUXN' 1.0 'RUYN' 0.0 'RETN' 0.0
  183.         'NATURE' 'DISCRET' ;
  184.    UNGY = 'MANUEL' 'CHPO' ('DOMA' $DOMTOT 'CENTRE') 4 'RN' 0.0
  185.       'RUXN' 0.0 'RUYN' 1.0 'RETN' 0.0
  186.         'NATURE' 'DISCRET' ;
  187.    UNRET = 'MANUEL' 'CHPO' ('DOMA' $DOMTOT 'CENTRE') 4 'RN' 0.0
  188.       'RUXN' 0.0 'RUYN' 0.0 'RETN' 1.0
  189.         'NATURE' 'DISCRET' ;
  190.  
  191.  
  192.    DDRHO = 'KOPS' IJACO 'MULT' UNRO ;
  193.  
  194. *   
  195. *
  196. ***** DDRHO contient: dRES_RN '/' dRN ('RN') ;
  197. *                     dRES_GXN '/' dRN ('RUXN') ; 
  198. *                     dRES_GYN '/' dRN ('RUYN') ;
  199. *                     dRES_RETN '/' dRN ('RETN') ;
  200. *    DDGX contient ...
  201. *
  202. *
  203.  
  204.    DDGX = 'KOPS' IJACO 'MULT' UNGX ;
  205.    DDGY = 'KOPS' IJACO 'MULT' UNGY ;    
  206.    DDRET = 'KOPS' IJACO 'MULT' UNRET ;
  207.  
  208.    DRR  = 'EXTRAIRE' DDRHO PCON  'RN' ;
  209.    DGXR = 'EXTRAIRE' DDRHO PCON 'RUXN' ;
  210.    DGYR = 'EXTRAIRE' DDRHO PCON 'RUYN' ;
  211.    DRETR = 'EXTRAIRE' DDRHO PCON 'RETN' ;
  212.  
  213.    DRGX   = 'EXTRAIRE' DDGX PCON  'RN' ;
  214.    DGXGX  = 'EXTRAIRE' DDGX PCON 'RUXN' ;
  215.    DGYGX  = 'EXTRAIRE' DDGX PCON 'RUYN' ;
  216.    DRETGX = 'EXTRAIRE' DDGX PCON 'RETN' ;
  217.  
  218.    DRGY   = 'EXTRAIRE' DDGY PCON  'RN' ;
  219.    DGXGY = 'EXTRAIRE' DDGY PCON 'RUXN' ;
  220.    DGYGY = 'EXTRAIRE' DDGY PCON 'RUYN' ;
  221.    DRETGY = 'EXTRAIRE' DDGY PCON 'RETN' ;
  222.  
  223.    DRRET   = 'EXTRAIRE' DDRET PCON  'RN' ;
  224.    DGXRET = 'EXTRAIRE' DDRET PCON 'RUXN' ;
  225.    DGYRET = 'EXTRAIRE' DDRET PCON 'RUYN' ;
  226.    DRETRET = 'EXTRAIRE' DDRET PCON 'RETN' ;
  227.  
  228. *********************************************************************
  229. *********************************************************************
  230. ***** On calcule les residues pour ro1 = ro0 * (1.'+' DELTA)  *******
  231. *********************************************************************
  232. *********************************************************************
  233.  
  234.     DELTA = 1.0D-4 ; 
  235.     ro1 = ro0 '*' (1.0 '+' DELTA) ;
  236.     gnx1 = gnx0 ;
  237.     gny1 = gny0 ;
  238.     ret1 = ret0 ;
  239.  
  240.  
  241.     GAMMAN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' gam ;
  242.  
  243.     RN  = 'KCHT' $DOMTOT   'SCAL' 'CENTRE' ro1 ;
  244.  
  245.     GN  = 'KCHT' $DOMTOT   'VECT' 'CENTRE' (gnx1 gny1) ;
  246.  
  247.     RETN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' ret1 ;
  248.  
  249.     VITESSE PRES   = 'PRIM' 'PERFMONO'  RN GN RETN GAMMAN ;
  250.  
  251.  
  252.     ROF VITF PF GAMF   =  'PRET' 'PERFMONO'  ORDESP ORDTEM
  253.                            $DOMTOT RN VITESSE PRES GAMMAN ;
  254.  
  255.     CHPRES DT = 'KONV' 'VF' 'PERFMONO' 'RESI' METO
  256.           $DOMTOT ROF VITF PF GAMF LISTINCO ;
  257.  
  258. ***********************
  259. ***** Le residu *******
  260. ***********************
  261.  
  262.    DEBRN1   = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 1) ;
  263.    DEBGNX1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 2) ;
  264.    DEBGNY1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 3) ;
  265.    DEBRETN1 = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 4) ;
  266.  
  267.  
  268. *
  269. **** On calcule le jacobien numeriquement
  270. *
  271.  
  272.    DRRN = (DEBRN1 '-' DEBRN0) '/' (DELTA '*' ro0) ;
  273.    DGXRN = (DEBGNX1 '-' DEBGNX0) '/' (DELTA '*' ro0) ;
  274.    DGYRN = (DEBGNY1 '-' DEBGNY0) '/' (DELTA '*' ro0) ;
  275.    DRETRN = (DEBRETN1 '-' DEBRETN0) '/' (DELTA '*' ro0) ;
  276.  
  277.  
  278.    ERR1 = ('ABS' (DRRN '-' DRR)) '*' (ro0 '/' (ro0 '*' cson)) ;
  279.  
  280.    'SI' (ERR1 > ERRTOL) ;
  281.         'MESSAGE' 'Problem 1';
  282.         'ERREUR' 5 ;
  283.    'FINSI' ;
  284.  
  285.    ERR1 = ('ABS' (DGXRN '-' DGXR)) '*'
  286.           (ro0 '/' (ro0 '*' cson '*' cson)) ;
  287.  
  288.    'SI' (ERR1 > ERRTOL) ;
  289.         'MESSAGE' 'Problem 2';
  290.         'ERREUR' 5 ;
  291.    'FINSI' ;
  292.  
  293.    ERR1 = ('ABS' (DGYRN '-' DGYR)) '*'
  294.           (ro0 '/' (ro0 '*' cson '*' cson)) ;
  295.  
  296.    'SI' (ERR1 > ERRTOL) ;
  297.         'MESSAGE' 'Problem 3';
  298.         'ERREUR' 5 ;
  299.    'FINSI' ;
  300.  
  301.  
  302.    ERR1 = ('ABS' (DRETRN '-' DRETR)) '*'
  303.           (ro0 '/' (ret0 '*' cson)) ; 
  304.  
  305.    'SI' (ERR1 > ERRTOL) ;
  306.         'MESSAGE' 'Problem 4';
  307.         'ERREUR' 5 ;
  308.    'FINSI' ;
  309.  
  310.  
  311.  
  312. ***************************************************************************
  313. ***************************************************************************
  314. ***** On calcule les residues pour gnx1 = gnx0 '+' (DELTA ro0 cson)  ******
  315. ***************************************************************************
  316. ***************************************************************************
  317.  
  318.  
  319.     DELTA = 1.0D-4 ;
  320.  
  321.     ro1 = ro0 ; 
  322.     gnx1 = gnx0 '+' (DELTA '*' ro0 '*' cson)  ;
  323.     gny1 = gny0 ;
  324.     ret1 = ret0 ;
  325.  
  326.  
  327.     GAMMAN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' gam ;
  328.  
  329.     RN  = 'KCHT' $DOMTOT   'SCAL' 'CENTRE' ro1 ;
  330.  
  331.     GN  = 'KCHT' $DOMTOT   'VECT' 'CENTRE' (gnx1 gny1) ;
  332.  
  333.     RETN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' ret1 ;
  334.  
  335.     VITESSE PRES   = 'PRIM' 'PERFMONO'  RN GN RETN GAMMAN ;
  336.  
  337.  
  338.     ROF VITF PF GAMF   =  'PRET' 'PERFMONO'  ORDESP ORDTEM
  339.                            $DOMTOT RN VITESSE PRES GAMMAN ;
  340.  
  341.     CHPRES DT = 'KONV' 'VF' 'PERFMONO' 'RESI' METO
  342.           $DOMTOT ROF VITF PF GAMF LISTINCO ;
  343.  
  344. ***********************
  345. ***** Le residu *******
  346. ***********************
  347.  
  348.    DEBRN1   = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 1) ;
  349.    DEBGNX1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 2) ;
  350.    DEBGNY1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 3) ;
  351.    DEBRETN1 = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 4) ;
  352.  
  353. *
  354. **** On calcule le jacobien numeriquement
  355. *
  356.  
  357.    DRGXN = (DEBRN1 '-' DEBRN0) '/' (DELTA '*' ro0 '*' cson) ;
  358.    DGXGXN = (DEBGNX1 '-' DEBGNX0) '/' (DELTA '*' ro0 '*' cson) ;
  359.    DGYGXN = (DEBGNY1 '-' DEBGNY0) '/'  (DELTA '*' ro0 '*' cson) ;
  360.    DRETGXN = (DEBRETN1 '-' DEBRETN0) '/'  (DELTA '*' ro0 '*' cson) ;
  361.  
  362.  
  363.    ERR1 = ('ABS' (DRGXN '-' DRGX))  ;
  364.  
  365.    'SI' (ERR1 > ERRTOL) ;
  366.         'MESSAGE' 'Problem 5';
  367.         'ERREUR' 5 ;
  368.    'FINSI' ;
  369.  
  370.    ERR1 = ('ABS' (DGXGXN '-' DGXGX)) '/' cson ;
  371.  
  372.    'SI' (ERR1 > ERRTOL) ;
  373.         'MESSAGE' 'Problem 6';
  374.         'ERREUR' 5 ;
  375.    'FINSI' ;
  376.  
  377.    ERR1 = ('ABS' (DGYGXN '-' DGYGX)) '/' cson ;
  378.  
  379.    'SI' (ERR1 > ERRTOL) ;
  380.         'MESSAGE' 'Problem 7';
  381.         'ERREUR' 5 ;
  382.    'FINSI' ;
  383.  
  384.  
  385.    ERR1 = ('ABS' (DRETGXN '-' DRETGX)) '*'
  386.         (ro0  / ret0) ;
  387.  
  388.    'SI' (ERR1 > ERRTOL) ;
  389.         'MESSAGE' 'Problem 8';
  390.         'ERREUR' 5 ;
  391.    'FINSI' ;
  392.  
  393.  
  394.  
  395.  
  396. ***************************************************************************
  397. ***************************************************************************
  398. ***** On calcule les residues pour gny1 = gny0 '+' (DELTA ro0 cson)  ******
  399. ***************************************************************************
  400. ***************************************************************************
  401.  
  402.  
  403.     DELTA = 1.0D-4 ;
  404.  
  405.     ro1 = ro0 ; 
  406.     gnx1 = gnx0 ; 
  407.     gny1 = gny0 '+' (DELTA '*' ro0 '*' cson) ;
  408.     ret1 = ret0 ;
  409.  
  410.  
  411.     GAMMAN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' gam ;
  412.  
  413.     RN  = 'KCHT' $DOMTOT   'SCAL' 'CENTRE' ro1 ;
  414.  
  415.     GN  = 'KCHT' $DOMTOT   'VECT' 'CENTRE' (gnx1 gny1) ;
  416.  
  417.     RETN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' ret1 ;
  418.  
  419.     VITESSE PRES   = 'PRIM' 'PERFMONO'  RN GN RETN GAMMAN ;
  420.  
  421.  
  422.     ROF VITF PF GAMF   =  'PRET' 'PERFMONO'  ORDESP ORDTEM
  423.                            $DOMTOT RN VITESSE PRES GAMMAN ;
  424.  
  425.     CHPRES DT = 'KONV' 'VF' 'PERFMONO' 'RESI' METO
  426.           $DOMTOT ROF VITF PF GAMF LISTINCO ;
  427.  
  428. ***********************
  429. ***** Le residu *******
  430. ***********************
  431.  
  432.    DEBRN1   = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 1) ;
  433.    DEBGNX1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 2) ;
  434.    DEBGNY1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 3) ;
  435.    DEBRETN1 = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 4) ;
  436.  
  437. *
  438. **** On calcule le jacobien numeriquement
  439. *
  440.  
  441.    DRGYN = (DEBRN1 '-' DEBRN0) '/' (DELTA '*' ro0 '*' cson) ;
  442.    DGXGYN = (DEBGNX1 '-' DEBGNX0) '/' (DELTA '*' ro0 '*' cson) ;
  443.    DGYGYN = (DEBGNY1 '-' DEBGNY0) '/'  (DELTA '*' ro0 '*' cson) ;
  444.    DRETGYN = (DEBRETN1 '-' DEBRETN0) '/'  (DELTA '*' ro0 '*' cson) ;
  445.  
  446.    ERR1 = ('ABS' (DRGYN '-' DRGY))  ;
  447.  
  448.    'SI' (ERR1 > ERRTOL) ;
  449.         'MESSAGE' 'Problem 9';
  450.         'ERREUR' 5 ;
  451.    'FINSI' ;
  452.  
  453.    ERR1 = ('ABS' (DGXGYN '-' DGXGY))  '/' cson ;
  454.  
  455.    'SI' (ERR1 > ERRTOL) ;
  456.         'MESSAGE' 'Problem 10';
  457.         'ERREUR' 5 ;
  458.    'FINSI' ;
  459.  
  460.    ERR1 = ('ABS' (DGYGYN '-' DGYGY)) '/' cson ;
  461.  
  462.    'SI' (ERR1 > ERRTOL) ;
  463.         'MESSAGE' 'Problem 11';
  464.         'ERREUR' 5 ;
  465.    'FINSI' ;
  466.  
  467.  
  468.    ERR1 = ('ABS' (DRETGYN '-' DRETGY))  '*'
  469.         (ro0  / ret0)  ;
  470.  
  471.    'SI' (ERR1 > ERRTOL) ;
  472.         'MESSAGE' 'Problem 12';
  473.         'ERREUR' 5 ;
  474.    'FINSI' ;
  475.  
  476.  
  477. ***************************************************************************
  478. ***************************************************************************
  479. ***** On calcule les residues pour ret1 = ret0 '*' (1. '+' DELTA)    ******
  480. ***************************************************************************
  481. ***************************************************************************
  482.  
  483.  
  484.     DELTA = 1.0D-4 ;
  485.  
  486.     ro1 = ro0 ; 
  487.     gnx1 = gnx0 ; 
  488.     gny1 = gny0 ; 
  489.     ret1 = ret0 * (1. '+' DELTA) ;
  490.  
  491.  
  492.     GAMMAN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' gam ;
  493.  
  494.     RN  = 'KCHT' $DOMTOT   'SCAL' 'CENTRE' ro1 ;
  495.  
  496.     GN  = 'KCHT' $DOMTOT   'VECT' 'CENTRE' (gnx1 gny1) ;
  497.  
  498.     RETN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' ret1 ;
  499.  
  500.     VITESSE PRES   = 'PRIM' 'PERFMONO'  RN GN RETN GAMMAN ;
  501.  
  502.  
  503.     ROF VITF PF GAMF   =  'PRET' 'PERFMONO'  ORDESP ORDTEM
  504.                            $DOMTOT RN VITESSE PRES GAMMAN ;
  505.  
  506.     CHPRES DT = 'KONV' 'VF' 'PERFMONO' 'RESI' METO
  507.           $DOMTOT ROF VITF PF GAMF LISTINCO ;
  508.  
  509. ***********************
  510. ***** Le residu *******
  511. ***********************
  512.  
  513.    DEBRN1   = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 1) ;
  514.    DEBGNX1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 2) ;
  515.    DEBGNY1  = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 3) ;
  516.    DEBRETN1 = 'EXTRAIRE' CHPRES PCON ('EXTRAIRE' LISTINCO 4) ;
  517.  
  518. *
  519. **** On calcule le jacobien numeriquement
  520. *
  521.  
  522.    DRRETN = (DEBRN1 '-' DEBRN0) '/' (DELTA '*' ret0) ;
  523.    DGXRETN = (DEBGNX1 '-' DEBGNX0) '/' (DELTA '*' ret0) ;
  524.    DGYRETN = (DEBGNY1 '-' DEBGNY0) '/'  (DELTA '*' ret0) ;
  525.    DRETRETN = (DEBRETN1 '-' DEBRETN0) '/'  (DELTA '*' ret0) ;
  526.  
  527.    ERR1 = ('ABS' (DRRETN '-' DRRET)) '*' (ret0 '/' (ro0 '*' cson))   ;
  528.  
  529.    'SI' (ERR1 > ERRTOL) ;
  530.         'MESSAGE' 'Problem 13';
  531.         'ERREUR' 5 ;
  532.    'FINSI' ;
  533.  
  534.    ERR1 = ('ABS' (DGXRETN '-' DGXRET)) '*'
  535.           (ret0 '/' (ro0 '*' cson * cson))   ;
  536.  
  537.    'SI' (ERR1 > ERRTOL) ;
  538.         'MESSAGE' 'Problem 14';
  539.         'ERREUR' 5 ;
  540.    'FINSI' ;
  541.  
  542.    ERR1 = ('ABS' (DGYRETN '-' DGYRET)) '*' 
  543.              (ret0 '/' (ro0 '*' cson * cson)) ;
  544.  
  545.    'SI' (ERR1 > ERRTOL) ;
  546.         'MESSAGE' 'Problem 15';
  547.         'ERREUR' 5 ;
  548.    'FINSI' ;
  549.  
  550.  
  551.    ERR1 = ('ABS' (DRETRETN '-' DRETRET)) '/' cson ;
  552.  
  553.    'SI' (ERR1 > ERRTOL) ;
  554.         'MESSAGE' 'Problem 16';
  555.         'ERREUR' 5 ;
  556.    'FINSI' ;
  557.  
  558.  
  559.  'FIN' BLOC ;
  560.  
  561.  'FIN' ;
  562.  
  563.  
  564.  
  565.  
  566.  
  567.  
  568.  
  569.  
  570.  

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