sglmts.f3

Usage

Functions/Subroutines

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

Function/Subroutine Documentation

subroutine sglmts (integer n, integer m, integer p, real, dimension( lda, * ) a, real, dimension( lda, * ) af, integer lda, real, dimension( ldb, * ) b, real, dimension( ldb, * ) bf, integer ldb, real, dimension( * ) d, real, dimension( * ) df, real, dimension( * ) x, real, dimension( * ) u, real, dimension( lwork ) work, integer lwork, real, dimension( * ) rwork, real result)

SGLMTS

Purpose:

!>
!> SGLMTS tests SGGGLM - 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 REAL array, dimension (LDA,M)
!>          The N-by-M matrix A.
!>

!>          AF is REAL array, dimension (LDA,M)
!>

LDA

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

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

!>          BF is REAL array, dimension (LDB,P)
!>

LDB

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

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

!>          DF is REAL array, dimension( N )
!>

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

!>          U is REAL array, dimension( P )
!>          solution vector U in the GLM problem.
!>

WORK

!>          WORK is REAL array, dimension (LWORK)
!>

LWORK

!>          LWORK is INTEGER
!>          The dimension of the array WORK.
!>

RWORK

!>          RWORK is REAL array, dimension (M)
!>

RESULT

!>          RESULT is REAL
!>          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 147 of file sglmts.f.