SUBROUTINE DEVLSF(N,A,LDA,EVAL) IMPLICIT DOUBLE PRECISION(A-H,O-Z) PARAMETER(LDD1SCR=10000) COMMON/DSCRB/D1SCRB(LDD1SCR) PARAMETER(LDD2SCR=100) COMMON/DSCRB/D2SCRB(LDD2SCR,LDD2SCR) DIMENSION A(LDA,*), EVAL(*),info(100) IF(N.GT.LDD1SCR) GO TO 100 IF(N.GT.LDD2SCR) GO TO 110 CALL DCRGRG(N,A,LDA,D2SCRB,LDD2SCR) CALL DSYEV('N','U',N,D2SCRB,LDD2SCR,EVAL,D1SCRB,LDD1SCR,INFO) RETURN 100 CALL ERR(N,'LDD1SCR','D1SCRBLK') RETURN 110 CALL ERR(N,'LDD2SCR','D2SCRBLK') RETURN END