Purpose
=======
DORMQL overwrites the general real M-by-N matrix C with
SIDE = 'L' SIDE = 'R'
TRANS = 'N': Q * C C * Q
TRANS = 'T': Q**T * C C * Q**T
where Q is a real orthogonal matrix defined as the product of k
elementary reflectors
Q = H(k) . . . H(2) H(1)
as returned by DGEQLF. Q is of order M if SIDE = 'L' and of order N
if SIDE = 'R'.
Namespace: DotNumerics.LinearAlgebra.CSLapackAssembly: DWSIM.MathOps.DotNumerics (in DWSIM.MathOps.DotNumerics.dll) Version: 1.0.0.0 (1.0.0.0)
Syntax public void Run(
string SIDE,
string TRANS,
int M,
int N,
int K,
ref double[] A,
int offset_a,
int LDA,
double[] TAU,
int offset_tau,
ref double[] C,
int offset_c,
int LDC,
ref double[] WORK,
int offset_work,
int LWORK,
ref int INFO
)
Public Sub Run (
SIDE As String,
TRANS As String,
M As Integer,
N As Integer,
K As Integer,
ByRef A As Double(),
offset_a As Integer,
LDA As Integer,
TAU As Double(),
offset_tau As Integer,
ByRef C As Double(),
offset_c As Integer,
LDC As Integer,
ByRef WORK As Double(),
offset_work As Integer,
LWORK As Integer,
ByRef INFO As Integer
)
Request Example
View SourceParameters
- SIDE String
-
= 'L' SIDE = 'R'
- TRANS String
-
(input) CHARACTER*1
= 'N': No transpose, apply Q;
= 'T': Transpose, apply Q**T.
- M Int32
-
(input) INTEGER
The number of rows of the matrix C. M .GE. 0.
- N Int32
-
(input) INTEGER
The number of columns of the matrix C. N .GE. 0.
- K Int32
-
(input) INTEGER
The number of elementary reflectors whose product defines
the matrix Q.
If SIDE = 'L', M .GE. K .GE. 0;
if SIDE = 'R', N .GE. K .GE. 0.
- A Double
-
(input) DOUBLE PRECISION array, dimension (LDA,K)
The i-th column must contain the vector which defines the
elementary reflector H(i), for i = 1,2,...,k, as returned by
DGEQLF in the last k columns of its array argument A.
A is modified by the routine but restored on exit.
- offset_a Int32
-
- LDA Int32
-
(input) INTEGER
The leading dimension of the array A.
If SIDE = 'L', LDA .GE. max(1,M);
if SIDE = 'R', LDA .GE. max(1,N).
- TAU Double
-
(input) DOUBLE PRECISION array, dimension (K)
TAU(i) must contain the scalar factor of the elementary
reflector H(i), as returned by DGEQLF.
- offset_tau Int32
-
- C Double
-
(input/output) DOUBLE PRECISION array, dimension (LDC,N)
On entry, the M-by-N matrix C.
On exit, C is overwritten by Q*C or Q**T*C or C*Q**T or C*Q.
- offset_c Int32
-
- LDC Int32
-
(input) INTEGER
The leading dimension of the array C. LDC .GE. max(1,M).
- WORK Double
-
(workspace/output) DOUBLE PRECISION array, dimension (MAX(1,LWORK))
On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
- offset_work Int32
-
- LWORK Int32
-
(input) INTEGER
The dimension of the array WORK.
If SIDE = 'L', LWORK .GE. max(1,N);
if SIDE = 'R', LWORK .GE. max(1,M).
For optimum performance LWORK .GE. N*NB if SIDE = 'L', and
LWORK .GE. M*NB if SIDE = 'R', 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 Int32
-
(output) INTEGER
= 0: successful exit
.LT. 0: if INFO = -i, the i-th argument had an illegal value
See Also