Sources Esope | Cast3M

heflp1
C HEFLP1    SOURCE    KK2000    14/04/09    21:15:26     8027
 
      SUBROUTINE HESSFLPOWDER1 (SIG,VAR,AMAT,NDIMA)
      IMPLICIT INTEGER(I-N)
      integer ndima,i,j
      real*8  sig(3),var,amat(ndima,ndima)
      real*8  aux2,x22(6,6),x33(6,6)
      real*8               nn1,nn2,eta0,aa1,aa2
      common /poder1_data/ nn1,nn2,eta0,aa1,aa2
      call DDINIJ(sig,x22,x33,3)
      aux2 = aa1/9.D0*2.D0
      do i=1,3
       do j=1,3
        amat(i,j)=2.D0*x22(i,j)+aux2
       enddo
      enddo
      return
      end