FNLMath : Eigensystem Analysis : The Base Class ARPACKBASE
The Base Class ARPACKBASE
The components of the derived type ARPACKBASE contain data used by the ARPACK routines. These will have initial or default values assigned. The default values can usually be left unchanged with a first use of our codes. They are used as arguments to the original routines of the ARPACK package. The more experienced user may wish to change the components marked with '=>' to new values, depending on their application. These can be changed prior to calling the ARPACK interface codes we provide. This base class can be extended to pass user data or procedure pointers for use within the array function.
Note that the derived type argument EXTYPE, is optional in all the ARPACK_ eigenvalue routines, but it is not optional for the user-written array functions. If EXTYPE is not included in the argument list of the ARPACK_ eigenvalue routine, an internally declared type is passed to the array functions as the argument, EXTYPE. Although the user may choose not to use this optional argument when calling our interface routines, they must include this argument in their user-supplied array function code. In this case, the array function code does not need to reference this argument.
 
TYPE, PUBLIC :: ARPACKBASE
INTEGER :: TASK = 0 ! Local store in Class for compute
! tasks to follow. Used in ARPACK_SVD.
INTEGER :: MROWS=0 ! Defines output vector size
INTEGER :: NCOLS=0 ! Defines input/output vector size
=> COMPLEX(DKIND) :: SHIFT=&
(0._DKIND, 0._DKIND) ! Shift factor
REAL(DKIND) :: TOL=EPSILON(0._DKIND) ! Error tolerance
INTEGER :: ISHFTS = 1 ! Number of shifts - don't change
=> INTEGER :: MAXITR = HUGE(1) ! Max number of iterations (many!)
=> INTEGER :: MAXMV = HUGE(1) ! Max number of matrix ops (many!)
=> INTEGER :: INFO = 0 ! ARPACK error flag, = 0 is OK
=> INTEGER :: NACC = 0 ! Number of accurate eigenvalues or
! singular values computed.
INTEGER :: IPARAM(11)=0 ! ARPACK array of direction and
! result flags
=> REAL(DKIND) :: FACTOR_MAXNCV = 2.5_DKIND ! Factor for the Number of Ritz
! vectors > the number of requested
! eigenvalues
LOGICAL :: RALEIGH_QUOTIENT = .TRUE. ! Compute eigenvalues using the
! computed eigenvectors and Raleigh
! quotients.
REAL(DKIND), ALLOCATABLE :: RESID(:) ! REAL Starting vector for Arnoldi
! iteration, if allocated; else
! random
COMPLEX(DKIND), ALLOCATABLE :: ZRESID(:) ! COMPLEX Starting vector for Arnoldi
! iteration, if allocated; else
! random
END TYPE