Esope source files

   1 : C WXWYSR    SOURCE    CHAT      05/01/13    04:13:08     5004           
   2 :       SUBROUTINE WXWYSR(WX,WY,DJAC,LREX,LREY,NST,REL)                                  
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C      CALCUL DE LA MUTUELLE ENTRE 2 TRIANGLES
   6 : C      VERSION Yann STEPHAN FEVRIER 97
   7 : C      ROUTINE FORTRAN PUR                                                      
   8 : C                                                                               
   9 : C=======================================================================        
  10 : C     WX(NST,LREX)=FONCTIONS DE FORME AU POINT X
  11 : C     WY(NST,LREY)=FONCTIONS DE FORME AU POINT Y
  12 : C     DJAC =JACOBIEN                                                         
  13 : C     REL=MATRICE MUTUELLE                                                         
  14 : C           LES MIJ SONT AJOUTES SUR REL                                        
  15 : C=======================================================================        
  16 :       IMPLICIT INTEGER(I-N)
  17 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  18 :       DIMENSION REL(LREX,1),WX(NST,1),WY(NST,1)
  19 :       DIMENSION WORK1(10)
  20 :       DO   300 IX=1,LREX                                                         
  21 : C                                                                               
  22 : C     BOUCLE SUR LES DDL                                                        
  23 : C                                                                               
  24 :       DO   400 INA=1,NST                                                        
  25 : C                                                                               
  26 : C     CALCUL DE JACOBIAN*POIDS*B   POUR LE DDL IX                               
  27 : C                                                                               
  28 :       WORK1(INA)=DJAC*WX(INA,IX)                                             
  29 :   400 CONTINUE                                                                  
  30 :       DO   700 IY=1,LREY                                                          
  31 : C                                                                               
  32 :       CC=0.D0                                                                   
  33 :       DO   800 IND=1,NST                                                        
  34 :       CC=CC+WY(IND,IY)*WORK1(IND)                                            
  35 :   800 CONTINUE                                                                  
  36 :       REL(IY,IX)=REL(IY,IX)+CC                                         
  37 :   700 CONTINUE                                                                  
  38 :   300 CONTINUE                                                                  
  39 :       RETURN                                                                    
  40 :       END                                                                       
  41 :  
  42 :  
  43 :