Esope source files

   1 : C DEINIJ    SOURCE    CHAT      05/01/12    22:41:11     5004           
   2 : 
   3 :       SUBROUTINE DERINVARIIJJ (SIG,I_S,J2_S,J3_S,N)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer n,i
   6 :       real*8    sig(n),I_s(n),J2_s(n),J3_s(n),aux1,J2
   6 :       call Desviador(sig,J2_s,n)
   6 :       call Invari_J2(sig,n,J2)
   7 :       J3_s(1)=J2_s(1)*J2_s(1)
   8 :       J3_s(2)=J2_s(2)*J2_s(2)
   9 :       J3_s(3)=J2_s(3)*J2_s(3)
  10 :       if (n.ge.4) then
  11 :        J3_s(1)=J3_s(1)+J2_s(4)*J2_s(4)
  12 :        J3_s(2)=J3_s(2)+J2_s(4)*J2_s(4)
  13 :        J3_s(4)=2.D0*(J2_s(1)+J2_s(2))*J2_s(4)
  14 :        if (n.eq.6) then
  15 :          J3_s(1)=J3_s(1)+J2_s(5)*J2_s(5)
  16 :          J3_s(2)=J3_s(2)+J2_s(6)*J2_s(6)
  17 :          J3_s(3)=J3_s(3)+J2_s(5)*J2_s(5)+J2_s(6)*J2_s(6)
  18 :          J3_s(4)=J3_s(4)+J2_s(5)*J2_s(6)
  19 :          J3_s(5)=(J2_s(1)+J2_s(3))*J2_s(5)+J2_s(4)*J2_s(6)
  20 :          J3_s(6)=(J2_s(2)+J2_s(3))*J2_s(6)+J2_s(4)*J2_s(5)
  21 :        endif
  22 :       endif
  23 :       aux1=-J2*2.D0/3.D0
  24 :       do i=1,3
  25 :          J3_s(i)=J3_s(i)+aux1
  26 :       enddo
  27 :       do i=1,3
  28 :          I_s(i)=1.D0
  29 :       enddo
  30 :       if (n.lt.4) return
  31 :       do i=4,n
  32 :          I_s(i)=0.D0
  33 :          J2_s(i)=2.D0*J2_s(i)
  34 :       enddo
  35 :       return
  36 :       end
  37 :  
  38 :  
  39 :