Sources Esope | Cast3M

cdipqt
C CDIPQT    SOURCE    KK2000    14/04/09    21:15:00     8027
 
      SUBROUTINE C_DERINVARIPQT (SIG,N,P_S,Q_S,T_S)
      IMPLICIT INTEGER(I-N)
      integer n,i
      complex*16 sig(n),p_s(n),q_s(n),t_s(n),j2,j3,aux0,aux1,aux2
      call C_DerInvariIJJ(sig,p_s,q_s,t_s,n)
      do i=1,3
         p_s(i)=-p_s(i)/(3.D0,0.D0)
      enddo
      call C_Invari_j2(sig,n,j2)
      call C_Invari_j3(sig,n,j3)
      aux0=(0.D0,0.D0)
      aux1=(0.D0,0.D0)
      aux2=(0.D0,0.D0)
      if (abs(j2).ge.1.D-10) then
        aux0=sqrt((3.D0,0.D0)/j2)/(2.D0,0.D0)
        aux1=aux0/j2
        aux2=(1.5D0,0.D0)*j3/j2
      endif
      do i=1,n
         t_s(i)=aux1*(aux2*q_s(i)-t_s(i))
         q_s(i)=aux0*q_s(i)
      enddo
      end