SUBROUTINE TSLD(A,X,R,M) INTEGER M DOUBLE PRECISION A(1),X(M),R(1) C C TSLD CALLS TSLD1 TO SOLVE THE DOUBLE PRECISION LINEAR SYSTEM C A * X = B C WITH THE T - MATRIX A . C C ON ENTRY C C A DOUBLE PRECISION(2*M - 1) C THE FIRST ROW OF THE T - MATRIX FOLLOWED BY ITS C FIRST COLUMN BEGINNING WITH THE SECOND ELEMENT . C ON RETURN A IS UNALTERED . C C X DOUBLE PRECISION(M) C THE RIGHT HAND SIDE VECTOR B . C C R DOUBLE PRECISION(2*M - 2) C A WORK VECTOR . C C M INTEGER C THE ORDER OF THE MATRIX A . C C ON RETURN C C X THE SOLUTION VECTOR . C C TOEPLITZ PACKAGE. THIS VERSION DATED 07/23/82 . C C SUBROUTINES AND FUNCTIONS C C TOEPLITZ PACKAGE ... TSLD1 C C CALL SUBROUTINE TSLD1 C CALL TSLD1(A,A(M+1),X,X,R,R(M),M) C RETURN END