Esope source files

   1 : C POUROG    SOURCE    CHAT      05/01/13    02:19:50     5004           
   2 :       SUBROUTINE POUROG(REL,LRE,P,RP)                                           
   3 : C-----------------------------------------------------------------------        
   4 : C  ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES        
   5 : C  GLOBAUX                                                                      
   6 : C  POUR UNE MATRICE ANTISYMETRIQUE                                                              
   7 : C                                                                               
   8 : C  ENTREE                                                                       
   9 : C         LRE    NOMBRE DE LIGNES DU TABLEAU REL                                
  10 : C         P      MATRICE DE PASSAGE                                             
  11 : C         RP     MATRICE ELEMENTAIRE 12X12 EN AXES LOCAUX                       
  12 : C                                                                               
  13 : C  SORTIE                                                                       
  14 : C         REL    MATRICE ELEMENTAIRE EN AXES GLOBAUX                            
  15 : C                                                                               
  16 : C 
  17 : C     Didier COMBESCURE Mars 2003                                                                
  18 : C-----------------------------------------------------------------------        
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  21 :       DIMENSION REL(12,12),P(3,3),RP(12,12)                                     
  22 :       DIMENSION NI(12),NJ(12)                                                   
  23 : C                                                                               
  24 :       DATA NI/1,1,1,4,4,4,7,7,7,10,10,10/                                       
  25 :       DATA NJ/1,2,3,1,2,3,1,2,3,1,2,3/ 
  26 : C
  27 :        DO 10 I=1,12  
  28 :         DO 10 J=1,12
  29 :          REL(I,J)=  0.D0                
  30 : C                                        
  31 :          DO 10 K=1,3  
  32 :           DO 10 L=1,3
  33 : C
  34 :            I1 = NI(I) + K - 1
  35 :            J1 = NI(J) + L - 1
  36 :            I2 = NJ(I)
  37 :            J2 = NJ(J)
  38 : C
  39 :            REL(I,J)=REL(I,J) + (P(K,I2)*RP(I1,J1)*P(L,J2))
  40 : C   
  41 :    10     CONTINUE    
  42 : C
  43 : C                                                               
  44 :       RETURN                                                                    
  45 :       END             
  46 :                                                    
  47 :  
  48 :  
  49 :  
  50 :