hbgvx#

Functions

void chbgvx(
    const char*          jobz,
    const char*          range,
    const char*          uplo,
    const INT            n,
    const INT            ka,
    const INT            kb,
          c64*  restrict AB,
    const INT            ldab,
          c64*  restrict BB,
    const INT            ldbb,
          c64*  restrict Q,
    const INT            ldq,
    const f32            vl,
    const f32            vu,
    const INT            il,
    const INT            iu,
    const f32            abstol,
          INT*           m,
          f32*  restrict W,
          c64*  restrict Z,
    const INT            ldz,
          c64*  restrict work,
          f32*  restrict rwork,
          INT*  restrict iwork,
          INT*  restrict ifail,
          INT*           info
);
void chbgvx(const char *jobz, const char *range, const char *uplo, const INT n, const INT ka, const INT kb, c64 *restrict AB, const INT ldab, c64 *restrict BB, const INT ldbb, c64 *restrict Q, const INT ldq, const f32 vl, const f32 vu, const INT il, const INT iu, const f32 abstol, INT *m, f32 *restrict W, c64 *restrict Z, const INT ldz, c64 *restrict work, f32 *restrict rwork, INT *restrict iwork, INT *restrict ifail, INT *info)#

CHBGVX computes all the eigenvalues, and optionally, the eigenvectors of a complex generalized Hermitian-definite banded eigenproblem, of the form A*x=(lambda)*B*x.

Here A and B are assumed to be Hermitian and banded, and B is also positive definite. Eigenvalues and eigenvectors can be selected by specifying either all eigenvalues, a range of values or a range of indices for the desired eigenvalues.

Parameters

in
jobz

= ‘N’: Compute eigenvalues only = ‘V’: Compute eigenvalues and eigenvectors

in
range

= ‘A’: all eigenvalues will be found = ‘V’: all eigenvalues in (vl,vu] will be found = ‘I’: the il-th through iu-th eigenvalues will be found

in
uplo

= ‘U’: Upper triangles of A and B are stored = ‘L’: Lower triangles of A and B are stored

in
n

The order of the matrices A and B. n >= 0.

in
ka

The number of super-/sub-diagonals of A. ka >= 0.

in
kb

The number of super-/sub-diagonals of B. kb >= 0.

inout
AB

Complex array, dimension (ldab, n). On entry, the Hermitian band matrix A. On exit, contents are destroyed.

in
ldab

The leading dimension of AB. ldab >= ka+1.

inout
BB

Complex array, dimension (ldbb, n). On entry, the Hermitian band matrix B. On exit, the split Cholesky factor S from cpbstf.

in
ldbb

The leading dimension of BB. ldbb >= kb+1.

out
Q

Complex array, dimension (ldq, n). If jobz=’V’, the transformation matrix.

in
ldq

The leading dimension of Q. ldq >= 1, and if jobz=’V’, ldq >= max(1,n).

in
vl

If range=’V’, the lower bound of the interval.

in
vu

If range=’V’, the upper bound of the interval. vl < vu.

in
il

If range=’I’, the index of the smallest eigenvalue.

in
iu

If range=’I’, the index of the largest eigenvalue.

in
abstol

The absolute error tolerance for eigenvalues.

out
m

The total number of eigenvalues found.

out
W

The eigenvalues in ascending order.

out
Z

Complex array, dimension (ldz, n). If jobz=’V’, the eigenvectors.

in
ldz

The leading dimension of Z. ldz >= 1, and if jobz=’V’, ldz >= max(1,n).

out
work

Complex workspace array of dimension (n).

out
rwork

Single precision workspace array of dimension (7*n).

out
iwork

Integer workspace array of dimension (5*n).

out
ifail

If jobz=’V’, indices of eigenvectors that failed to converge. Array of dimension (n).

out
info

  • = 0: successful exit

  • < 0: if info = -i, the i-th argument had an illegal value

  • <= n: i eigenvectors failed to converge

  • > n: cpbstf returned info = i (B not positive definite)

