SUBROUTINE DLSLDS(N,A,LDA,B,X) IMPLICIT DOUBLE PRECISION (A-H,O-Z) parameter(ldd2scr=100) common/dscrb/d2scrb(ldd2scr,ldd2scr) DIMENSION A(LDA,*),B(*),X(*) if(n.gt.ldd2scr) go to 100 DO 10 I=1,N DO 20 J=1,N d2scrb(I,J)=A(I,J) 20 CONTINUE 10 CONTINUE DO 30 I=1,N X(I)=B(I) 30 CONTINUE CALL DPOTRF('U',N,d2scrb,ldd2scr,INFO) CALL DPOTRS('U',N,1,d2scrb,ldd2scr,X,N,INFO) return 100 call err(n,'ldd2scr','d2scrblk') RETURN END