Esope source files

   1 : C CDDIIJ    SOURCE    CHAT      05/01/12    21:51:20     5004           
   2 : 
   3 :       SUBROUTINE C_DERDERINVARIIJJ (SIG,PIJKL,WIJKL,N)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer  n,i,j
   6 :       complex*16 sig(n),pijkl(6,6),wijkl(6,6),desv(6)
   7 :       do i=1,3
   8 :          do j=1,3
   9 :             pijkl(i,j)=-1.D0/3.D0
  10 :          enddo
  11 :          pijkl(i,i)=2.D0/3.D0
  12 :       enddo
  12 :       call C_Desviador(sig,desv,n)
  13 :       wijkl(1,1)=desv(1)
  14 :       wijkl(2,2)=desv(2)
  15 :       wijkl(3,3)=desv(3)
  16 :       wijkl(2,1)=desv(3)
  17 :       wijkl(3,1)=desv(2)
  18 :       wijkl(3,2)=desv(1)
  19 :       if (n.ge.4) then
  20 :        do i=4,n
  21 :          pijkl(i,i)=2.D0
  22 :        enddo
  23 :        wijkl(4,1)=desv(4)
  24 :        wijkl(4,2)=desv(4)
  25 :        wijkl(4,3)=-2.D0*desv(4)
  26 :        wijkl(4,4)=-3.D0*desv(3)
  27 :       endif
  28 :       do i=1,n
  29 :          do j=1,i
  30 :             wijkl(i,j)=2.D0/3.D0*wijkl(i,j)
  31 :          enddo
  32 :       enddo
  33 :       do i=2,n
  34 :          do j=1,i-1
  35 :             wijkl(j,i)=wijkl(i,j)
  36 :          enddo
  37 :       enddo
  38 :       return
  39 :       end
  40 :  
  41 :  
  42 :