Functions

void zhbgvx(
    const char*          jobz,
    const char*          range,
    const char*          uplo,
    const INT            n,
    const INT            ka,
    const INT            kb,
          c128* restrict AB,
    const INT            ldab,
          c128* restrict BB,
    const INT            ldbb,
          c128* restrict Q,
    const INT            ldq,
    const f64            vl,
    const f64            vu,
    const INT            il,
    const INT            iu,
    const f64            abstol,
          INT*           m,
          f64*  restrict W,
          c128* restrict Z,
    const INT            ldz,
          c128* restrict work,
          f64*  restrict rwork,
          INT*  restrict iwork,
          INT*  restrict ifail,
          INT*           info
);
void zhbgvx(const char *jobz, const char *range, const char *uplo, const INT n, const INT ka, const INT kb, c128 *restrict AB, const INT ldab, c128 *restrict BB, const INT ldbb, c128 *restrict Q, const INT ldq, const f64 vl, const f64 vu, const INT il, const INT iu, const f64 abstol, INT *m, f64 *restrict W, c128 *restrict Z, const INT ldz, c128 *restrict work, f64 *restrict rwork, INT *restrict iwork, INT *restrict ifail, INT *info)#

ZHBGVX computes all the eigenvalues, and optionally, the eigenvectors of a complex generalized Hermitian-definite banded eigenproblem, of the form A*x=(lambda)*B*x.

Here A and B are assumed to be Hermitian and banded, and B is also positive definite. Eigenvalues and eigenvectors can be selected by specifying either all eigenvalues, a range of values or a range of indices for the desired eigenvalues.

Parameters

in
jobz

= ‘N’: Compute eigenvalues only = ‘V’: Compute eigenvalues and eigenvectors

in
range

= ‘A’: all eigenvalues will be found = ‘V’: all eigenvalues in (vl,vu] will be found = ‘I’: the il-th through iu-th eigenvalues will be found

in
uplo

= ‘U’: Upper triangles of A and B are stored = ‘L’: Lower triangles of A and B are stored

in
n

The order of the matrices A and B. n >= 0.

in
ka

The number of super-/sub-diagonals of A. ka >= 0.

in
kb

The number of super-/sub-diagonals of B. kb >= 0.

inout
AB

Complex array, dimension (ldab, n). On entry, the Hermitian band matrix A. On exit, contents are destroyed.

in
ldab

The leading dimension of AB. ldab >= ka+1.

inout
BB

Complex array, dimension (ldbb, n). On entry, the Hermitian band matrix B. On exit, the split Cholesky factor S from zpbstf.

in
ldbb

The leading dimension of BB. ldbb >= kb+1.

out
Q

Complex array, dimension (ldq, n). If jobz=’V’, the transformation matrix.

in
ldq

The leading dimension of Q. ldq >= 1, and if jobz=’V’, ldq >= max(1,n).

in
vl

If range=’V’, the lower bound of the interval.

in
vu

If range=’V’, the upper bound of the interval. vl < vu.

in
il

If range=’I’, the index of the smallest eigenvalue.

in
iu

If range=’I’, the index of the largest eigenvalue.

in
abstol

The absolute error tolerance for eigenvalues.

out
m

The total number of eigenvalues found.

out
W

The eigenvalues in ascending order.

out
Z

Complex array, dimension (ldz, n). If jobz=’V’, the eigenvectors.

in
ldz

The leading dimension of Z. ldz >= 1, and if jobz=’V’, ldz >= max(1,n).

out
work

Complex workspace array of dimension (n).

out
rwork

Double precision workspace array of dimension (7*n).

out
iwork

Integer workspace array of dimension (5*n).

out
ifail

If jobz=’V’, indices of eigenvectors that failed to converge. Array of dimension (n).

out
info

  • = 0: successful exit

  • < 0: if info = -i, the i-th argument had an illegal value

  • <= n: i eigenvectors failed to converge

  • > n: zpbstf returned info = i (B not positive definite)

On this page

This Page