NAG Library Function Document

nag_complex_sparse_eigensystem_sol (f12aqc)

1

Purpose

2

Specification

3

Description

4

References

5

Arguments

1:    $\mathbf{nconv}$ – Integer *Output

On exit: the number of converged eigenvalues as found by nag_complex_sparse_eigensystem_option (f12arc).

2:    $\mathbf{d}\left[\mathit{dim}\right]$ – ComplexOutput

Note: the dimension, dim, of the array d must be at least $\mathbf{ncv}$ (see nag_complex_sparse_eigensystem_init (f12anc)).

On exit: the first nconv locations of the array d contain the converged approximate eigenvalues.

3:    $\mathbf{z}\left[\mathbf{n}\times \mathbf{ncv}\right]$ – ComplexOutput

On exit: if the default option $\mathbf{Vectors}=\mathrm{RITZ}$ (see nag_real_sparse_eigensystem_option (f12adc)) has been selected then z contains the final set of eigenvectors corresponding to the eigenvalues held in d. The complex eigenvector associated with an eigenvalue is stored in the corresponding array section of z.

4:    $\mathbf{sigma}$ – ComplexInput

On entry: if one of the $\mathbf{Shifted\; Inverse}$ (see nag_complex_sparse_eigensystem_option (f12arc)) modes has been selected then sigma contains the shift used; otherwise sigma is not referenced.

5:    $\mathbf{resid}\left[\mathit{dim}\right]$ – const ComplexInput

Note: the dimension, dim, of the array resid must be at least $\mathbf{n}$ (see nag_complex_sparse_eigensystem_init (f12anc)).

On entry: must not be modified following a call to nag_complex_sparse_eigensystem_iter (f12apc) since it contains data required by nag_complex_sparse_eigensystem_sol (f12aqc).

6:    $\mathbf{v}\left[\mathbf{n}\times \mathbf{ncv}\right]$ – ComplexInput/Output

The $\mathit{i}$th element of the $\mathit{j}$th basis vector is stored in location $\mathbf{v}\left[\mathbf{n}\times \left(\mathit{j}-1\right)+\mathit{i}-1\right]$, for $\mathit{i}=1,2,\dots ,\mathbf{n}$ and $\mathit{j}=1,2,\dots ,\mathbf{ncv}$.

On entry: the ncv sections of v, of length $n$, contain the Arnoldi basis vectors for $\mathrm{OP}$ as constructed by nag_complex_sparse_eigensystem_iter (f12apc).

On exit: if the option $\mathbf{Vectors}=\mathrm{SCHUR}$ or $\mathrm{RITZ}$ has been set and a separate array z has been passed (i.e., z does not equal v), then the first nconv sections of v, of length $n$, will contain approximate Schur vectors that span the desired invariant subspace.

7:    $\mathbf{comm}\left[\mathit{dim}\right]$ – ComplexCommunication Array

Note: the dimension, dim, of the array comm must be at least $\max \left(1,\mathbf{lcomm}\right)$ (see nag_complex_sparse_eigensystem_init (f12anc)).

On initial entry: must remain unchanged from the prior call to nag_complex_sparse_eigensystem_init (f12anc).

On exit: contains data on the current state of the solution.

8:    $\mathbf{icomm}\left[\mathit{dim}\right]$ – IntegerCommunication Array

Note: the dimension, dim, of the array icomm must be at least $\max \left(1,\mathbf{licomm}\right)$ (see nag_complex_sparse_eigensystem_init (f12anc)).

On initial entry: must remain unchanged from the prior call to nag_complex_sparse_eigensystem_init (f12anc).

On exit: contains data on the current state of the solution.

9:    $\mathbf{fail}$ – NagError *Input/Output

The NAG error argument (see Section 3.7 in How to Use the NAG Library and its Documentation).

6

Error Indicators and Warnings

NE_ALLOC_FAIL

Dynamic memory allocation failed.
See Section 2.3.1.2 in How to Use the NAG Library and its Documentation for further information.

NE_BAD_PARAM

On entry, argument $〈\mathit{\text{value}}〉$ had an illegal value.

NE_INTERNAL_EIGVEC_FAIL

In calculating eigenvectors, an internal call returned with an error. The function returned with $\mathbf{fail}\mathbf{.}\mathbf{code}=〈\mathit{\text{value}}〉$. Please contact NAG.

NE_INTERNAL_ERROR

An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.
See Section 2.7.6 in How to Use the NAG Library and its Documentation for further information.

NE_INVALID_OPTION

On entry, $\mathbf{Vectors}=\text{Select}$, but this is not yet implemented.

NE_NO_LICENCE

Your licence key may have expired or may not have been installed correctly.
See Section 2.7.5 in How to Use the NAG Library and its Documentation for further information.

NE_RITZ_COUNT

Got a different count of the number of converged Ritz values than the value passed to it through the argument icomm: number counted $=〈\mathit{\text{value}}〉$, number expected $=〈\mathit{\text{value}}〉$. This usually indicates that a communication array has been altered or has become corrupted between calls to nag_complex_sparse_eigensystem_iter (f12apc) and nag_complex_sparse_eigensystem_sol (f12aqc).

NE_SCHUR_EIG_FAIL

During calculation of a Schur form, there was a failure to compute $〈\mathit{\text{value}}〉$ eigenvalues in a total of $〈\mathit{\text{value}}〉$ iterations.

NE_SCHUR_REORDER

The computed Schur form could not be reordered by an internal call. This function returned with $\mathbf{fail}\mathbf{.}\mathbf{code}=〈\mathit{\text{value}}〉$. Please contact NAG.

NE_ZERO_EIGS_FOUND

The number of eigenvalues found to sufficient accuracy, as communicated through the argument icomm, is zero. You should experiment with different values of nev and ncv, or select a different computational mode or increase the maximum number of iterations prior to calling nag_complex_sparse_eigensystem_iter (f12apc).

7

Accuracy

8

Parallelism and Performance

9

Further Comments

10

Example

10.1

Program Text

Program Text (f12aqce.c)

10.2

Program Data

Program Data (f12aqce.d)

10.3

Program Results

Program Results (f12aqce.r)