File: clatzm.l

package info (click to toggle)
lapack 3.0.20000531a-28
  • links: PTS
  • area: main
  • in suites: sarge
  • size: 61,920 kB
  • ctags: 46,200
  • sloc: fortran: 584,835; perl: 8,226; makefile: 2,331; awk: 71; sh: 45
file content (105 lines) | stat: -rwxr-xr-x 2,212 bytes parent folder | download | duplicates (4)
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
100
101
102
103
104
105
.TH CLATZM l "15 June 2000" "LAPACK version 3.0" ")"
.SH NAME
CLATZM - routine is deprecated and has been replaced by routine CUNMRZ
.SH SYNOPSIS
.TP 19
SUBROUTINE CLATZM(
SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK )
.TP 19
.ti +4
CHARACTER
SIDE
.TP 19
.ti +4
INTEGER
INCV, LDC, M, N
.TP 19
.ti +4
COMPLEX
TAU
.TP 19
.ti +4
COMPLEX
C1( LDC, * ), C2( LDC, * ), V( * ), WORK( * )
.SH PURPOSE
This routine is deprecated and has been replaced by routine CUNMRZ. 
CLATZM applies a Householder matrix generated by CTZRQF to a matrix.

Let P = I - tau*u*u',   u = ( 1 ),
.br
                            ( v )
.br
where v is an (m-1) vector if SIDE = 'L', or a (n-1) vector if
SIDE = 'R'.
.br

If SIDE equals 'L', let
.br
       C = [ C1 ] 1
.br
           [ C2 ] m-1
.br
             n
.br
Then C is overwritten by P*C.
.br

If SIDE equals 'R', let
.br
       C = [ C1, C2 ] m
.br
              1  n-1
.br
Then C is overwritten by C*P.
.br

.SH ARGUMENTS
.TP 8
SIDE    (input) CHARACTER*1
= 'L': form P * C
.br
= 'R': form C * P
.TP 8
M       (input) INTEGER
The number of rows of the matrix C.
.TP 8
N       (input) INTEGER
The number of columns of the matrix C.
.TP 8
V       (input) COMPLEX array, dimension
(1 + (M-1)*abs(INCV)) if SIDE = 'L'
(1 + (N-1)*abs(INCV)) if SIDE = 'R'
The vector v in the representation of P. V is not used
if TAU = 0.
.TP 8
INCV    (input) INTEGER
The increment between elements of v. INCV <> 0
.TP 8
TAU     (input) COMPLEX
The value tau in the representation of P.
.TP 8
C1      (input/output) COMPLEX array, dimension
(LDC,N) if SIDE = 'L'
(M,1)   if SIDE = 'R'
On entry, the n-vector C1 if SIDE = 'L', or the m-vector C1
if SIDE = 'R'.

On exit, the first row of P*C if SIDE = 'L', or the first
column of C*P if SIDE = 'R'.
.TP 8
C2      (input/output) COMPLEX array, dimension
(LDC, N)   if SIDE = 'L'
(LDC, N-1) if SIDE = 'R'
On entry, the (m - 1) x n matrix C2 if SIDE = 'L', or the
m x (n - 1) matrix C2 if SIDE = 'R'.

On exit, rows 2:m of P*C if SIDE = 'L', or columns 2:m of C*P
if SIDE = 'R'.
.TP 8
LDC     (input) INTEGER
The leading dimension of the arrays C1 and C2.
LDC >= max(1,M).
.TP 8
WORK    (workspace) COMPLEX array, dimension
(N) if SIDE = 'L'
(M) if SIDE = 'R'