Sources Esope | Cast3M

courot
C COUROT    SOURCE    PV        18/06/20    21:15:05     9860           
      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 1 J=1,10
      DO 1 I=1,10
1     RPP(I,J)=0.
C
      DO 2 I=1,10
         DO 2 J=1,10
            DO 3 K=1,10
3           RPP(I,J)=RPP(I,J)+RP(I,K)*P(K,J)
2     CONTINUE
C
      DO 4 I=1,10
         DO 4 J=1,10
            DO 5 K=1,10
5           REL(I,J)=REL(I,J)+P(K,I)*RPP(K,J)
4     CONTINUE
C
      RETURN
      END