Sources Esope | Cast3M

cinvat
C CINVAT    SOURCE    KK2000    14/04/09    21:15:08     8027
 
      SUBROUTINE CINVAT(SIG,N,T)
      IMPLICIT INTEGER(I-N)
      integer n
      complex*16 sig(n),t,J2,J3
      call C_Invari_J2(sig,n,J2)
      t=(0.D0,0.D0)
      if (abs(J2).gt.1.D-17) then
         call C_Invari_J3(sig,n,J3)
         t=CMPLX(3.D0*sqrt(3.D0))*J3/(CMPLX(2.D0)*J2**CMPLX(1.5D0))
c         if (abs(t).lt.-1.D0) t=-1.D0
c         if (abs(t).gt.1.D0)  t= 1.D0
         t=-(0.D0,1.D0)*log(t+sqrt((1.D0,0.D0)-t*t)*(0.D0,1.D0))/
     &    (3.D0,0.D0)
      endif
      return
      end