12.07.2015 Views

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

C. Programas 141w(ll,j-ndumb)=e(j)sigw(j-ndumb)=sige(j)do i=1,natdLx(1,i,j-ndumb)=FFcc(3*i-2,j)dLx(2,i,j-ndumb)=FFcc(3*i-1,j)dLx(3,i,j-ndumb)=FFcc(3*i,j)enddoenddoCif (vlts.eqv..true.) thenC w(ll,ndumb+1)=-dtemporw(ll,1)=-dtemporendifC write(6,*) ’w’,(w(ll,j),j=1,9)CC l1=1C l2=3C33 do i=1,natC if(i.eq.1) thenCC write(6,5) l1,l2CC write(6,10) (sigw(j)*w(ll,j),j=l1,l2)CC write(6,*)C endifCC write(6,20) ((dLx(ic,i,j),ic=1,3),j=l1,l2)C enddoC l1=l2+1C l2=l2+3C if(l1.le.n) thenC if(l2.le.n) thenC goto 33C else if(l2.gt.n) thenC l2=nC goto 33C elseC continueC endifC endif5 format(//,1x,’w(’,i3,’) to w(’,i3,’):’,/)10 format(2x,3(10x,f8.1,5x))20 format(4x,3(2x,3(1x,f6.3)))C copy copy copyCreturnstopendC======================================================================SUBROUTINE HQRII(A,N,M,E,V)IMPLICIT DOUBLE PRECISION (A-H,O-Z)inclu<strong>de</strong> ’data.h’REAL*8 A((3*natoms)**2), E(3*natoms), V(3*natoms,3*natoms)C REAL*8 A((300)**2), E(300), V(300,300)*************************************************************** HQRII IS A DIAGONALISATION ROUTINE, WRITTEN BY YOSHITAKA BEPPU OF

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!