1 : C DIMEN3 SOURCE CHAT 09/10/09 21:17:14 6519 2 : SUBROUTINE DIMEN3 (IPRIGI,NBRINC) 3 : ************************************************************************ 4 : * 5 : * D I M E N 3 6 : * ----------- 7 : * 8 : * FONCTION: 9 : * --------- 10 : * 11 : * DONNER LE NOMBRE D'INCONNUES DU PROBLEME PHYSIQUE ASSOCIE A UN 12 : * OBJET 'RIGIDITE'. 13 : * 14 : * MODE D'APPEL: 15 : * ------------- 16 : * 16 : * CALL DIMEN3 (IPRIGI,NBRINC) 17 : * 18 : * PARAMETRES: (E)=ENTREE (S)=SORTIE 19 : * ----------- 20 : * 21 : * IPRIGI ENTIER (E) POINTEUR D'UNE 'RIGIDITE'. 22 : * NBRINC ENTIER (S) NOMBRE D'INCONNUES DE LA 'RIGIDITE'. 23 : * 24 : * REMARQUES: 25 : * ---------- 26 : * 27 : * CE SOUS-PROGRAMME FACTORISE L'OBJET 'RIGIDITE', SI CE N'EST DEJA 28 : * FAIT. 29 : * 30 : * AUTEUR, DATE DE CREATION: 31 : * ------------------------- 32 : * 33 : * PASCAL MANIGOT 26 AOUT 1985 34 : * 35 : * LANGAGE: 36 : * -------- 37 : * 38 : * ESOPE + FORTRAN77 39 : * 40 : ************************************************************************ 41 : * 42 : IMPLICIT INTEGER(I-N) 43 : -INC CCOPTIO 44 : -INC SMMATRI 45 : -INC SMRIGID 46 : * 47 : integer insym 48 : insym = 0 49 : MRIGID = IPRIGI 50 : SEGACT,MRIGID 51 : NRG = IRIGEL(/1) 52 : NBR = IRIGEL(/2) 53 : IPMATR = ICHOLE 54 : c 55 : * 56 : IF(NORINC.GT.0 .AND. NORIND.GT.0) THEN 57 : INSYM = 1 58 : ENDIF 59 : IF (NRG.GE.7) THEN 60 : DO 9 IN = 1,NBR 61 : IANTI=IRIGEL(7,IN) 62 : IF(IANTI.GT.0) THEN 63 : INSYM = 1 64 : ENDIF 65 : 9 CONTINUE 66 : ENDIF 67 : C 68 : SEGDES,MRIGID 69 : IF (IPMATR .EQ. 0) THEN 70 : IF (INSYM .EQ. 0) THEN 70 : CALL TRIANG (IPRIGI,1D-18) 71 : ELSE 71 : CALL ldmt1(IPRIGI,1d-18) 72 : ENDIF 73 : IF (IERR .NE. 0) RETURN 74 : MRIGID = IPRIGI 75 : SEGACT,MRIGID 76 : IPMATR = ICHOLE 77 : SEGDES,MRIGID 78 : END IF 79 : * 80 : MMATRI = IPMATR 81 : SEGACT,MMATRI 82 : MILIGN = IILIGN 83 : SEGACT,MILIGN 84 : NBRINC = IPNO(/1) 85 : SEGDES,MILIGN 86 : SEGDES,MMATRI 87 : * 88 : END 89 : 90 : 91 : 92 : 93 : 94 : 95 : 96 : 97 :
© Cast3M 2003 - All rights reserved.
Disclaimer