Esope source files

   1 : C DIACL     SOURCE    CHAT      05/01/12    22:49:34     5004           
   2 :         SUBROUTINE DIACL(EPS,SIN,COS)
   3 :       IMPLICIT INTEGER(I-N)
   4 :         IMPLICIT REAL*8(A-H,O-Z)
   5 : C
   6 : C       EPS = TABLEAU DES DEFORMATIONS (EPS11,EPS22,SQRT(2)*EPS12)
   7 : C       EPSL1 = DEFORMATIONS PRINCIPALES
   8 : C       SIN ET COS = SINUS ET COSINUS DE L"ANGLE DE ROTATION
   9 :         DIMENSION  EPS(3)
  10 : C
  11 :         EPS12=EPS(3)/2.D0
  12 :         X1=2.D0*EPS12
  13 :         X2=EPS(1)-EPS(2)
  14 :         X3=SQRT(X1**2+X2**2)
  15 :         X8=X2*1.D+10
  16 :         IF (DABS(X1).GT.DABS(X8).OR.X2.EQ.0.D0) THEN
  17 :           COS=0.5D0*SQRT(2.D0)
  18 :           SIN=SIGN(COS,X1)
  19 :         ELSE
  20 :           X4=X1/X2
  21 :           X5=1.D0+X4**2
  22 :           X5=SIGN(1.D0,X2)/SQRT(X5)
  23 :           X6=(1.D0+X5)*0.5D0
  24 :           X6=SQRT(X6)
  25 :           X7=(1.D0-X5)*0.5D0
  26 :           X7=SIGN(1.D0,X1)*SQRT(X7)
  27 :           COS=X6
  28 :           SIN=X7
  29 :         END IF
  30 : C
  31 :         RETURN
  32 :         END
  33 :