Sources Esope | Cast3M

maxro3
C MAXRO3    SOURCE    CHAT      05/01/13    01:37:53     5004
      SUBROUTINE MAXRO3(ICAS,IDEFO,S,A,R,RT,TRA,NDIM,IDIM,IFOUR)
*
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*
      DIMENSION R(NDIM,*),RT(NDIM,*),A(NDIM,*),TRA(NDIM,*)
      DIMENSION S(*)
*
 
*  ICAS = 1  :   RT A R
*  ICAS = 2  :   R A RT
 
 
 
 
            A(1,1) =  S(1)
            A(2,2) =  S(2)
            A(1,2) =  S(4)
            IF (IDEFO.EQ.1) A(1,2)=A(1,2)/2.D0
            A(2,1)=A(1,2)
*
            IF (IFOUR.LT.1) GOTO 666
*
            A(3,3) = S(3)
            A(3,1) = S(5)
            A(3,2) = S(6)
*
            IF (IDEFO.EQ.1) A(3,1)=A(3,1)/2.D0
            IF (IDEFO.EQ.1) A(3,2)=A(3,2)/2.D0
            A(1,3)=A(3,1)
            A(2,3)=A(3,2)
            A(3,3) = S(3)
*
 666        CONTINUE
*
 
            IF(ICAS.EQ.1) THEN
 
*                                           t
*              >>>  Rotation du tenseur : A = R  A  R  &lt;&lt;&lt;
*
               CALL MULMAT(TRA,A,R,NDIM,NDIM,NDIM)
               CALL MULMAT(A,RT,TRA,NDIM,NDIM,NDIM)
 
 
            ELSE IF(ICAS.EQ.2) THEN
 
*                                                t
*              >>>  Rotation du tenseur : A =  R A R   &lt;&lt;&lt;
*
               CALL MULMAT(TRA,A,RT,NDIM,NDIM,NDIM)
               CALL MULMAT(A,R,TRA,NDIM,NDIM,NDIM)
 
            ENDIF
 
            S(1) = A(1,1)
            S(2) = A(2,2)
*
            IF (IDEFO.EQ.1) A(1,2)=A(1,2)*2.D0
            S(4) = A(1,2)
*
            IF (IFOUR.GE.1) THEN
*
               S(3) = A(3,3)
*
               IF (IDEFO.EQ.1) A(3,1)=A(3,1)*2.D0
               IF (IDEFO.EQ.1) A(3,2)=A(3,2)*2.D0
               S(5) = A(3,1)
               S(6) = A(3,2)
*
            ENDIF
*
*
      RETURN
      END