ScaLAPACKFX
|
Solves generalized Hermitian eigenvalue problem by the divide and conquer algorithm. More...
Private Member Functions | |
subroutine | scalafx_phegvd_complex (aa, desca, bb, descb, ww, zz, descz, jobz, uplo, ia, ja, ib, jb, iz, jz, work, rwork, iwork, allocfix, skipchol, info) |
Solves Hermitian eigenvalue problem by the divide and conquer algorithm. | |
subroutine | scalafx_phegvd_dcomplex (aa, desca, bb, descb, ww, zz, descz, jobz, uplo, ia, ja, ib, jb, iz, jz, work, rwork, iwork, allocfix, skipchol, info) |
Solves Hermitian eigenvalue problem by the divide and conquer algorithm. | |
Solves generalized Hermitian eigenvalue problem by the divide and conquer algorithm.
|
private |
Solves Hermitian eigenvalue problem by the divide and conquer algorithm.
Invokes SCALAPACK routines p?potrf, p?hegst, p?heevd, p?trsm in order to transform the general eigenvalue problem to the standard form and transform the eigenvectors back.
aa | Matrix to diagonalize (A), transformed matrix on exit. |
desca | Descriptor of matrix A. |
bb | Matrix on the right hand side (B), transformed matrix on exit. |
descb | Descriptor of matrix B. |
ww | Eigenvalues on exit. |
zz | Eigenvectors on exit (Z). |
descz | Descriptor of the eigenvector matrix. |
jobz | Job type (default: "V") |
uplo | Upper or lower diagonal matrix (default: "L") |
ia | First row of the submatrix A (default: 1) |
ja | First column of the submatrix A (default: 1) |
ib | First row of the submatrix B (default: 1) |
jb | First column of the submatrix B (default: 1) |
iz | First row of the submatrix Z (default: 1) |
jz | First column of the submatrix Z (default: 1) |
work | Complex working array (if not specified, allocated automatically) |
rwork | Real working array (if not specified, allocated automatically) |
iwork | Integer working array (if not specified, allocated automatically) |
allocfix | If yes, the routine tries to enlarge the workspace size as returned by the appropriate p?syevd() routine by some empirical values. See the scalafx_pheevd_complex() routine for details. |
skipchol | If true, the Cholesky transformation will be skipped. Array bb must have the Cholesky transformed form. |
info | Info flag. If not specified and SCALAPACK calls returns nozero, subroutine stops. |
|
private |
Solves Hermitian eigenvalue problem by the divide and conquer algorithm.
Invokes SCALAPACK routines p?potrf, p?hegst, p?heevd, p?trsm in order to transform the general eigenvalue problem to the standard form and transform the eigenvectors back.
aa | Matrix to diagonalize (A), transformed matrix on exit. |
desca | Descriptor of matrix A. |
bb | Matrix on the right hand side (B), transformed matrix on exit. |
descb | Descriptor of matrix B. |
ww | Eigenvalues on exit. |
zz | Eigenvectors on exit (Z). |
descz | Descriptor of the eigenvector matrix. |
jobz | Job type (default: "V") |
uplo | Upper or lower diagonal matrix (default: "L") |
ia | First row of the submatrix A (default: 1) |
ja | First column of the submatrix A (default: 1) |
ib | First row of the submatrix B (default: 1) |
jb | First column of the submatrix B (default: 1) |
iz | First row of the submatrix Z (default: 1) |
jz | First column of the submatrix Z (default: 1) |
work | Complex working array (if not specified, allocated automatically) |
rwork | Real working array (if not specified, allocated automatically) |
iwork | Integer working array (if not specified, allocated automatically) |
allocfix | If yes, the routine tries to enlarge the workspace size as returned by the appropriate p?syevd() routine by some empirical values. See the scalafx_pheevd_dcomplex() routine for details. |
skipchol | If true, the Cholesky transformation will be skipped. Array bb must have the Cholesky transformed form. |
info | Info flag. If not specified and SCALAPACK calls returns nozero, subroutine stops. |