ZUNG2R(l) ) ZUNG2R(l)
NAME
ZUNG2R - generate an m by n complex matrix Q with orthonormal columns,
SYNOPSIS
SUBROUTINE ZUNG2R( M, N, K, A, LDA, TAU, WORK, INFO )
INTEGER INFO, K, LDA, M, N
COMPLEX*16 A( LDA, * ), TAU( * ), WORK( * )
PURPOSE
ZUNG2R generates an m by n complex matrix Q with orthonormal columns, which is defined as the first n columns of a product of k elementary
reflectors of order m
Q = H(1)H(2) . . . H(k)
as returned by ZGEQRF.
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*16 array, dimension (LDA,N)
On entry, the i-th column must contain the vector which defines the elementary reflector H(i), for i = 1,2,...,k, as returned by
ZGEQRF in the first 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*16 array, dimension (K)
TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by ZGEQRF.
WORK (workspace) COMPLEX*16 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 ZUNG2R(l)
Check Out this Related Man Page
ZUNGQR(l) ) ZUNGQR(l)
NAME
ZUNGQR - generate an M-by-N complex matrix Q with orthonormal columns,
SYNOPSIS
SUBROUTINE ZUNGQR( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
INTEGER INFO, K, LDA, LWORK, M, N
COMPLEX*16 A( LDA, * ), TAU( * ), WORK( * )
PURPOSE
ZUNGQR generates an M-by-N complex matrix Q with orthonormal columns, which is defined as the first N columns of a product of K elementary
reflectors of order M
Q = H(1)H(2) . . . H(k)
as returned by ZGEQRF.
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*16 array, dimension (LDA,N)
On entry, the i-th column must contain the vector which defines the elementary reflector H(i), for i = 1,2,...,k, as returned by
ZGEQRF in the first 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*16 array, dimension (K)
TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by ZGEQRF.
WORK (workspace/output) COMPLEX*16 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 ZUNGQR(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)
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)