
PROGRAM TEST

IMPLICIT NONE

C

INTEGER LENGTA, N

PARAMETER (N = 3)

PARAMETER (LENGTA = (N * N + N) / 2)

C

DOUBLE PRECISION ALPHA, AP(LENGTA), X(N)

INTEGER I

C

EXTERNAL DSPR

C

C Initialize the array A to store in symmetric packed form the

C matrix A shown below. Initialize the array X to store the

C vector x shown below.

C

C 100 200 300 1

C A = 200 400 500 x = 2

C 300 500 600 3

C

DATA AP / 100.0, 200.0, 300.0, 400.0, 500.0, 600.0 /

DATA X / 1.0, 2.0, 3.0 /

C

PRINT 1000

PRINT 1010, AP(1), AP(2), AP(3)

PRINT 1010, AP(2), AP(4), AP(5)

PRINT 1010, AP(3), AP(5), AP(6)

PRINT 1020

PRINT 1030, AP

PRINT 1040

PRINT 1050, (X(I), I = 1, N)

ALPHA = 1.0D0

CALL DSPR ('LOWER TRIANGULAR A', N, ALPHA, X, 1, AP)

PRINT 1060

PRINT 1010, AP(1), AP(2), AP(3)

PRINT 1010, AP(2), AP(4), AP(5)

PRINT 1010, AP(3), AP(5), AP(6)

C

1000 FORMAT (1X, 'A in full form:')

1010 FORMAT (1X, 3(2X, F5.1))

1020 FORMAT (/1X, 'A in packed form:')

1030 FORMAT (1X, 6(2X, F5.1))

1040 FORMAT (/1X, 'x:')

1050 FORMAT (1X, 2X, F5.1)

1060 FORMAT (/1X, 'Axx'':')

C

END

