zlsets.f3
Usage
Functions/Subroutines
subroutine zlsets (m, p, n, a, af, lda, b, bf, ldb, c, cf, d, df, x, work, lwork, rwork, result)
ZLSETS
Function/Subroutine Documentation
subroutine zlsets (integer m, integer p, integer n, complex*16, dimension( lda, * ) a, complex*16, dimension( lda, * ) af, integer lda, complex*16, dimension( ldb, * ) b, complex*16, dimension( ldb, * ) bf, integer ldb, complex*16, dimension( * ) c, complex*16, dimension( * ) cf, complex*16, dimension( * ) d, complex*16, dimension( * ) df, complex*16, dimension( * ) x, complex*16, dimension( lwork ) work, integer lwork, double precision, dimension( * ) rwork, double precision, dimension( 2 ) result)
ZLSETS
Purpose:
!> !> ZLSETS tests ZGGLSE - 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 COMPLEX*16 array, dimension (LDA,N) !> The M-by-N matrix A. !>
AF
!> AF is COMPLEX*16 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 COMPLEX*16 array, dimension (LDB,N) !> The P-by-N matrix A. !>
BF
!> BF is COMPLEX*16 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 COMPLEX*16 array, dimension( M ) !> the vector C in the LSE problem. !>
CF
!> CF is COMPLEX*16 array, dimension( M ) !>
D
!> D is COMPLEX*16 array, dimension( P ) !> the vector D in the LSE problem. !>
DF
!> DF is COMPLEX*16 array, dimension( P ) !>
X
!> X is COMPLEX*16 array, dimension( N ) !> solution vector X in the LSE problem. !>
WORK
!> WORK is COMPLEX*16 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 zlsets.f.