Sources Esope | Cast3M

pourog
C POUROG    SOURCE    CB215821  26/08/24    21:17:52     12622          
      SUBROUTINE POUROG(REL,LRE,P,RP)
C-----------------------------------------------------------------------
C  ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES
C  GLOBAUX
C  POUR UNE MATRICE ANTISYMETRIQUE
C
C  ENTREE
C         LRE    NOMBRE DE LIGNES DU TABLEAU REL
C         P      MATRICE DE PASSAGE
C         RP     MATRICE ELEMENTAIRE 12X12 EN AXES LOCAUX
C
C  SORTIE
C         REL    MATRICE ELEMENTAIRE EN AXES GLOBAUX
C
C
C     Didier COMBESCURE Mars 2003
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION REL(12,12),P(3,3),RP(12,12)
      DIMENSION NI(12),NJ(12)
C
      DATA NI/1,1,1,4,4,4,7,7,7,10,10,10/
      DATA NJ/1,2,3,1,2,3,1,2,3,1,2,3/
C
       DO 13 I=1,12
        DO 12 J=1,12
         REL(I,J)=  0.D0
C
         DO 11 K=1,3
          DO 10 L=1,3
C
           I1 = NI(I) + K - 1
           J1 = NI(J) + L - 1
           I2 = NJ(I)
           J2 = NJ(J)
C
           REL(I,J)=REL(I,J) + (P(K,I2)*RP(I1,J1)*P(L,J2))
C
   10     CONTINUE
   11    CONTINUE
   12   CONTINUE
   13  CONTINUE
C
C
      RETURN
      END