1 : C DYNERE SOURCE CHAT 05/01/12 23:18:54 5004 2 : SUBROUTINE DYNERE(ICHT12,IPZ,ICHT13,IDDL) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : ************************************************************************ 6 : * 7 : * D Y N E R E 8 : * ----------- 9 : * 10 : * FONCTION: 11 : * --------- 12 : * UTILITAIRE DE L'OPERATEUR "DYNE" 13 : * RESOUD LE SYSTEME : Z L1 = C1 14 : * 15 : * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN) 16 : * ----------- 17 : * 18 : * VEC (E) VECTEUR CONTENANT LES MULTIPLICATEURS 19 : * IPZ (E) POINTEUR SUR LE SEGMENT MZZZ ( ACTIF ) 20 : * VEC1 (S) VECTEUR RESULTAT 21 : * IDDL (E) NOMBRE DE DDL LIES 22 : * 23 : * REMARQUE: 24 : * --------- 25 : * L'OBJET RIGIDITE Z QUI EST FABRIQUE DANS DYNE02 POURRA ETRE 26 : * UTILISE A LA PLACE DE MZZZ LORSQU'ON INTRODUIRA LES RELATIONS 27 : * 28 : * AUTEUR, DATE DE CREATION: 29 : * ------------------------- 30 : * DENIS ROBERT, LE 1ER DECEMBRE 1988. 31 : * REVU ET CORRIGE DECEMBRE 89 MP 32 : * 33 : * LANGAGE: 34 : * -------- 35 : * ESOPE + FORTRAN77 36 : * 37 : ************************************************************************ 38 : * 39 : PARAMETER( UN=1.D0 ) 40 : * 41 : SEGMENT,MZZZ 42 : REAL*8 ZZZZ(LONZ) 43 : ENDSEGMENT 44 : * 45 : SEGMENT,MCHT1 46 : REAL*8 VCC1(IDDL) 47 : ENDSEGMENT 48 : POINTEUR MCHT12.MCHT1,MCHT13.MCHT1 49 : * 50 : * REAL*8 VEC(*),VEC1(*) 51 : * 52 : MCHT12=ICHT12 53 : MCHT13=ICHT13 54 : MZZZ=IPZ 55 : * 56 : DO 10 I=1,IDDL 57 : IKO=I*(I+1)/2 58 : MCHT13.VCC1(I)=(UN/ZZZZ(IKO))*MCHT12.VCC1 59 : 10 CONTINUE 60 : * END DO 61 : * 62 : END 63 :
© Cast3M 2003 - All rights reserved.
Disclaimer