Esope source files

   1 : C HEFLRH    SOURCE    CHAT      05/01/13    00:23:54     5004           
   2 : 
   3 :       SUBROUTINE HESSFLSIGRHMC (SIG,NDIMS,VAR,NDIMV,AMAT,NDIMX)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       integer ndims,ndimv,ndimx,i,j
   6 :       real*8  sig(ndims),Amat(ndimx,ndimx),var(ndimv)
   7 :       real*8  x1(6),x2(6),x3(6),x22(6,6),x33(6,6),siginv(3)
   8 :       real*8  xqq,xtt,xqt,xpp,xq,xt,xp
   8 :       call InvariantesIJJ (sig,ndims,siginv)
   8 :       call DerInvariIJJ   (sig,x1,x2,x3,ndims)
   8 :       call DerDerInvariIJJ(sig,x22,x33,ndims)
   8 :       call vflrhmcinv(siginv,xp,xq,xt)
   8 :       call Dvflrhmcinv(siginv,xpp,xqq,xqt,xtt)
   9 :       do i=1,ndims
  10 :        do j=1,ndims
  11 :           Amat(i,j)=xqq*x2(i)*x2(j)+xq*x22(i,j)+
  12 :      .              xtt*x3(i)*x3(j)+xt*x33(i,j)+
  13 :      .              xqt*(x2(i)*x3(j)+x2(j)*x3(i))+xpp*x1(i)*x1(j)
  14 :        enddo
  15 :       enddo
  16 :       return
  17 :       end
  18 :  
  19 :  
  20 :