Esope source files

   1 : C POUROT    SOURCE    CHAT      05/01/13    02:20:19     5004           
   2 :       SUBROUTINE POUROT(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  D APRES TEDEL                                                                
   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     VIENT DE TEDEL AOUT 84                                                    
  17 : C-----------------------------------------------------------------------        
  18 :       IMPLICIT INTEGER(I-N)
  19 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  20 :       DIMENSION REL(LRE,1),P(3,3),RP(12,12)                                     
  21 :       DIMENSION NI(12),NJ(12)                                                   
  22 : C                                                                               
  23 :       DATA NI/1,1,1,4,4,4,7,7,7,10,10,10/                                       
  24 :       DATA NJ/1,2,3,1,2,3,1,2,3,1,2,3/                                          
  25 :        DO 10 K=1,12                                                             
  26 :       I3=NI(K)                                                                  
  27 :        I4=I3+2                                                                  
  28 :       LK=NJ(K)                                                                  
  29 :       DO 11 L=K,12                                                              
  30 :       REL(K,L)=  0.                                                             
  31 :        I1=NI(L)                                                                 
  32 :       I2=I1+2                                                                   
  33 :       LL=NJ(L)                                                                  
  34 :       LM=0                                                                      
  35 :        DO 12 M=I3,I4                                                            
  36 :       LM=LM+1                                                                   
  37 :       S=0.                                                                      
  38 :       II=0                                                                      
  39 :       DO 13 N=I1,I2                                                             
  40 :       II=II+1                                                                   
  41 :    13 S=S+  RP(M,N)*P(II,LL)                                                    
  42 :    12 REL(K,L)= REL(K,L)+ S*P( LM,LK)                                           
  43 :    11 REL(L,K)= REL(K,L)                                                        
  44 :    10 CONTINUE                                                                  
  45 :       RETURN                                                                    
  46 :       END                                                                       
  47 :