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