ScaLAPACKFX
|
Solves generalized symmetric eigenvalue problem by the divide and conquer algorithm. More...
Private Member Functions | |
subroutine | scalafx_psygvd_real (aa, desca, bb, descb, ww, zz, descz, jobz, uplo, ia, ja, ib, jb, iz, jz, work, iwork, allocfix, skipchol, info) |
Solves real generalized eigenvalue problem by the divide and conquer algorithm. | |
subroutine | scalafx_psygvd_dreal (aa, desca, bb, descb, ww, zz, descz, jobz, uplo, ia, ja, ib, jb, iz, jz, work, iwork, allocfix, skipchol, info) |
Solves real generalized eigenvalue problem by the divide and conquer algorithm. | |
Solves generalized symmetric eigenvalue problem by the divide and conquer algorithm.
|
private |
Solves real generalized eigenvalue problem by the divide and conquer algorithm.
Invokes SCALAPACK routines p?potrf, p?sygst, p?syevd, 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 | 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_psyevd_dreal() 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 real generalized eigenvalue problem by the divide and conquer algorithm.
Invokes SCALAPACK routines p?potrf, p?sygst, p?syevd, 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 | 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_psyevd_real() 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. |