Sources Esope | Cast3M

chrep
C CHREP     SOURCE    CHAT      05/01/12    22:02:01     5004
      SUBROUTINE CHREP(ENT,SOR,ANG)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION ENT(3),SOR(3)
C
C----------------------------------------------
C     CETTE ROUTINE PERMET DE CHANGER DE      |
C     REPERE POUR CONTRAINTE ET DEFORMATION   |
C----------------------------------------------
C
      S=SIN(ANG)
      C=COS(ANG)
      CC=C*C
      SS=S*S
      CS=S*C
      SOR(1)=CC*ENT(1)+SS*ENT(2)+2.D0*CS*ENT(3)
      SOR(2)=CC*ENT(2)+SS*ENT(1)-2.D0*CS*ENT(3)
      SOR(3)=CS*(ENT(2)-ENT(1))+(CC-SS)*ENT(3)
      RETURN
      END