LAPACK  3.5.0
LAPACK: Linear Algebra PACKage
 All Classes Files Functions Variables Typedefs Macros
dglmts.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine dglmts (N, M, P, A, AF, LDA, B, BF, LDB, D, DF, X, U, WORK, LWORK, RWORK, RESULT)
 DGLMTS More...
 

Function/Subroutine Documentation

subroutine dglmts ( integer  N,
integer  M,
integer  P,
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( * )  D,
double precision, dimension( * )  DF,
double precision, dimension( * )  X,
double precision, dimension( * )  U,
double precision, dimension( lwork )  WORK,
integer  LWORK,
double precision, dimension( * )  RWORK,
double precision  RESULT 
)

DGLMTS

Purpose:
 DGLMTS tests DGGGLM - a subroutine for solving the generalized
 linear model problem.
Parameters
[in]N
          N is INTEGER
          The number of rows of the matrices A and B.  N >= 0.
[in]M
          M is INTEGER
          The number of columns of the matrix A.  M >= 0.
[in]P
          P is INTEGER
          The number of columns of the matrix B.  P >= 0.
[in]A
          A is DOUBLE PRECISION array, dimension (LDA,M)
          The N-by-M matrix A.
[out]AF
          AF is DOUBLE PRECISION array, dimension (LDA,M)
[in]LDA
          LDA is INTEGER
          The leading dimension of the arrays A, AF. LDA >= max(M,N).
[in]B
          B is DOUBLE PRECISION array, dimension (LDB,P)
          The N-by-P matrix A.
[out]BF
          BF is DOUBLE PRECISION array, dimension (LDB,P)
[in]LDB
          LDB is INTEGER
          The leading dimension of the arrays B, BF. LDB >= max(P,N).
[in]D
          D is DOUBLE PRECISION array, dimension( N )
          On input, the left hand side of the GLM.
[out]DF
          DF is DOUBLE PRECISION array, dimension( N )
[out]X
          X is DOUBLE PRECISION array, dimension( M )
          solution vector X in the GLM problem.
[out]U
          U is DOUBLE PRECISION array, dimension( P )
          solution vector U in the GLM problem.
[out]WORK
          WORK is DOUBLE PRECISION array, dimension (LWORK)
[in]LWORK
          LWORK is INTEGER
          The dimension of the array WORK.
[out]RWORK
          RWORK is DOUBLE PRECISION array, dimension (M)
[out]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.
Date
November 2011

Definition at line 146 of file dglmts.f.

Here is the call graph for this function:

Here is the caller graph for this function: