NAG Library Routine Document

F12FEF

Note: this routine uses optional parameters to define choices in the problem specification. If you wish to use default settings for all of the optional parameters, then the option setting routine F12FDF need not be called. If, however, you wish to reset some or all of the settings please refer to Section 10 in F12FDF for a detailed description of the specification of the optional parameters.

+− Contents

1 Purpose

2 Specification

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Further Comments

+− 9 Example

9.1 Program Text

9.2 Program Data

9.3 Program Results

1 Purpose

F12FEF can be used to return additional monitoring information during computation. It is in a suite of routines which includes F12FAF, F12FBF, F12FCF and F12FDF.

2 Specification

SUBROUTINE F12FEF (

NITER, NCONV, RITZ, RZEST, ICOMM, COMM)

INTEGER	NITER, NCONV, ICOMM(*)
REAL (KIND=nag_wp)	RITZ(), RZEST(), COMM(*)

3 Description

The suite of routines is designed to calculate some of the eigenvalues,

λ

, (and optionally the corresponding eigenvectors,

x

) of a standard eigenvalue problem

A x = λ x

, or of a generalized eigenvalue problem

A x = λ B x

of order

n

, where

n

is large and the coefficient matrices

A

and

B

are sparse, real and symmetric. The suite can also be used to find selected eigenvalues/eigenvectors of smaller scale dense, real and symmetric problems.

On an intermediate exit from F12FBF with

IREVCM = 4

, F12FEF may be called to return monitoring information on the progress of the Arnoldi iterative process. The information returned by F12FEF is:

–	the number of the current Arnoldi iteration;
–	the number of converged eigenvalues at this point;
–	the real and imaginary parts of the converged eigenvalues;
–	the error bounds on the converged eigenvalues.

F12FEF does not have an equivalent routine from the ARPACK package which prints various levels of detail of monitoring information through an output channel controlled via a parameter value (see Lehoucq et al. (1998) for details of ARPACK routines). F12FEF should not be called at any time other than immediately following an

IREVCM = 4

return from F12FBF.

4 References

Lehoucq R B (2001) Implicitly restarted Arnoldi methods and subspace iteration SIAM Journal on Matrix Analysis and Applications 23 551–562

Lehoucq R B and Scott J A (1996) An evaluation of software for computing eigenvalues of sparse nonsymmetric matrices Preprint MCS-P547-1195 Argonne National Laboratory

Lehoucq R B and Sorensen D C (1996) Deflation techniques for an implicitly restarted Arnoldi iteration SIAM Journal on Matrix Analysis and Applications 17 789–821

Lehoucq R B, Sorensen D C and Yang C (1998) ARPACK Users' Guide: Solution of Large-scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods SIAM, Philidelphia

5 Parameters

1: NITER – INTEGEROutput: On exit: the number of the current Arnoldi iteration.
2: NCONV – INTEGEROutput: On exit: the number of converged eigenvalues so far.
3: RITZ( $*$ ) – REAL (KIND=nag_wp) arrayOutput: Note: the dimension of the array RITZ must be at least $NCV$ (see F12FAF).
On exit: the first NCONV locations of the array RITZ contain the real converged approximate eigenvalues.
4: RZEST( $*$ ) – REAL (KIND=nag_wp) arrayOutput: Note: the dimension of the array RZEST must be at least $NCV$ (see F12FAF).
On exit: the first NCONV locations of the array RZEST contain the Ritz estimates (error bounds) on the real NCONV converged approximate eigenvalues.
5: ICOMM( $*$ ) – INTEGER arrayCommunication Array: Note: the dimension of the array ICOMM must be at least $\max (1, LICOMM)$ , where LICOMM is passed to the setup routine (see F12FAF).
On entry: the array ICOMM output by the preceding call to F12FBF.
6: COMM( $*$ ) – REAL (KIND=nag_wp) arrayCommunication Array: Note: the dimension of the array COMM must be at least $\max (1, LCOMM)$ , where LCOMM is passed to the setup routine (see F12FAF).
On entry: the array COMM output by the preceding call to F12FBF.

6 Error Indicators and Warnings

None.

7 Accuracy

A Ritz value,

λ

, is deemed to have converged if its Ritz estimate

\leq Tolerance \times |λ|

. The default Tolerance used is the machine precision given by X02AJF.

8 Further Comments

None.

9 Example

This example solves

K x = λ K_{G} x

using the Buckling option (see F12FDF, where

K

and

K_{G}

are obtained by the finite element method applied to the one-dimensional discrete Laplacian operator

\frac{\partial^{2} u}{\partial x^{2}}

[0, 1]

, with zero Dirichlet boundary conditions using piecewise linear elements. The shift,

σ

, is a real number, and the operator used in the Buckling iterative process is

OP = inv (K - σ K_{G}) \times K

and

B = K

NAG Library Routine DocumentF12FEF