matvec
C MATVEC SOURCE CHAT 05/01/13 01:35:46 5004 C======================================================================= C ROUTINE DE PASSAGE DU VECTEUR DEPLACEMENT GLOBAL AU LOCAL C C C ENTREE C U(3*NB) = DEPLACEMENTS GLOBAUX AUX 2 NOEUDS C P(3,3) = MATRICE DE PASSAGE C NB = NOMBRE DE VECTEURS A 3 LIGNES CONTENUS DANS U ET V C SORTIE C V(3*NB) = DEPLACEMENTS LOCAUX AUX 2 NOEUDS C======================================================================= IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION U(*),V(*),P(3,*) DATA XZER/0.D0/ C DO 600 J=1,3 KINC=0 DO 700 K=1,NB V(J+KINC)=XZER C DO 500 I=1,3 V(J+KINC)=V(J+KINC)+P(J,I)*U(I+KINC) 500 CONTINUE KINC=KINC+3 700 CONTINUE 600 CONTINUE RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales