Sources Esope | Cast3M

cinvj2
C CINVJ2    SOURCE    KK2000    14/04/10    21:15:03     8031
 
      SUBROUTINE C_INVARI_J2 (SIG,N,J2)
      IMPLICIT INTEGER(I-N)
      integer n,i
      complex*16 sig(n),J2,desv(6)
      call CDSVIA(sig,desv,n)
      J2=(0.D0,0.D0)
      do i=1,3
         J2=J2+desv(i)*desv(i)
      enddo
      J2=J2/(2.D0,0.D0)
      if (n.lt.4) return
      do i=4,n
         J2=J2+desv(i)*desv(i)
      enddo
      return
      end