1 : C VCH1 SOURCE PV 10/03/31 21:15:00 6654 2 : SUBROUTINE VCH1 (MMATRX,MVECTX,ISOLU,KRIGI) 3 : C 4 : C **** A PARTIR D UN OBJET DE TYPE MATRICE TRANSFORME UN VECTEUR EN 5 : C **** CHPOIN DE TYPE PREMIER MEMBRE 6 : C **** KRIGI EST LE POINTEUR DE LA RIGIDITE.SERT A SAVOIRE SI LA GEOM 7 : C **** EXISTE DEJA. 8 : C 9 : IMPLICIT INTEGER(I-N) 10 : -INC SMMATRI 11 : -INC SMCHPOI 12 : -INC SMELEME 13 : -INC CCOPTIO 14 : -INC SMVECTD 15 : -INC SMRIGID 16 : -INC TMTRAV 17 : C 18 : IF(IIMPI.EQ.3) WRITE(IOIMP,3000) MMATRX,MVECTX 19 : 3000 FORMAT(' OPERATEUR VCH1 : POINTEUR DE LA MATRICE=',I5, 20 : 1 ' POINTEUR DU VECTEUR=',I5) 21 : C 22 : C **** ACTIVATION DES SEGMENTS 23 : C 24 : MRIGID=KRIGI 25 : SEGACT,MRIGID 26 : IMGEOQ=IMGEO1 27 : 28 : MMATRI=MMATRX 29 : SEGACT,MMATRI 30 : MILIGN=IILIGN 31 : SEGACT,MILIGN 32 : INC=IPNO(/1) 33 : SEGDES,MILIGN 34 : MINCPO=IINCPO 35 : SEGACT,MINCPO 36 : MIMIK=IIMIK 37 : SEGACT,MIMIK 38 : MHARK=IHARK 39 : SEGACT,MHARK 40 : NNIN=IMIK(/2) 41 : midua=iidua 42 : segact midua 43 : MVECTD=MVECTX 44 : SEGACT MVECTD 45 : C 46 : C **** CREATION D'UN SEGMENT DE TYPE MCHPOI(VOIR SMCHPOI) 47 : C 48 : 49 : if( ivecri.eq.0) then 50 : NND=VECTBB(/1) 51 : MELEME=IGEOMA 52 : SEGACT,MELEME 53 : NNNOE=NUM(/2) 54 : SEGINI,MTRAV 55 : DO 40 ITYU=1,NNNOE 56 : IGEO(ITYU)=NUM(1,ITYU) 57 : 40 CONTINUE 58 : SEGDES,MELEME 59 : DO 41 ITYU=1,NNIN 60 : NHAR(ITYU)=IHAR(ITYU) 61 : 41 INCO(ITYU)=IMIK(ITYU) 62 : DO 23 I=1,NNNOE 63 : DO 23 J=1,NNIN 64 : IK=INCPO(J,I) 65 : IF(IK.EQ.0) GO TO 23 66 : IBIN(J,I)=ik 67 : BB(J,I)=VECTBB(IK) 68 : 23 CONTINUE 69 : segact mrigid*mod 70 : ichoa=ichole 71 : if(ichole.eq.0) then 72 : ichole = mmatri 73 : endif 73 : CALL CRECH3(MTRAV,ISOLU,nnd,mrigid) 74 : ichole=ichoa 75 : SEGSUP,MTRAV 76 : else 77 : lvecri=ivecri 77 : CALL CRECH2(ISOLU,mvectx,lvecri,1) 78 : endif 79 : 80 : SEGDES,MINCPO 81 : SEGDES,MIMIK 82 : segdes midua 83 : SEGDES,MMATRI 84 : SEGDES,MHARK 85 : * SEGDES MVECTD 86 : SEGDES,MRIGID 87 : IF(IIMPI.EQ.3) WRITE(IOIMP,3001) ISOLU 88 : 3001 FORMAT(' SUBROUTINE VCH1 : POINTEUR DE L'' OBJET CHPOINT=',I5) 89 : RETURN 90 : END 91 : 92 : 93 : 94 : 95 : 96 : 97 : 98 : 99 : 100 : 101 : 102 :
© Cast3M 2003 - All rights reserved.
Disclaimer