1 : C HEFLP1 SOURCE CHAT 05/01/13 00:23:45 5004 2 : 3 : SUBROUTINE HESSFLPOWDER1 (SIG,VAR,AMAT,NDIMA) 4 : IMPLICIT INTEGER(I-N) 5 : integer ndima,m,i,j 6 : real*8 sig(3),var,amat(ndima,ndima) 7 : real*8 vecm(3),aux1,aux2,xj2,x22(6,6),x33(6,6) 8 : real*8 nn1,nn2,eta0,aa1,aa2 9 : common /poder1_data/ nn1,nn2,eta0,aa1,aa2 9 : call DerDerInvariIJJ(sig,x22,x33,3) 10 : aux2 = aa1/9.D0*2.D0 11 : do i=1,3 12 : do j=1,3 13 : amat(i,j)=2.D0*x22(i,j)+aux2 14 : enddo 15 : enddo 16 : return 17 : end 18 : 19 : 20 :
© Cast3M 2003 - All rights reserved.
Disclaimer