SUBROUTINE DLFTDS(N,A,LDA,FAC,LDFAC) IMPLICIT DOUBLE PRECISION(A-H,O-Z) DIMENSION A(LDA,*),FAC(LDFAC,*) CALL DCRGRG(N,A,LDA,FAC,LDFAC) CALL DPOTRF('U',N,FAC,LDFAC,INFO) RETURN END