1DSYGST(1)                LAPACK routine (version 3.2)                DSYGST(1)
2
3
4

NAME

6       DSYGST  - reduces a real symmetric-definite generalized eigenproblem to
7       standard form
8

SYNOPSIS

10       SUBROUTINE DSYGST( ITYPE, UPLO, N, A, LDA, B, LDB, INFO )
11
12           CHARACTER      UPLO
13
14           INTEGER        INFO, ITYPE, LDA, LDB, N
15
16           DOUBLE         PRECISION A( LDA, * ), B( LDB, * )
17

PURPOSE

19       DSYGST reduces a real symmetric-definite  generalized  eigenproblem  to
20       standard form.  If ITYPE = 1, the problem is A*x = lambda*B*x,
21       and  A  is  overwritten  by inv(U**T)*A*inv(U) or inv(L)*A*inv(L**T) If
22       ITYPE = 2 or 3, the problem is A*B*x = lambda*x or
23       B*A*x = lambda*x, and A is overwritten by U*A*U**T or L**T*A*L.  B must
24       have been previously factorized as U**T*U or L*L**T by DPOTRF.
25

ARGUMENTS

27       ITYPE   (input) INTEGER
28               = 1: compute inv(U**T)*A*inv(U) or inv(L)*A*inv(L**T);
29               = 2 or 3: compute U*A*U**T or L**T*A*L.
30
31       UPLO    (input) CHARACTER*1
32               =  'U':   Upper  triangle  of  A is stored and B is factored as
33               U**T*U; = 'L':  Lower triangle of A is stored and B is factored
34               as L*L**T.
35
36       N       (input) INTEGER
37               The order of the matrices A and B.  N >= 0.
38
39       A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
40               On  entry,  the symmetric matrix A.  If UPLO = 'U', the leading
41               N-by-N upper triangular part of A contains the upper triangular
42               part of the matrix A, and the strictly lower triangular part of
43               A is not referenced.  If UPLO = 'L', the leading  N-by-N  lower
44               triangular  part of A contains the lower triangular part of the
45               matrix A, and the strictly upper triangular part of  A  is  not
46               referenced.   On  exit,  if  INFO  = 0, the transformed matrix,
47               stored in the same format as A.
48
49       LDA     (input) INTEGER
50               The leading dimension of the array A.  LDA >= max(1,N).
51
52       B       (input) DOUBLE PRECISION array, dimension (LDB,N)
53               The triangular factor from the Cholesky factorization of B,  as
54               returned by DPOTRF.
55
56       LDB     (input) INTEGER
57               The leading dimension of the array B.  LDB >= max(1,N).
58
59       INFO    (output) INTEGER
60               = 0:  successful exit
61               < 0:  if INFO = -i, the i-th argument had an illegal value
62
63
64
65 LAPACK routine (version 3.2)    November 2008                       DSYGST(1)
Impressum