Sources Esope | Cast3M

heflj2
C HEFLJ2    SOURCE    KK2000    14/04/09    21:15:26     8027
 
      SUBROUTINE HESSFLSIGJ2 (SIG,NDIMS,VAR,NDIMV,AMAT,NDIMX)
      IMPLICIT INTEGER(I-N)
      integer ndims,ndimv,ndimx,i,j
      real*8  sig(ndims),Amat(ndimx,ndimx),var(ndimv)
      real*8  x1(6),x2(6),x3(6),x22(6,6),x33(6,6),siginv(3)
      real*8  xqq,xtt,xqt,xpp,xq,xt,xp
      call InvariantesIJJ (sig,ndims,siginv)
      call DerInvariIJJ   (sig,x1,x2,x3,ndims)
      call DDINIJ(sig,x22,x33,ndims)
      call vflvminv(siginv,xp,xq,xt)
      call Dvflvminv(siginv,xpp,xqq,xqt,xtt)
      do i=1,ndims
       do j=1,ndims
          Amat(i,j)=xqq*x2(i)*x2(j)+xq*x22(i,j)+
     .              xtt*x3(i)*x3(j)+xt*x33(i,j)+
     .              xqt*(x2(i)*x3(j)+x2(j)*x3(i))+xpp*x1(i)*x1(j)
       enddo
      enddo
      return
      end