C MATVEC    SOURCE    CHAT      05/01/13    01:35:46     5004
      SUBROUTINE MATVEC(U,V,P,NB)
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