CUNG2L(l) ) CUNG2L(l)
NAME
CUNG2L - generate an m by n complex matrix Q with orthonormal columns,
SYNOPSIS
SUBROUTINE CUNG2L( M, N, K, A, LDA, TAU, WORK, INFO )
INTEGER INFO, K, LDA, M, N
COMPLEX A( LDA, * ), TAU( * ), WORK( * )
PURPOSE
CUNG2L generates an m by n complex matrix Q with orthonormal columns, which is defined as the last n columns of a product of k elementary
reflectors of order m
Q = H(k) . . . H(2)H(1)
as returned by CGEQLF.
ARGUMENTS
M (input) INTEGER
The number of rows of the matrix Q. M >= 0.
N (input) INTEGER
The number of columns of the matrix Q. M >= N >= 0.
K (input) INTEGER
The number of elementary reflectors whose product defines the matrix Q. N >= K >= 0.
A (input/output) COMPLEX array, dimension (LDA,N)
On entry, the (n-k+i)-th column must contain the vector which defines the elementary reflector H(i), for i = 1,2,...,k, as returned
by CGEQLF in the last k columns of its array argument A. On exit, the m-by-n matrix Q.
LDA (input) INTEGER
The first dimension of the array A. LDA >= max(1,M).
TAU (input) COMPLEX array, dimension (K)
TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by CGEQLF.
WORK (workspace) COMPLEX array, dimension (N)
INFO (output) INTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument has an illegal value
LAPACK version 3.0 15 June 2000 CUNG2L(l)
Check Out this Related Man Page
CUNGQL(l) ) CUNGQL(l)
NAME
CUNGQL - generate an M-by-N complex matrix Q with orthonormal columns,
SYNOPSIS
SUBROUTINE CUNGQL( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
INTEGER INFO, K, LDA, LWORK, M, N
COMPLEX A( LDA, * ), TAU( * ), WORK( * )
PURPOSE
CUNGQL generates an M-by-N complex matrix Q with orthonormal columns, which is defined as the last N columns of a product of K elementary
reflectors of order M
Q = H(k) . . . H(2)H(1)
as returned by CGEQLF.
ARGUMENTS
M (input) INTEGER
The number of rows of the matrix Q. M >= 0.
N (input) INTEGER
The number of columns of the matrix Q. M >= N >= 0.
K (input) INTEGER
The number of elementary reflectors whose product defines the matrix Q. N >= K >= 0.
A (input/output) COMPLEX array, dimension (LDA,N)
On entry, the (n-k+i)-th column must contain the vector which defines the elementary reflector H(i), for i = 1,2,...,k, as returned
by CGEQLF in the last k columns of its array argument A. On exit, the M-by-N matrix Q.
LDA (input) INTEGER
The first dimension of the array A. LDA >= max(1,M).
TAU (input) COMPLEX array, dimension (K)
TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by CGEQLF.
WORK (workspace/output) COMPLEX array, dimension (LWORK)
On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
LWORK (input) INTEGER
The dimension of the array WORK. LWORK >= max(1,N). For optimum performance LWORK >= N*NB, where NB is the optimal blocksize.
If LWORK = -1, then a workspace query is assumed; the routine only calculates the optimal size of the WORK array, returns this
value as the first entry of the WORK array, and no error message related to LWORK is issued by XERBLA.
INFO (output) INTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument has an illegal value
LAPACK version 3.0 15 June 2000 CUNGQL(l)
Hi,
I'm having a problem printing two consecutive columns, as I iterate through a large matrix by twenty columns and I was looking for a solution.
My input file looks something like this
1 id1 A1 A2 A3 A4 A5 A6....A20 A21 A22 A23....A4001 A4002
2 id2 B1 B2 B3 B4 B5 B6...
3 id3 ...
4 id4... (8 Replies)
Hi All,
I have a huge matrix file consisting some some millions rows and 6000 columns. The contents are just floating point numbers in the matrix. I want to extract each column (i.e. 6000 of them) and store each column in a separate file.
For example, 1.dat will consist of elements from column... (4 Replies)
hello all,
i have an input file like that
A A X0
A B X1
A C X2
...
A Z Xx
B A X1
B B X3
....
Z A Xx
Z B X4
and i want to have an output like that
A B C D
A X0 X1 X2 Xy
B X1 X3 X4 (4 Replies)
This is my first post, I apologize if I have broken rules.
Some assistance with the following will be very helpful.
I have a couple of files, both should ultimately have common columns only, arranged in the same order.
This file needs to be transposed, to bring the rows to columns
... (2 Replies)
The following code transform the matrix to columns. Is it possible to do it other way around ( get the input from the output) ?
input
y1 y2 y3 y4 y5
x1 0.3 0.5 2.3 3.1 5.1
x2 1.2 4.1 3.5 1.7 1.2
x3 3.1 2.1 1.0 4.1 2.1
x4 5.0 4.0 6.0 7.0 1.1
output
x1 y1 0.3
x2 y1 1.2
x3... (1 Reply)