Purpose
=======
DORMTR 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 of order nq, with nq = m if
SIDE = 'L' and nq = n if SIDE = 'R'. Q is defined as the product of
nq-1 elementary reflectors, as returned by DSYTRD:
if UPLO = 'U', Q = H(nq-1) . . . H(2) H(1);
if UPLO = 'L', Q = H(1) H(2) . . . H(nq-1).
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 UPLO,
string TRANS,
int M,
int N,
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,
UPLO As String,
TRANS As String,
M As Integer,
N 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'
- UPLO String
-
(input) CHARACTER*1
= 'U': Upper triangle of A contains elementary reflectors
from DSYTRD;
= 'L': Lower triangle of A contains elementary reflectors
from DSYTRD.
- 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.
- A Double
-
(input) DOUBLE PRECISION array, dimension
(LDA,M) if SIDE = 'L'
(LDA,N) if SIDE = 'R'
The vectors which define the elementary reflectors, as
returned by DSYTRD.
- offset_a Int32
-
- LDA Int32
-
(input) INTEGER
The leading dimension of the array A.
LDA .GE. max(1,M) if SIDE = 'L'; LDA .GE. max(1,N) if SIDE = 'R'.
- TAU Double
-
(input) DOUBLE PRECISION array, dimension
(M-1) if SIDE = 'L'
(N-1) if SIDE = 'R'
TAU(i) must contain the scalar factor of the elementary
reflector H(i), as returned by DSYTRD.
- 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