C CDDIIJ SOURCE CHAT 05/01/12 21:51:20 5004 SUBROUTINE C_DERDERINVARIIJJ (SIG,PIJKL,WIJKL,N) IMPLICIT INTEGER(I-N) integer n,i,j complex*16 sig(n),pijkl(6,6),wijkl(6,6),desv(6) do i=1,3 do j=1,3 pijkl(i,j)=-1.D0/3.D0 enddo pijkl(i,i)=2.D0/3.D0 enddo call C_Desviador(sig,desv,n) wijkl(1,1)=desv(1) wijkl(2,2)=desv(2) wijkl(3,3)=desv(3) wijkl(2,1)=desv(3) wijkl(3,1)=desv(2) wijkl(3,2)=desv(1) if (n.ge.4) then do i=4,n pijkl(i,i)=2.D0 enddo wijkl(4,1)=desv(4) wijkl(4,2)=desv(4) wijkl(4,3)=-2.D0*desv(4) wijkl(4,4)=-3.D0*desv(3) endif do i=1,n do j=1,i wijkl(i,j)=2.D0/3.D0*wijkl(i,j) enddo enddo do i=2,n do j=1,i-1 wijkl(j,i)=wijkl(i,j) enddo enddo return end