p?laevswp

Moves the eigenvectors from where they are computed to ScaLAPACK standard block cyclic array.

Syntax

call pslaevswp(n, zin, ldzi, z, iz, jz, descz, nvs, key, rwork, lrwork)

call pdlaevswp(n, zin, ldzi, z, iz, jz, descz, nvs, key, rwork, lrwork)

call pclaevswp(n, zin, ldzi, z, iz, jz, descz, nvs, key, rwork, lrwork)

call pzlaevswp(n, zin, ldzi, z, iz, jz, descz, nvs, key, rwork, lrwork)

Include Files

C: mkl_scalapack.h

Description

The p?laevswp routine moves the eigenvectors (potentially unsorted) from where they are computed, to a ScaLAPACK standard block cyclic array, sorted so that the corresponding eigenvalues are sorted.

Input Parameters

np = the number of rows local to a given process.

nq = the number of columns local to a given process.

n

(global). INTEGER.

The order of the matrix A. n ≥ 0.

zin

(local).

REAL for pslaevswp

DOUBLE PRECISION for pdlaevswp

COMPLEX for pclaevswp

COMPLEX*16 for pzlaevswp. Array, DIMENSION (ldzi, nvs(iam) ). The eigenvectors on input. Each eigenvector resides entirely in one process. Each process holds a contiguous set of nvs(iam) eigenvectors. The first eigenvector which the process holds is: sum for i=[0, iam-1) of nvs(i).

ldzi

(local)

INTEGER.The leading dimension of the zin array.

iz, jz

(global) INTEGER.The row and column indices in the global array Z indicating the first row and the first column of the submatrix Z, respectively.

descz

(global and local) INTEGER array, DIMENSION (dlen_). The array descriptor for the distributed matrix Z.

nvs

(global) INTEGER.

Array, DIMENSION( nprocs+1)

nvs(i) = number of processes number of eigenvectors held by processes [0, i-1)

nvs(1) = number of eigen vectors held by[0, 1 -1) = 0

nvs(nprocs+1)= number of eigen vectors held by [0, nprocs)= total number of eigenvectors.

key

(global) INTEGER.

Array, DIMENSION (n). Indicates the actual index (after sorting) for each of the eigenvectors.

rwork

(local).

REAL for pslaevswp

DOUBLE PRECISIONfor pdlaevswp

COMPLEX for pclaevswp

COMPLEX*16 for pzlaevswp. Array, DIMENSION (lrwork).

lrwork

(local)

INTEGER. Dimension of work.

Output Parameters

z

(local).

REAL for pslaevswp

DOUBLE PRECISION for pdlaevswp

COMPLEX for pclaevswp

COMPLEX*16 for pzlaevswp.

Array, global DIMENSION (n, n), local DIMENSION (descz(dlen_), nq). The eigenvectors on output. The eigenvectors are distributed in a block cyclic manner in both dimensions, with a block size of nb.