Sources Esope | Cast3M

chrepe
C CHREPE    SOURCE    CHAT      05/01/12    22:01:58     5004
      SUBROUTINE CHREPE(CO,SI,ADON,FDON)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION ADON(*),FDON(*)
C -------------------------------------------------------------
C      CHANGEMENT DE REPERE (CONTRAINTES OU DEFORMATIONS)
C -------------------------------------------------------------
C
      FDON(1)=CO*CO*ADON(1)+SI*SI*ADON(2)+2.D0*CO*SI*ADON(4)
      FDON(2)=SI*SI*ADON(1)+CO*CO*ADON(2)-2.D0*CO*SI*ADON(4)
      FDON(3)=ADON(3)
      FDON(4)=SI*CO*(-ADON(1)+ADON(2))+(CO*CO-SI*SI)*ADON(4)
      END