Sources Esope | Cast3M

ottopr
C OTTOPR    SOURCE    CHAT      05/01/13    02:07:47     5004
      SUBROUTINE OTTOPR(D,DDEF,DSIGMA)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*
      PARAMETER (DEUX=2.D0)
*
      DIMENSION DDEF(6),DSIGMA(6),D(18)
*
      DSIGMA(1)=D(1)*DDEF(1)+D(6)*DDEF(2)+D(5)*DDEF(3)
      DSIGMA(2)=D(6)*DDEF(1)+D(2)*DDEF(2)+D(4)*DDEF(3)
      DSIGMA(3)=D(5)*DDEF(1)+D(4)*DDEF(2)+D(3)*DDEF(3)
      DSIGMA(4)=D(10)*DDEF(1)+D(11)*DDEF(2)+D(12)*DDEF(3)
     &          +DEUX*D(7)*DDEF(4)
      DSIGMA(5)=D(13)*DDEF(1)+D(14)*DDEF(2)+D(15)*DDEF(3)
     &          +DEUX*D(8)*DDEF(5)
      DSIGMA(6)=D(16)*DDEF(1)+D(17)*DDEF(2)+D(18)*DDEF(3)
     &          +DEUX*D(9)*DDEF(6)
*
      RETURN
      END