Esope source files

   1 : C MATVEC    SOURCE    CHAT      05/01/13    01:35:46     5004           
   2 :       SUBROUTINE MATVEC(U,V,P,NB)                                               
   3 : C=======================================================================        
   4 : C  ROUTINE DE PASSAGE DU VECTEUR DEPLACEMENT GLOBAL AU LOCAL                    
   5 : C                                                                               
   6 : C                                                                               
   7 : C  ENTREE                                                                       
   8 : C         U(3*NB) =  DEPLACEMENTS GLOBAUX AUX 2 NOEUDS                          
   9 : C         P(3,3)  =  MATRICE DE PASSAGE                                         
  10 : C         NB      =  NOMBRE DE VECTEURS A 3 LIGNES CONTENUS DANS U ET V         
  11 : C  SORTIE                                                                       
  12 : C         V(3*NB) =  DEPLACEMENTS LOCAUX AUX 2 NOEUDS                           
  13 : C=======================================================================        
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  16 :       DIMENSION U(*),V(*),P(3,*)                                                
  17 :       DATA XZER/0.D0/                                                           
  18 : C                                                                               
  19 :       DO 600 J=1,3                                                              
  20 :       KINC=0                                                                    
  21 :       DO 700 K=1,NB                                                             
  22 :       V(J+KINC)=XZER                                                            
  23 : C                                                                               
  24 :       DO 500 I=1,3                                                              
  25 :       V(J+KINC)=V(J+KINC)+P(J,I)*U(I+KINC)                                      
  26 :   500 CONTINUE                                                                  
  27 :       KINC=KINC+3                                                               
  28 :   700 CONTINUE                                                                  
  29 :   600 CONTINUE                                                                  
  30 :       RETURN                                                                    
  31 :       END                                                                       
  32 :