Esope source files

   1 : C INVAJ3    SOURCE    CHAT      05/01/13    00:42:12     5004           
   2 : 
   3 :       SUBROUTINE INVARI_J3 (SIG,N,J3)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer n,i
   6 :       real*8 sig(n),J3,desv(6),desvdesv(6)
   6 :       call Desviador(sig,desv,n)
   7 :       desvdesv(1)=desv(1)*desv(1)
   8 :       desvdesv(2)=desv(2)*desv(2)
   9 :       desvdesv(3)=desv(3)*desv(3)
  10 :       if (n.ge.4) then
  11 :        desvdesv(1)=desvdesv(1)+desv(4)*desv(4)
  12 :        desvdesv(2)=desvdesv(2)+desv(4)*desv(4)
  13 :        desvdesv(4)=(desv(1)+desv(2))*desv(4)
  14 :        if (n.eq.6) then
  15 :          desvdesv(1)=desvdesv(1)+desv(5)*desv(5)
  16 :          desvdesv(2)=desvdesv(2)+desv(6)*desv(6)
  17 :          desvdesv(3)=desvdesv(3)+desv(5)*desv(5)+desv(6)*desv(6)
  18 :          desvdesv(4)=desvdesv(4)+desv(5)*desv(6)
  19 :          desvdesv(5)=(desv(1)+desv(3))*desv(5)+desv(4)*desv(6)
  20 :          desvdesv(6)=(desv(2)+desv(3))*desv(6)+desv(4)*desv(5)
  21 :        endif
  22 :       endif
  23 :       J3=0.D0
  24 :       do i=1,3
  25 :          J3=J3+desvdesv(i)*desv(i)
  26 :       enddo
  27 :       J3=J3/3.D0
  28 :       if (n.lt.4) return
  29 :       do i=4,n
  30 :          J3=J3+2.D0*desvdesv(i)*desv(i)/3.D0
  31 :       enddo
  32 :       return
  33 :       end
  34 :  
  35 :  
  36 :