zglmts.f3

Usage

Functions/Subroutines

subroutine zglmts (n, m, p, a, af, lda, b, bf, ldb, d, df, x, u, work, lwork, rwork, result)
ZGLMTS

Function/Subroutine Documentation

subroutine zglmts (integer n, integer m, integer p, complex16, dimension( lda, ) a, complex16, dimension( lda, ) af, integer lda, complex16, dimension( ldb, ) b, complex16, dimension( ldb, ) bf, integer ldb, complex16, dimension( ) d, complex16, dimension( ) df, complex16, dimension( ) x, complex16, dimension( ) u, complex16, dimension( lwork ) work, integer lwork, double precision, dimension( ) rwork, double precision result)

ZGLMTS

Purpose:

!>
!> ZGLMTS tests ZGGGLM - a subroutine for solving the generalized
!> linear model problem.
!>

Parameters

!>          N is INTEGER
!>          The number of rows of the matrices A and B.  N >= 0.
!>

!>          M is INTEGER
!>          The number of columns of the matrix A.  M >= 0.
!>

!>          P is INTEGER
!>          The number of columns of the matrix B.  P >= 0.
!>

!>          A is COMPLEX*16 array, dimension (LDA,M)
!>          The N-by-M matrix A.
!>

!>          AF is COMPLEX*16 array, dimension (LDA,M)
!>

LDA

!>          LDA is INTEGER
!>          The leading dimension of the arrays A, AF. LDA >= max(M,N).
!>

!>          B is COMPLEX*16 array, dimension (LDB,P)
!>          The N-by-P matrix A.
!>

!>          BF is COMPLEX*16 array, dimension (LDB,P)
!>

LDB

!>          LDB is INTEGER
!>          The leading dimension of the arrays B, BF. LDB >= max(P,N).
!>

!>          D is COMPLEX*16 array, dimension( N )
!>          On input, the left hand side of the GLM.
!>

!>          DF is COMPLEX*16 array, dimension( N )
!>

!>          X is COMPLEX*16 array, dimension( M )
!>          solution vector X in the GLM problem.
!>

!>          U is COMPLEX*16 array, dimension( P )
!>          solution vector U in the GLM 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
!>          The test ratio:
!>                           norm( d - A*x - B*u )
!>            RESULT = -----------------------------------------
!>                     (norm(A)+norm(B))*(norm(x)+norm(u))*EPS
!>

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 144 of file zglmts.f.