C C ------------------------------------------------------------------ C SUBROUTINE SMVPC(NBLOCK, BET, MAXJ, J, S, NUMBER, RESNRM, * ORTHCF, RV) C INTEGER NBLOCK, MAXJ, J, NUMBER REAL BET(NBLOCK,1), S(MAXJ,1), RESNRM(1), * ORTHCF(1), RV(1) C C THIS SUBROUTINE COMPUTES THE NORM AND THE SMALLEST ELEMENT C (IN ABSOLUTE VALUE) OF THE VECTOR BET*SJI, WHERE SJI C IS AN NBLOCK VECTOR OF THE LAST NBLOCK ELEMENTS OF THE ITH C EIGENVECTOR OF T. THESE QUANTITIES ARE THE RESIDUAL NORM C AND THE ORTHOGONALITY COEFFICIENT RESPECTIVELY FOR THE C CORRESPONDING RITZ PAIR. THE ORTHOGONALITY COEFFICIENT IS C NORMALIZED TO ACCOUNT FOR THE LOCAL REORTHOGONALIZATION. C INTEGER I, K, M REAL SDOT, SNRM2, ABS, AMIN1 C M = J - NBLOCK + 1 DO 20 I=1,NUMBER DO 10 K=1,NBLOCK RV(K) = SDOT(NBLOCK,S(M,I),1,BET(K,1),NBLOCK) IF (K.EQ.1) ORTHCF(I) = ABS(RV(K)) ORTHCF(I) = AMIN1(ORTHCF(I),ABS(RV(K))) 10 CONTINUE RESNRM(I) = SNRM2(NBLOCK,RV,1) 20 CONTINUE RETURN END .