Esope source files

   1 : C COUROT    SOURCE    CHAT      05/01/12    22:25:39     5004           
   2 :       SUBROUTINE COUROT(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
   7 : C  REL = TRANSPOSE(P).RP.P
   8 : C                                                                               
   9 : C  ENTREE                                                                       
  10 : C         LRE    NOMBRE DE LIGNES DU TABLEAU REL                                
  11 : C         P      MATRICE DE PASSAGE                                             
  12 : C         RP     MATRICE ELEMENTAIRE 10X10 EN AXES LOCAUX                       
  13 : C                                                                               
  14 : C  SORTIE                                                                       
  15 : C         REL    MATRICE ELEMENTAIRE EN AXES GLOBAUX                            
  16 : C-----------------------------------------------------------------------        
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  19 :       DIMENSION REL(LRE,LRE),P(LRE,LRE),RP(LRE,LRE),RPP(10,10)
  20 : C
  21 :       DO 1 I=1,100
  22 : 1     RPP(I,1)=0.
  23 : C
  24 :       DO 2 I=1,10
  25 :          DO 2 J=1,10
  26 :             DO 3 K=1,10
  27 : 3           RPP(I,J)=RPP(I,J)+RP(I,K)*P(K,J)                        
  28 : 2     CONTINUE
  29 : C
  30 :       DO 4 I=1,10
  31 :          DO 4 J=1,10
  32 :             DO 5 K=1,10
  33 : 5           REL(I,J)=REL(I,J)+P(K,I)*RPP(K,J)                        
  34 : 4     CONTINUE
  35 : C
  36 :       RETURN
  37 :       END
  38 : 
  39 :