1 : C POUVEC SOURCE CHAT 05/01/13 02:20:50 5004 2 : SUBROUTINE POUVEC(U,V,P,ICOD) 3 : C----------------------------------------------------------------------- 4 : C ROUTINE DE PASSAGE DU VECTEUR DEPLACEMENT GLOBAL AU LOCAL 5 : C OU DU VECTEUR FORCE LOCAL AU GLOBAL 6 : C 7 : C ENTREE 8 : C ------ U DEPLACEMENTS GLOBAUX AUX 2 NOEUDS 9 : C OU FORCES LOCALES AUX 2 NOEUDS 10 : C P MATRICE DE PASSAGE 11 : C ICOD 1 POUR Q-LOC =P*Q-GLOB 12 : C 2 POUR F-GLOB=PT*F-LOC 13 : C 14 : C SORTIE 15 : C ------ V DEPLACEMENTS LOCAUX AUX 2 NOEUDS 16 : C OU FORCES GLOBALES AUX 2 NOEUDS 17 : C----------------------------------------------------------------------- 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : DIMENSION U(1),V(1),P(3,3) 21 : C 22 : GO TO (100,200),ICOD 23 : 100 CONTINUE 24 : DO 6 K=1,4 25 : IDEC=(K-1)*3 26 : DO 6 J=1,3 27 : JJ=J+IDEC 28 : V( JJ)=0.D00 29 : DO 6 I=1,3 30 : II=I+IDEC 31 : V(JJ)=V(JJ)+P(J,I)*U(II) 32 : 6 CONTINUE 33 : RETURN 34 : 200 CONTINUE 35 : DO 7 K=1,4 36 : IDEC=(K-1)*3 37 : DO 7 J=1,3 38 : JJ=J+IDEC 39 : V(JJ)=0.D00 40 : DO 7 I=1,3 41 : II=I+IDEC 42 : V(JJ)=V(JJ)+P(I,J)*U(II) 43 : 7 CONTINUE 44 : RETURN 45 : END 46 :
© Cast3M 2003 - All rights reserved.
Disclaimer