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 145CCIF(W(3,J).EQ.ZERO) W(3,J)=DELW(4,J+1)=0.D0GOTO 240230 IF(W(5,J).EQ.ZERO) W(5,J)=DELW(2,J)=-W(1,J)/W(5,J)240 W(4,J+1)=W(3,J)*W(2,J)+W(4,J+1)250 W(5,J+1)=W(4,J)*W(2,J)+W(5,J+1)IF(ABS(W(5,N)) .LT. EPS3) W(5,N)=DELDO 310 ITERE=1,5IF(ITERE.EQ.1) GOTO 270DO 260 J=1,NM1IF(W(3,J).EQ.ZERO) GOTO 260T=V(J,I)V(J,I)=V(J+1,I)V(J+1,I)=T260 V(J+1,I)=V(J,I)*W(2,J)+V(J+1,I)270 V(N,I)=V(N,I)/W(5,N)V(NM1,I)=(V(NM1,I)-V(N,I)*W(4,NM1))/W(5,NM1)VN=MAX(ABS(V(N,I)),ABS(V(NM1,I)),1.D-20)IF(N.EQ.2) GOTO 290K=NM2280 V(K,I)=(V(K,I)-V(K+1,I)*W(4,K)-V(K+2,I)*W(3,K))/W(5,K)VN=MAX(ABS(V(K,I)),VN,1.D-20)K=K-1IF(K.GE.1) GOTO 280290 S=EPS1/VNDO 300 J=1,N300 V(J,I)=V(J,I)*SIF(ITERE.GT.1 .AND. VN.GT.1) GOTO 320310 CONTINUETRANSFORMATION OF EIGENVECTORS320 IF(N.EQ.2) GOTO 360DO 350 J=1,NM2K=N-J-1IF(A((K*(K+1))/2).EQ.ZERO) GOTO 350KP1=K+1SUM=0.D0DO 330 KK=KP1,N330 SUM=SUM+A(K+(KK*(KK-1))/2)*V(KK,I)S=-SUM/A((K*(K+1))/2)DO 340 KK=KP1,N340 V(KK,I)=A(K+(KK*(KK-1))/2)*S+V(KK,I)350 CONTINUE360 DO 370 J=IG,IIF(ABS(E(J)-E(I)) .LT. EPS2) GOTO 380370 CONTINUEJ=I380 IG=JIF(IG .EQ. I) GOTO 410RE-ORTHOGONALISATIONDO 400 K=IG,IM1SUM=0.D0DO 390 J=1,N390 SUM=V(J,K)*V(J,I)+SUMS=-SUM

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

Saved successfully!

Ooh no, something went wrong!