Sources Esope | Cast3M

cdipqt
C CDIPQT    SOURCE    CB215821  25/04/24    21:15:02     12248          
 
      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)/CMPLX(3.D0,0.D0)
      enddo
      call C_Invari_j2(sig,n,j2)
      call C_Invari_j3(sig,n,j3)
      aux0=CMPLX(0.D0,0.D0)
      aux1=CMPLX(0.D0,0.D0)
      aux2=CMPLX(0.D0,0.D0)
      if (abs(j2).ge.1.D-10) then
        aux0=sqrt(CMPLX(3.D0,0.D0)/j2)/CMPLX(2.D0,0.D0)
        aux1=aux0/j2
        aux2=CMPLX(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