Click or drag to resize

DLAQP2Run Method

Purpose ======= DLAQP2 computes a QR factorization with column pivoting of the block A(OFFSET+1:M,1:N). The block A(1:OFFSET,1:N) is accordingly pivoted, but not factorized.

Namespace: DotNumerics.LinearAlgebra.CSLapack
Assembly: DWSIM.MathOps.DotNumerics (in DWSIM.MathOps.DotNumerics.dll) Version: 1.0.0.0 (1.0.0.0)
Syntax
public void Run(
	int M,
	int N,
	int OFFSET,
	ref double[] A,
	int offset_a,
	int LDA,
	ref int[] JPVT,
	int offset_jpvt,
	ref double[] TAU,
	int offset_tau,
	ref double[] VN1,
	int offset_vn1,
	ref double[] VN2,
	int offset_vn2,
	ref double[] WORK,
	int offset_work
)
Request Example View Source

Parameters

M  Int32
(input) INTEGER The number of rows of the matrix A. M .GE. 0.
N  Int32
(input) INTEGER The number of columns of the matrix A. N .GE. 0.
OFFSET  Int32
(input) INTEGER The number of rows of the matrix A that must be pivoted but no factorized. OFFSET .GE. 0.
A  Double
(input/output) DOUBLE PRECISION array, dimension (LDA,N) On entry, the M-by-N matrix A. On exit, the upper triangle of block A(OFFSET+1:M,1:N) is the triangular factor obtained; the elements in block A(OFFSET+1:M,1:N) below the diagonal, together with the array TAU, represent the orthogonal matrix Q as a product of elementary reflectors. Block A(1:OFFSET,1:N) has been accordingly pivoted, but no factorized.
offset_a  Int32
 
LDA  Int32
(input) INTEGER The leading dimension of the array A. LDA .GE. max(1,M).
JPVT  Int32
(input/output) INTEGER array, dimension (N) On entry, if JPVT(i) .ne. 0, the i-th column of A is permuted to the front of A*P (a leading column); if JPVT(i) = 0, the i-th column of A is a free column. On exit, if JPVT(i) = k, then the i-th column of A*P was the k-th column of A.
offset_jpvt  Int32
 
TAU  Double
(output) DOUBLE PRECISION array, dimension (min(M,N)) The scalar factors of the elementary reflectors.
offset_tau  Int32
 
VN1  Double
(input/output) DOUBLE PRECISION array, dimension (N) The vector with the partial column norms.
offset_vn1  Int32
 
VN2  Double
(input/output) DOUBLE PRECISION array, dimension (N) The vector with the exact column norms.
offset_vn2  Int32
 
WORK  Double
(workspace) DOUBLE PRECISION array, dimension (N)
offset_work  Int32
 
See Also