Esope source files

   1 : C CDIPQT    SOURCE    CHAT      05/01/12    21:51:31     5004           
   2 : 
   3 :       SUBROUTINE C_DERINVARIPQT (SIG,N,P_S,Q_S,T_S)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer n,i
   6 :       complex*16 sig(n),p_s(n),q_s(n),t_s(n),j2,j3,aux0,aux1,aux2
   6 :       call C_DerInvariIJJ(sig,p_s,q_s,t_s,n)
   7 :       do i=1,3
   8 :          p_s(i)=-p_s(i)/3.D0
   9 :       enddo
   9 :       call C_Invari_j2(sig,n,j2)
   9 :       call C_Invari_j3(sig,n,j3)
  10 :       aux0=0.D0
  11 :       aux1=0.D0
  12 :       aux2=0.D0
  13 :       if (abs(j2).ge.1.D-10) then
  14 :         aux0=sqrt(3.D0/j2)/2.D0
  15 :         aux1=aux0/j2
  16 :         aux2=1.5D0*j3/j2
  17 :       endif
  18 :       do i=1,n
  19 :          t_s(i)=aux1*(aux2*q_s(i)-t_s(i))
  20 :          q_s(i)=aux0*q_s(i)
  21 :       enddo
  22 :       end
  23 :  
  24 :  
  25 :