SORG2L(l) ) SORG2L(l)
NAME
SORG2L - generate an m by n real matrix Q with orthonormal columns,
SYNOPSIS
SUBROUTINE SORG2L( M, N, K, A, LDA, TAU, WORK, INFO )
INTEGER INFO, K, LDA, M, N
REAL A( LDA, * ), TAU( * ), WORK( * )
PURPOSE
SORG2L generates an m by n real 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 SGEQLF.
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) REAL 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 SGEQLF 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) REAL array, dimension (K)
TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by SGEQLF.
WORK (workspace) REAL 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 SORG2L(l)
Check Out this Related Man Page
SORGQL(l) ) SORGQL(l)
NAME
SORGQL - generate an M-by-N real matrix Q with orthonormal columns,
SYNOPSIS
SUBROUTINE SORGQL( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
INTEGER INFO, K, LDA, LWORK, M, N
REAL A( LDA, * ), TAU( * ), WORK( * )
PURPOSE
SORGQL generates an M-by-N real 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 SGEQLF.
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) REAL 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 SGEQLF 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) REAL array, dimension (K)
TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by SGEQLF.
WORK (workspace/output) REAL 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 SORGQL(l)
Hello all,
I have a text file that is arranged:
name 3 7 2 9 5
jim a d e g k
max d g u x g
rob f w v k o
This is just an example as my real file has >1000 individuals and >64,000 columns. I need to rearrange the file so that the columns appear in numerical order so that
name... (3 Replies)
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)