Esope source files

   1 : C CINVAT    SOURCE    CHAT      05/01/12    22:02:47     5004           
   2 : 
   3 :       SUBROUTINE C_INVARI_T (SIG,N,T)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer n
   6 :       complex*16 sig(n),t,J2,J3
   6 :       call C_Invari_J2(sig,n,J2)
   7 :       t=0.D0
   8 :       if (abs(J2).gt.1.D-17) then
   8 :          call C_Invari_J3(sig,n,J3)
   9 :          t=3.D0*sqrt(3.D0)*J3/(2.D0*J2**1.5D0)
  10 : c         if (abs(t).lt.-1.D0) t=-1.D0
  11 : c         if (abs(t).gt.1.D0)  t= 1.D0         
  12 :          t=-(0,1.D0)*log(t+sqrt(1.D0-t*t)*(0,1.D0))/3.D0
  13 :       endif
  14 :       return
  15 :       end
  16 :  
  17 :  
  18 :