Sources Esope | Cast3M

courot
C COUROT    SOURCE    CB215821  26/08/24    21:15:49     12622          
      SUBROUTINE COUROT(REL,LRE,P,RP)
C-----------------------------------------------------------------------
C  ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES
C  GLOBAUX
C
C  REL = TRANSPOSE(P).RP.P
C
C  ENTREE
C         LRE    NOMBRE DE LIGNES DU TABLEAU REL
C         P      MATRICE DE PASSAGE
C         RP     MATRICE ELEMENTAIRE 10X10 EN AXES LOCAUX
C
C  SORTIE
C         REL    MATRICE ELEMENTAIRE EN AXES GLOBAUX
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION REL(LRE,LRE),P(LRE,LRE),RP(LRE,LRE),RPP(10,10)
C
      DO 6 J=1,10
      DO 1 I=1,10
      RPP(I,J)=0.D0
    1 CONTINUE
    6 CONTINUE
C
      DO 7 I=1,10
         DO 2 J=1,10
            DO 3 K=1,10
            RPP(I,J)=RPP(I,J)+RP(I,K)*P(K,J)
    3       CONTINUE
    2    CONTINUE
    7 CONTINUE
C
      DO 8 I=1,10
         DO 4 J=1,10
            DO 5 K=1,10
            REL(I,J)=REL(I,J)+P(K,I)*RPP(K,J)
    5       CONTINUE
    4    CONTINUE
    8 CONTINUE
C
      RETURN
      END