Sources Esope | Cast3M

cinvat
C CINVAT    SOURCE    CB215821  25/04/24    21:15:02     12248          
 
      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=CMPLX(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=-CMPLX(0.D0,1.D0)*log(t+sqrt(CMPLX(1.D0,0.D0)-t*t)*
     &    CMPLX(0.D0,1.D0))/CMPLX(3.D0,0.D0)
      endif
      return
      end