1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
SUBROUTINE MA02ED( UPLO, N, A, LDA )
C
C SLICOT RELEASE 5.0.
C
C Copyright (c) 2002-2009 NICONET e.V.
C
C This program is free software: you can redistribute it and/or
C modify it under the terms of the GNU General Public License as
C published by the Free Software Foundation, either version 2 of
C the License, or (at your option) any later version.
C
C This program is distributed in the hope that it will be useful,
C but WITHOUT ANY WARRANTY; without even the implied warranty of
C MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
C GNU General Public License for more details.
C
C You should have received a copy of the GNU General Public License
C along with this program. If not, see
C <http://www.gnu.org/licenses/>.
C
C PURPOSE
C
C To store by symmetry the upper or lower triangle of a symmetric
C matrix, given the other triangle.
C
C ARGUMENTS
C
C Mode Parameters
C
C UPLO CHARACTER*1
C Specifies which part of the matrix is given as follows:
C = 'U': Upper triangular part;
C = 'L': Lower triangular part.
C For all other values, the array A is not referenced.
C
C Input/Output Parameters
C
C N (input) INTEGER
C The order of the matrix A. N >= 0.
C
C A (input/output) DOUBLE PRECISION array, dimension (LDA,N)
C On entry, the leading N-by-N upper triangular part
C (if UPLO = 'U'), or lower triangular part (if UPLO = 'L'),
C of this array must contain the corresponding upper or
C lower triangle of the symmetric matrix A.
C On exit, the leading N-by-N part of this array contains
C the symmetric matrix A with all elements stored.
C
C LDA INTEGER
C The leading dimension of the array A. LDA >= max(1,N).
C
C CONTRIBUTOR
C
C V. Sima, Research Institute for Informatics, Bucharest, Romania,
C Oct. 1998.
C
C REVISIONS
C
C -
C
C ******************************************************************
C
C .. Scalar Arguments ..
CHARACTER UPLO
INTEGER LDA, N
C .. Array Arguments ..
DOUBLE PRECISION A(LDA,*)
C .. Local Scalars ..
INTEGER J
C .. External Functions ..
LOGICAL LSAME
EXTERNAL LSAME
C .. External Subroutines ..
EXTERNAL DCOPY
C
C .. Executable Statements ..
C
C For efficiency reasons, the parameters are not checked for errors.
C
IF( LSAME( UPLO, 'L' ) ) THEN
C
C Construct the upper triangle of A.
C
DO 20 J = 2, N
CALL DCOPY( J-1, A(J,1), LDA, A(1,J), 1 )
20 CONTINUE
C
ELSE IF( LSAME( UPLO, 'U' ) ) THEN
C
C Construct the lower triangle of A.
C
DO 40 J = 2, N
CALL DCOPY( J-1, A(1,J), 1, A(J,1), LDA )
40 CONTINUE
C
END IF
RETURN
C *** Last line of MA02ED ***
END
|
