1 : C XDISP SOURCE MAGN 09/07/30 21:21:15 6458 2 : SUBROUTINE XDISP(MATRIK,INCX,MCHSOL,IMPR,IRET) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : C*********************************************************************** 6 : C NOM : XDISP 7 : C DESCRIPTION : 8 : C 9 : C Ce sp transforme le vecteur des inconnues primales 10 : C (i.e. x de Ax=b) en un chpoint qui N'EST PAS au standard 11 : C CASTEM 2000 (en vue de l'exploitation par l'opérateur 12 : C GIBIANE 'TCRR') 13 : C 14 : C Ce sp reprend la structure de mexini.eso 15 : C C'est une interface à vec2ch. 16 : C 17 : C 18 : C LANGAGE : ESOPE 19 : C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/TTMF) 20 : C mél : gounand@semt2.smts.cea.fr 21 : C*********************************************************************** 22 : C APPELE : VEC2CH 23 : C*********************************************************************** 24 : C ENTREES : MATRIK, INCX, IMPR 25 : C ENTREES/SORTIES : - 26 : C SORTIES : MCHSOL, IRET 27 : C CODE RETOUR (IRET) : 0 si ok 28 : C <0 si problème 29 : C MATRIK : pointeur sur segment MATRIK de l'include SMMATRIK 30 : C on pioche dedans les informations nécessaires 31 : C (numérotations, inconnues, support géométriques) 32 : C pour effectuer la conversion MCHINI->INCX. 33 : C INCX : pointeur sur segment IZA de l'include SMMATRIK 34 : C vecteur généralement issu de la résolution 35 : C d'un système de type Ax=b. 36 : C IMPR : niveau d'impression 37 : C MCHTOT : pointeur sur segment MCHPOI de l'include SMCHPOI 38 : C chpoint résultat de la transformation 39 : C du vecteur INCX. 40 : C*********************************************************************** 41 : C VERSION : v1, 01/04/98, version initiale 42 : C HISTORIQUE : v1, 01/04/98, création 43 : C HISTORIQUE : 44 : C HISTORIQUE : 45 : C*********************************************************************** 46 : C Prière de PRENDRE LE TEMPS de compléter les commentaires 47 : C en cas de modification de ce sous-programme afin de faciliter 48 : C la maintenance ! 49 : C*********************************************************************** 50 : -INC CCOPTIO 51 : -INC SMMATRIK 52 : POINTEUR INCX.IZA 53 : -INC SMCHPOI 54 : POINTEUR MCHSOL.MCHPOI 55 : -INC SMELEME 56 : POINTEUR KISPG.MELEME 57 : C 58 : IRET=0 59 : SEGACT MATRIK 60 : C 61 : C On transforme le vecteur en chpoint 62 : C 63 : C Entrées : 64 : C INCX 65 : MINC =KMINCP 66 : KISPG =KISPGP 67 : IDMAT =KIDMAT(1) 68 : C 68 : CALL VEC2CH(INCX,MINC,KISPG,IDMAT,MCHPOI,IMPR,IRET) 69 : IF (IRET.NE.0) GOTO 9999 70 : C Sortie : 71 : MCHSOL= MCHPOI 72 : C 73 : IF (IMPR.GT.6) THEN 74 : WRITE(IOIMP,*) ' xdisp : on a créé le chpoint',MCHSOL 75 : ENDIF 76 : C 77 : SEGDES MATRIK 78 : * 79 : * Normal termination 80 : * 81 : RETURN 82 : * 83 : * Error handling 84 : * 85 : 9999 CONTINUE 86 : WRITE(IOIMP,*) 'An error was detected in xdisp.eso' 87 : RETURN 88 : * 89 : * End of XDISP 90 : * 91 : END 92 : 93 : 94 : 95 : 96 : 97 : 98 : 99 :
© Cast3M 2003 - All rights reserved.
Disclaimer