Sources Esope | Cast3M

deinij
C DEINIJ    SOURCE    CHAT      05/01/12    22:41:11     5004
 
      SUBROUTINE DERINVARIIJJ (SIG,I_S,J2_S,J3_S,N)
      IMPLICIT INTEGER(I-N)
      integer n,i
      real*8    sig(n),I_s(n),J2_s(n),J3_s(n),aux1,J2
      call Desviador(sig,J2_s,n)
      call Invari_J2(sig,n,J2)
      J3_s(1)=J2_s(1)*J2_s(1)
      J3_s(2)=J2_s(2)*J2_s(2)
      J3_s(3)=J2_s(3)*J2_s(3)
      if (n.ge.4) then
       J3_s(1)=J3_s(1)+J2_s(4)*J2_s(4)
       J3_s(2)=J3_s(2)+J2_s(4)*J2_s(4)
       J3_s(4)=2.D0*(J2_s(1)+J2_s(2))*J2_s(4)
       if (n.eq.6) then
         J3_s(1)=J3_s(1)+J2_s(5)*J2_s(5)
         J3_s(2)=J3_s(2)+J2_s(6)*J2_s(6)
         J3_s(3)=J3_s(3)+J2_s(5)*J2_s(5)+J2_s(6)*J2_s(6)
         J3_s(4)=J3_s(4)+J2_s(5)*J2_s(6)
         J3_s(5)=(J2_s(1)+J2_s(3))*J2_s(5)+J2_s(4)*J2_s(6)
         J3_s(6)=(J2_s(2)+J2_s(3))*J2_s(6)+J2_s(4)*J2_s(5)
       endif
      endif
      aux1=-J2*2.D0/3.D0
      do i=1,3
         J3_s(i)=J3_s(i)+aux1
      enddo
      do i=1,3
         I_s(i)=1.D0
      enddo
      if (n.lt.4) return
      do i=4,n
         I_s(i)=0.D0
         J2_s(i)=2.D0*J2_s(i)
      enddo
      return
      end