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
|
<HTML>
<HEAD><TITLE>MA02CD - SLICOT Library Routine Documentation</TITLE>
</HEAD>
<BODY>
<H2><A Name="MA02CD">MA02CD</A></H2>
<H3>
Pertransposing the central band of a square matrix
</H3>
<A HREF ="#Specification"><B>[Specification]</B></A>
<A HREF ="#Arguments"><B>[Arguments]</B></A>
<A HREF ="#Method"><B>[Method]</B></A>
<A HREF ="#References"><B>[References]</B></A>
<A HREF ="#Comments"><B>[Comments]</B></A>
<A HREF ="#Example"><B>[Example]</B></A>
<P>
<B><FONT SIZE="+1">Purpose</FONT></B>
<PRE>
To compute the pertranspose of a central band of a square matrix.
</PRE>
<A name="Specification"><B><FONT SIZE="+1">Specification</FONT></B></A>
<PRE>
SUBROUTINE MA02CD( N, KL, KU, A, LDA )
C .. Scalar Arguments ..
INTEGER KL, KU, LDA, N
C .. Array Arguments ..
DOUBLE PRECISION A(LDA,*)
</PRE>
<A name="Arguments"><B><FONT SIZE="+1">Arguments</FONT></B></A>
<P>
</PRE>
<B>Input/Output Parameters</B>
<PRE>
N (input) INTEGER
The order of the square matrix A. N >= 0.
KL (input) INTEGER
The number of subdiagonals of A to be pertransposed.
0 <= KL <= N-1.
KU (input) INTEGER
The number of superdiagonals of A to be pertransposed.
0 <= KU <= N-1.
A (input/output) DOUBLE PRECISION array, dimension (LDA,N)
On entry, the leading N-by-N part of this array must
contain a square matrix whose central band formed from
the KL subdiagonals, the main diagonal and the KU
superdiagonals will be pertransposed.
On exit, the leading N-by-N part of this array contains
the matrix A with its central band (the KL subdiagonals,
the main diagonal and the KU superdiagonals) pertransposed
(that is the elements of each antidiagonal appear in
reversed order). This is equivalent to forming P*B'*P,
where B is the matrix formed from the central band of A
and P is a permutation matrix with ones down the secondary
diagonal.
LDA INTEGER
The leading dimension of the array A. LDA >= max(1,N).
</PRE>
<A name="Comments"><B><FONT SIZE="+1">Further Comments</FONT></B></A>
<PRE>
None
</PRE>
<A name="Example"><B><FONT SIZE="+1">Example</FONT></B></A>
<P>
<B>Program Text</B>
<PRE>
None
</PRE>
<B>Program Data</B>
<PRE>
None
</PRE>
<B>Program Results</B>
<PRE>
None
</PRE>
<HR>
<A HREF=support.html><B>Return to Supporting Routines index</B></A></BODY>
</HTML>
|
