Esope source files

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