1 : C MTREDU SOURCE CHAT 09/10/09 21:21:16 6519 2 : C SUBROUTINE MTREDU ( IPX, IPRIGI, IPMASS, IPRES) 3 : ************************************************************************ 4 : * 5 : * MTREDU 6 : * ----------- 7 : * 8 : * FONCTION: 9 : * --------- 10 : * 11 : * EFFECTUE L'OPERATION IPRES = (K**-1)M IPX 12 : * 13 : * MODE D'APPEL: 14 : * ------------- 15 : * 15 : * CALL MTREDU ( IPX, IPRIGI, IPMASS, IPRES) 16 : * 17 : * PARAMETRES: (E)=ENTREE (S)=SORTIE 18 : * ----------- 19 : * 20 : * IPX ENTIER (E) POINTEUR DE L'OBJET CHPOIN CONTENANT 21 : * LE VECTEUR A MULTIPLIER 22 : * IPRIGI ENTIER (E) POINTEUR SUR L'OBJET 'RIGIDITE' K 23 : * IPMASS ENTIER (E) POINTEUR SUR L'OBJET 'RIGIDITE' M 24 : * IPRES ENTIER (S) POINTEUR SUR L'OBJET CHPOIN CONTENANT 25 : * LE VECTEUR SOLUTION 26 : * 27 : * 28 : * AUTEUR, DATE DE CREATION: 29 : * ------------------------- 30 : * 31 : * C. LE BIDEAU 09 / 2001 ( FORTRAN + ESOPE ) 32 : * 33 : *********************************************************************** 34 : 35 : SUBROUTINE MTREDU ( IPX, IPRIGI, IPMASS, IPRES) 36 : IMPLICIT INTEGER(I-N) 37 : IMPLICIT REAL*8 (A-H, O-Z) 38 : 39 : -INC CCOPTIO 40 : -INC SMLMOTS 41 : -INC SMCHPOI 42 : -INC SMRIGID 43 : -INC CCHAMP 44 : 45 : SEGMENT IDEMEN (0) 46 : 47 : 48 : 49 : IPRES = 0 49 : CALL MUCPRI (IPX, IPMASS, IPBV) 50 : SEGINI IDEMEN 51 : IDEMEN(**)=IPBV 52 : NOID=0 53 : NOEN=1 53 : CALL LDMT (IPRIGI,IDEMEN,NOID,NOEN,1.D-15) 54 : SEGACT IDEMEN 55 : IPRES = IDEMEN(1) 56 : SEGSUP IDEMEN 57 : IF (IERR .NE. 0) RETURN 58 : 59 : END 60 : 61 : 62 : 63 : 64 : 65 : 66 :
© Cast3M 2003 - All rights reserved.
Disclaimer