dlsets.f3
Usage
Functions/Subroutines
subroutine dlsets (m, p, n, a, af, lda, b, bf, ldb, c, cf, d, df, x, work, lwork, rwork, result)
DLSETS
Function/Subroutine Documentation
subroutine dlsets (integer m, integer p, integer n, double precision, dimension( lda, * ) a, double precision, dimension( lda, * ) af, integer lda, double precision, dimension( ldb, * ) b, double precision, dimension( ldb, * ) bf, integer ldb, double precision, dimension( * ) c, double precision, dimension( * ) cf, double precision, dimension( * ) d, double precision, dimension( * ) df, double precision, dimension( * ) x, double precision, dimension( lwork ) work, integer lwork, double precision, dimension( * ) rwork, double precision, dimension( 2 ) result)
DLSETS
Purpose:
!> !> DLSETS tests DGGLSE - a subroutine for solving linear equality !> constrained least square problem (LSE). !>
Parameters
M
!> M is INTEGER !> The number of rows of the matrix A. M >= 0. !>
P
!> P is INTEGER !> The number of rows of the matrix B. P >= 0. !>
N
!> N is INTEGER !> The number of columns of the matrices A and B. N >= 0. !>
A
!> A is DOUBLE PRECISION array, dimension (LDA,N) !> The M-by-N matrix A. !>
AF
!> AF is DOUBLE PRECISION array, dimension (LDA,N) !>
LDA
!> LDA is INTEGER !> The leading dimension of the arrays A, AF, Q and R. !> LDA >= max(M,N). !>
B
!> B is DOUBLE PRECISION array, dimension (LDB,N) !> The P-by-N matrix A. !>
BF
!> BF is DOUBLE PRECISION array, dimension (LDB,N) !>
LDB
!> LDB is INTEGER !> The leading dimension of the arrays B, BF, V and S. !> LDB >= max(P,N). !>
C
!> C is DOUBLE PRECISION array, dimension( M ) !> the vector C in the LSE problem. !>
CF
!> CF is DOUBLE PRECISION array, dimension( M ) !>
D
!> D is DOUBLE PRECISION array, dimension( P ) !> the vector D in the LSE problem. !>
DF
!> DF is DOUBLE PRECISION array, dimension( P ) !>
X
!> X is DOUBLE PRECISION array, dimension( N ) !> solution vector X in the LSE problem. !>
WORK
!> WORK is DOUBLE PRECISION array, dimension (LWORK) !>
LWORK
!> LWORK is INTEGER !> The dimension of the array WORK. !>
RWORK
!> RWORK is DOUBLE PRECISION array, dimension (M) !>
RESULT
!> RESULT is DOUBLE PRECISION array, dimension (2) !> The test ratios: !> RESULT(1) = norm( A*x - c )/ norm(A)*norm(X)*EPS !> RESULT(2) = norm( B*x - d )/ norm(B)*norm(X)*EPS !>
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
Definition at line 149 of file dlsets.f.