octave-nkf: libcruft/arpack/src/ssaupd.f annotate

annotate libcruft/arpack/src/ssaupd.f @ 12274:9f5d2ef078e8 release-3-4-x

import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1

author	John W. Eaton <jwe@octave.org>
date	Fri, 28 Jan 2011 14:04:33 -0500
parents
children

rev	line source
12274 9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1 c-----------------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	2 c\BeginDoc
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	3 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	4 c\Name: ssaupd
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	5 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	6 c\Description:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	7 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	8 c Reverse communication interface for the Implicitly Restarted Arnoldi
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	9 c Iteration. For symmetric problems this reduces to a variant of the Lanczos
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	10 c method. This method has been designed to compute approximations to a
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	11 c few eigenpairs of a linear operator OP that is real and symmetric
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	12 c with respect to a real positive semi-definite symmetric matrix B,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	13 c i.e.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	14 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	15 c BOP = (OP`)B.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	16 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	17 c Another way to express this condition is
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	18 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	19 c < x,OPy > = < OPx,y > where < z,w > = z`Bw .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	20 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	21 c In the standard eigenproblem B is the identity matrix.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	22 c ( A` denotes transpose of A)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	23 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	24 c The computed approximate eigenvalues are called Ritz values and
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	25 c the corresponding approximate eigenvectors are called Ritz vectors.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	26 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	27 c ssaupd is usually called iteratively to solve one of the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	28 c following problems:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	29 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	30 c Mode 1: Ax = lambdax, A symmetric
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	31 c ===> OP = A and B = I.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	32 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	33 c Mode 2: Ax = lambdaM*x, A symmetric, M symmetric positive definite
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	34 c ===> OP = inv[M]*A and B = M.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	35 c ===> (If M can be factored see remark 3 below)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	36 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	37 c Mode 3: Kx = lambdaM*x, K symmetric, M symmetric positive semi-definite
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	38 c ===> OP = (inv[K - sigmaM])M and B = M.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	39 c ===> Shift-and-Invert mode
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	40 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	41 c Mode 4: Kx = lambdaKG*x, K symmetric positive semi-definite,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	42 c KG symmetric indefinite
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	43 c ===> OP = (inv[K - sigmaKG])K and B = K.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	44 c ===> Buckling mode
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	45 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	46 c Mode 5: Ax = lambdaM*x, A symmetric, M symmetric positive semi-definite
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	47 c ===> OP = inv[A - sigmaM][A + sigma*M] and B = M.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	48 c ===> Cayley transformed mode
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	49 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	50 c NOTE: The action of w <- inv[A - sigmaM]v or w <- inv[M]*v
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	51 c should be accomplished either by a direct method
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	52 c using a sparse matrix factorization and solving
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	53 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	54 c [A - sigmaM]w = v or M*w = v,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	55 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	56 c or through an iterative method for solving these
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	57 c systems. If an iterative method is used, the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	58 c convergence test must be more stringent than
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	59 c the accuracy requirements for the eigenvalue
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	60 c approximations.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	61 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	62 c\Usage:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	63 c call ssaupd
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	64 c ( IDO, BMAT, N, WHICH, NEV, TOL, RESID, NCV, V, LDV, IPARAM,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	65 c IPNTR, WORKD, WORKL, LWORKL, INFO )
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	66 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	67 c\Arguments
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	68 c IDO Integer. (INPUT/OUTPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	69 c Reverse communication flag. IDO must be zero on the first
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	70 c call to ssaupd. IDO will be set internally to
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	71 c indicate the type of operation to be performed. Control is
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	72 c then given back to the calling routine which has the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	73 c responsibility to carry out the requested operation and call
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	74 c ssaupd with the result. The operand is given in
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	75 c WORKD(IPNTR(1)), the result must be put in WORKD(IPNTR(2)).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	76 c (If Mode = 2 see remark 5 below)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	77 c -------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	78 c IDO = 0: first call to the reverse communication interface
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	79 c IDO = -1: compute Y = OP * X where
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	80 c IPNTR(1) is the pointer into WORKD for X,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	81 c IPNTR(2) is the pointer into WORKD for Y.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	82 c This is for the initialization phase to force the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	83 c starting vector into the range of OP.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	84 c IDO = 1: compute Y = OP * X where
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	85 c IPNTR(1) is the pointer into WORKD for X,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	86 c IPNTR(2) is the pointer into WORKD for Y.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	87 c In mode 3,4 and 5, the vector B * X is already
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	88 c available in WORKD(ipntr(3)). It does not
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	89 c need to be recomputed in forming OP * X.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	90 c IDO = 2: compute Y = B * X where
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	91 c IPNTR(1) is the pointer into WORKD for X,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	92 c IPNTR(2) is the pointer into WORKD for Y.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	93 c IDO = 3: compute the IPARAM(8) shifts where
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	94 c IPNTR(11) is the pointer into WORKL for
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	95 c placing the shifts. See remark 6 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	96 c IDO = 99: done
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	97 c -------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	98 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	99 c BMAT Character*1. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	100 c BMAT specifies the type of the matrix B that defines the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	101 c semi-inner product for the operator OP.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	102 c B = 'I' -> standard eigenvalue problem Ax = lambdax
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	103 c B = 'G' -> generalized eigenvalue problem Ax = lambdaB*x
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	104 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	105 c N Integer. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	106 c Dimension of the eigenproblem.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	107 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	108 c WHICH Character*2. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	109 c Specify which of the Ritz values of OP to compute.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	110 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	111 c 'LA' - compute the NEV largest (algebraic) eigenvalues.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	112 c 'SA' - compute the NEV smallest (algebraic) eigenvalues.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	113 c 'LM' - compute the NEV largest (in magnitude) eigenvalues.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	114 c 'SM' - compute the NEV smallest (in magnitude) eigenvalues.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	115 c 'BE' - compute NEV eigenvalues, half from each end of the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	116 c spectrum. When NEV is odd, compute one more from the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	117 c high end than from the low end.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	118 c (see remark 1 below)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	119 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	120 c NEV Integer. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	121 c Number of eigenvalues of OP to be computed. 0 < NEV < N.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	122 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	123 c TOL Real scalar. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	124 c Stopping criterion: the relative accuracy of the Ritz value
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	125 c is considered acceptable if BOUNDS(I) .LE. TOL*ABS(RITZ(I)).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	126 c If TOL .LE. 0. is passed a default is set:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	127 c DEFAULT = SLAMCH('EPS') (machine precision as computed
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	128 c by the LAPACK auxiliary subroutine SLAMCH).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	129 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	130 c RESID Real array of length N. (INPUT/OUTPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	131 c On INPUT:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	132 c If INFO .EQ. 0, a random initial residual vector is used.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	133 c If INFO .NE. 0, RESID contains the initial residual vector,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	134 c possibly from a previous run.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	135 c On OUTPUT:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	136 c RESID contains the final residual vector.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	137 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	138 c NCV Integer. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	139 c Number of columns of the matrix V (less than or equal to N).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	140 c This will indicate how many Lanczos vectors are generated
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	141 c at each iteration. After the startup phase in which NEV
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	142 c Lanczos vectors are generated, the algorithm generates
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	143 c NCV-NEV Lanczos vectors at each subsequent update iteration.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	144 c Most of the cost in generating each Lanczos vector is in the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	145 c matrix-vector product OP*x. (See remark 4 below).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	146 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	147 c V Real N by NCV array. (OUTPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	148 c The NCV columns of V contain the Lanczos basis vectors.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	149 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	150 c LDV Integer. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	151 c Leading dimension of V exactly as declared in the calling
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	152 c program.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	153 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	154 c IPARAM Integer array of length 11. (INPUT/OUTPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	155 c IPARAM(1) = ISHIFT: method for selecting the implicit shifts.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	156 c The shifts selected at each iteration are used to restart
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	157 c the Arnoldi iteration in an implicit fashion.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	158 c -------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	159 c ISHIFT = 0: the shifts are provided by the user via
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	160 c reverse communication. The NCV eigenvalues of
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	161 c the current tridiagonal matrix T are returned in
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	162 c the part of WORKL array corresponding to RITZ.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	163 c See remark 6 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	164 c ISHIFT = 1: exact shifts with respect to the reduced
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	165 c tridiagonal matrix T. This is equivalent to
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	166 c restarting the iteration with a starting vector
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	167 c that is a linear combination of Ritz vectors
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	168 c associated with the "wanted" Ritz values.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	169 c -------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	170 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	171 c IPARAM(2) = LEVEC
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	172 c No longer referenced. See remark 2 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	173 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	174 c IPARAM(3) = MXITER
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	175 c On INPUT: maximum number of Arnoldi update iterations allowed.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	176 c On OUTPUT: actual number of Arnoldi update iterations taken.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	177 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	178 c IPARAM(4) = NB: blocksize to be used in the recurrence.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	179 c The code currently works only for NB = 1.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	180 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	181 c IPARAM(5) = NCONV: number of "converged" Ritz values.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	182 c This represents the number of Ritz values that satisfy
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	183 c the convergence criterion.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	184 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	185 c IPARAM(6) = IUPD
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	186 c No longer referenced. Implicit restarting is ALWAYS used.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	187 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	188 c IPARAM(7) = MODE
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	189 c On INPUT determines what type of eigenproblem is being solved.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	190 c Must be 1,2,3,4,5; See under \Description of ssaupd for the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	191 c five modes available.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	192 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	193 c IPARAM(8) = NP
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	194 c When ido = 3 and the user provides shifts through reverse
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	195 c communication (IPARAM(1)=0), ssaupd returns NP, the number
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	196 c of shifts the user is to provide. 0 < NP <=NCV-NEV. See Remark
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	197 c 6 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	198 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	199 c IPARAM(9) = NUMOP, IPARAM(10) = NUMOPB, IPARAM(11) = NUMREO,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	200 c OUTPUT: NUMOP = total number of OP*x operations,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	201 c NUMOPB = total number of B*x operations if BMAT='G',
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	202 c NUMREO = total number of steps of re-orthogonalization.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	203 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	204 c IPNTR Integer array of length 11. (OUTPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	205 c Pointer to mark the starting locations in the WORKD and WORKL
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	206 c arrays for matrices/vectors used by the Lanczos iteration.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	207 c -------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	208 c IPNTR(1): pointer to the current operand vector X in WORKD.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	209 c IPNTR(2): pointer to the current result vector Y in WORKD.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	210 c IPNTR(3): pointer to the vector B * X in WORKD when used in
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	211 c the shift-and-invert mode.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	212 c IPNTR(4): pointer to the next available location in WORKL
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	213 c that is untouched by the program.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	214 c IPNTR(5): pointer to the NCV by 2 tridiagonal matrix T in WORKL.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	215 c IPNTR(6): pointer to the NCV RITZ values array in WORKL.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	216 c IPNTR(7): pointer to the Ritz estimates in array WORKL associated
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	217 c with the Ritz values located in RITZ in WORKL.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	218 c IPNTR(11): pointer to the NP shifts in WORKL. See Remark 6 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	219 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	220 c Note: IPNTR(8:10) is only referenced by sseupd. See Remark 2.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	221 c IPNTR(8): pointer to the NCV RITZ values of the original system.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	222 c IPNTR(9): pointer to the NCV corresponding error bounds.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	223 c IPNTR(10): pointer to the NCV by NCV matrix of eigenvectors
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	224 c of the tridiagonal matrix T. Only referenced by
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	225 c sseupd if RVEC = .TRUE. See Remarks.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	226 c -------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	227 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	228 c WORKD Real work array of length 3*N. (REVERSE COMMUNICATION)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	229 c Distributed array to be used in the basic Arnoldi iteration
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	230 c for reverse communication. The user should not use WORKD
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	231 c as temporary workspace during the iteration. Upon termination
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	232 c WORKD(1:N) contains B*RESID(1:N). If the Ritz vectors are desired
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	233 c subroutine sseupd uses this output.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	234 c See Data Distribution Note below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	235 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	236 c WORKL Real work array of length LWORKL. (OUTPUT/WORKSPACE)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	237 c Private (replicated) array on each PE or array allocated on
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	238 c the front end. See Data Distribution Note below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	239 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	240 c LWORKL Integer. (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	241 c LWORKL must be at least NCV*2 + 8NCV .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	242 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	243 c INFO Integer. (INPUT/OUTPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	244 c If INFO .EQ. 0, a randomly initial residual vector is used.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	245 c If INFO .NE. 0, RESID contains the initial residual vector,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	246 c possibly from a previous run.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	247 c Error flag on output.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	248 c = 0: Normal exit.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	249 c = 1: Maximum number of iterations taken.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	250 c All possible eigenvalues of OP has been found. IPARAM(5)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	251 c returns the number of wanted converged Ritz values.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	252 c = 2: No longer an informational error. Deprecated starting
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	253 c with release 2 of ARPACK.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	254 c = 3: No shifts could be applied during a cycle of the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	255 c Implicitly restarted Arnoldi iteration. One possibility
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	256 c is to increase the size of NCV relative to NEV.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	257 c See remark 4 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	258 c = -1: N must be positive.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	259 c = -2: NEV must be positive.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	260 c = -3: NCV must be greater than NEV and less than or equal to N.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	261 c = -4: The maximum number of Arnoldi update iterations allowed
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	262 c must be greater than zero.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	263 c = -5: WHICH must be one of 'LM', 'SM', 'LA', 'SA' or 'BE'.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	264 c = -6: BMAT must be one of 'I' or 'G'.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	265 c = -7: Length of private work array WORKL is not sufficient.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	266 c = -8: Error return from trid. eigenvalue calculation;
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	267 c Informatinal error from LAPACK routine ssteqr.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	268 c = -9: Starting vector is zero.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	269 c = -10: IPARAM(7) must be 1,2,3,4,5.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	270 c = -11: IPARAM(7) = 1 and BMAT = 'G' are incompatable.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	271 c = -12: IPARAM(1) must be equal to 0 or 1.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	272 c = -13: NEV and WHICH = 'BE' are incompatable.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	273 c = -9999: Could not build an Arnoldi factorization.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	274 c IPARAM(5) returns the size of the current Arnoldi
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	275 c factorization. The user is advised to check that
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	276 c enough workspace and array storage has been allocated.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	277 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	278 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	279 c\Remarks
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	280 c 1. The converged Ritz values are always returned in ascending
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	281 c algebraic order. The computed Ritz values are approximate
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	282 c eigenvalues of OP. The selection of WHICH should be made
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	283 c with this in mind when Mode = 3,4,5. After convergence,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	284 c approximate eigenvalues of the original problem may be obtained
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	285 c with the ARPACK subroutine sseupd.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	286 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	287 c 2. If the Ritz vectors corresponding to the converged Ritz values
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	288 c are needed, the user must call sseupd immediately following completion
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	289 c of ssaupd. This is new starting with version 2.1 of ARPACK.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	290 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	291 c 3. If M can be factored into a Cholesky factorization M = LL`
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	292 c then Mode = 2 should not be selected. Instead one should use
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	293 c Mode = 1 with OP = inv(L)Ainv(L`). Appropriate triangular
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	294 c linear systems should be solved with L and L` rather
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	295 c than computing inverses. After convergence, an approximate
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	296 c eigenvector z of the original problem is recovered by solving
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	297 c L`z = x where x is a Ritz vector of OP.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	298 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	299 c 4. At present there is no a-priori analysis to guide the selection
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	300 c of NCV relative to NEV. The only formal requrement is that NCV > NEV.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	301 c However, it is recommended that NCV .ge. 2*NEV. If many problems of
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	302 c the same type are to be solved, one should experiment with increasing
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	303 c NCV while keeping NEV fixed for a given test problem. This will
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	304 c usually decrease the required number of OP*x operations but it
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	305 c also increases the work and storage required to maintain the orthogonal
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	306 c basis vectors. The optimal "cross-over" with respect to CPU time
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	307 c is problem dependent and must be determined empirically.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	308 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	309 c 5. If IPARAM(7) = 2 then in the Reverse commuication interface the user
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	310 c must do the following. When IDO = 1, Y = OP * X is to be computed.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	311 c When IPARAM(7) = 2 OP = inv(B)A. After computing AX the user
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	312 c must overwrite X with A*X. Y is then the solution to the linear set
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	313 c of equations BY = AX.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	314 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	315 c 6. When IPARAM(1) = 0, and IDO = 3, the user needs to provide the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	316 c NP = IPARAM(8) shifts in locations:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	317 c 1 WORKL(IPNTR(11))
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	318 c 2 WORKL(IPNTR(11)+1)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	319 c .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	320 c .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	321 c .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	322 c NP WORKL(IPNTR(11)+NP-1).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	323 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	324 c The eigenvalues of the current tridiagonal matrix are located in
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	325 c WORKL(IPNTR(6)) through WORKL(IPNTR(6)+NCV-1). They are in the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	326 c order defined by WHICH. The associated Ritz estimates are located in
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	327 c WORKL(IPNTR(8)), WORKL(IPNTR(8)+1), ... , WORKL(IPNTR(8)+NCV-1).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	328 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	329 c-----------------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	330 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	331 c\Data Distribution Note:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	332 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	333 c Fortran-D syntax:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	334 c ================
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	335 c REAL RESID(N), V(LDV,NCV), WORKD(3*N), WORKL(LWORKL)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	336 c DECOMPOSE D1(N), D2(N,NCV)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	337 c ALIGN RESID(I) with D1(I)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	338 c ALIGN V(I,J) with D2(I,J)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	339 c ALIGN WORKD(I) with D1(I) range (1:N)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	340 c ALIGN WORKD(I) with D1(I-N) range (N+1:2*N)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	341 c ALIGN WORKD(I) with D1(I-2N) range (2N+1:3*N)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	342 c DISTRIBUTE D1(BLOCK), D2(BLOCK,:)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	343 c REPLICATED WORKL(LWORKL)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	344 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	345 c Cray MPP syntax:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	346 c ===============
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	347 c REAL RESID(N), V(LDV,NCV), WORKD(N,3), WORKL(LWORKL)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	348 c SHARED RESID(BLOCK), V(BLOCK,:), WORKD(BLOCK,:)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	349 c REPLICATED WORKL(LWORKL)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	350 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	351 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	352 c\BeginLib
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	353 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	354 c\References:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	355 c 1. D.C. Sorensen, "Implicit Application of Polynomial Filters in
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	356 c a k-Step Arnoldi Method", SIAM J. Matr. Anal. Apps., 13 (1992),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	357 c pp 357-385.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	358 c 2. R.B. Lehoucq, "Analysis and Implementation of an Implicitly
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	359 c Restarted Arnoldi Iteration", Rice University Technical Report
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	360 c TR95-13, Department of Computational and Applied Mathematics.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	361 c 3. B.N. Parlett, "The Symmetric Eigenvalue Problem". Prentice-Hall,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	362 c 1980.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	363 c 4. B.N. Parlett, B. Nour-Omid, "Towards a Black Box Lanczos Program",
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	364 c Computer Physics Communications, 53 (1989), pp 169-179.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	365 c 5. B. Nour-Omid, B.N. Parlett, T. Ericson, P.S. Jensen, "How to
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	366 c Implement the Spectral Transformation", Math. Comp., 48 (1987),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	367 c pp 663-673.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	368 c 6. R.G. Grimes, J.G. Lewis and H.D. Simon, "A Shifted Block Lanczos
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	369 c Algorithm for Solving Sparse Symmetric Generalized Eigenproblems",
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	370 c SIAM J. Matr. Anal. Apps., January (1993).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	371 c 7. L. Reichel, W.B. Gragg, "Algorithm 686: FORTRAN Subroutines
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	372 c for Updating the QR decomposition", ACM TOMS, December 1990,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	373 c Volume 16 Number 4, pp 369-377.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	374 c 8. R.B. Lehoucq, D.C. Sorensen, "Implementation of Some Spectral
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	375 c Transformations in a k-Step Arnoldi Method". In Preparation.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	376 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	377 c\Routines called:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	378 c ssaup2 ARPACK routine that implements the Implicitly Restarted
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	379 c Arnoldi Iteration.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	380 c sstats ARPACK routine that initialize timing and other statistics
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	381 c variables.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	382 c ivout ARPACK utility routine that prints integers.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	383 c arscnd ARPACK utility routine for timing.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	384 c svout ARPACK utility routine that prints vectors.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	385 c slamch LAPACK routine that determines machine constants.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	386 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	387 c\Authors
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	388 c Danny Sorensen Phuong Vu
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	389 c Richard Lehoucq CRPC / Rice University
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	390 c Dept. of Computational & Houston, Texas
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	391 c Applied Mathematics
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	392 c Rice University
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	393 c Houston, Texas
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	394 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	395 c\Revision history:
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	396 c 12/15/93: Version ' 2.4'
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	397 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	398 c\SCCS Information: @(#)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	399 c FILE: saupd.F SID: 2.8 DATE OF SID: 04/10/01 RELEASE: 2
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	400 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	401 c\Remarks
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	402 c 1. None
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	403 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	404 c\EndLib
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	405 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	406 c-----------------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	407 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	408 subroutine ssaupd
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	409 & ( ido, bmat, n, which, nev, tol, resid, ncv, v, ldv, iparam,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	410 & ipntr, workd, workl, lworkl, info )
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	411 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	412 c %----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	413 c \| Include files for debugging and timing information \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	414 c %----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	415 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	416 include 'debug.h'
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	417 include 'stat.h'
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	418 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	419 c %------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	420 c \| Scalar Arguments \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	421 c %------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	422 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	423 character bmat1, which2
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	424 integer ido, info, ldv, lworkl, n, ncv, nev
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	425 Real
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	426 & tol
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	427 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	428 c %-----------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	429 c \| Array Arguments \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	430 c %-----------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	431 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	432 integer iparam(11), ipntr(11)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	433 Real
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	434 & resid(n), v(ldv,ncv), workd(3*n), workl(lworkl)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	435 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	436 c %------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	437 c \| Parameters \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	438 c %------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	439 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	440 Real
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	441 & one, zero
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	442 parameter (one = 1.0E+0 , zero = 0.0E+0 )
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	443 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	444 c %---------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	445 c \| Local Scalars \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	446 c %---------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	447 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	448 integer bounds, ierr, ih, iq, ishift, iupd, iw,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	449 & ldh, ldq, msglvl, mxiter, mode, nb,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	450 & nev0, next, np, ritz, j
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	451 save bounds, ierr, ih, iq, ishift, iupd, iw,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	452 & ldh, ldq, msglvl, mxiter, mode, nb,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	453 & nev0, next, np, ritz
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	454 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	455 c %----------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	456 c \| External Subroutines \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	457 c %----------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	458 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	459 external ssaup2, svout, ivout, arscnd, sstats
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	460 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	461 c %--------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	462 c \| External Functions \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	463 c %--------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	464 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	465 Real
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	466 & slamch
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	467 external slamch
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	468 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	469 c %-----------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	470 c \| Executable Statements \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	471 c %-----------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	472 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	473 if (ido .eq. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	474 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	475 c %-------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	476 c \| Initialize timing statistics \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	477 c \| & message level for debugging \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	478 c %-------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	479 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	480 call sstats
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	481 call arscnd (t0)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	482 msglvl = msaupd
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	483 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	484 ierr = 0
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	485 ishift = iparam(1)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	486 mxiter = iparam(3)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	487 c nb = iparam(4)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	488 nb = 1
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	489 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	490 c %--------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	491 c \| Revision 2 performs only implicit restart. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	492 c %--------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	493 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	494 iupd = 1
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	495 mode = iparam(7)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	496 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	497 c %----------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	498 c \| Error checking \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	499 c %----------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	500 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	501 if (n .le. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	502 ierr = -1
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	503 else if (nev .le. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	504 ierr = -2
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	505 else if (ncv .le. nev .or. ncv .gt. n) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	506 ierr = -3
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	507 end if
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	508 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	509 c %----------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	510 c \| NP is the number of additional steps to \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	511 c \| extend the length NEV Lanczos factorization. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	512 c %----------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	513 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	514 np = ncv - nev
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	515 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	516 if (mxiter .le. 0) ierr = -4
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	517 if (which .ne. 'LM' .and.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	518 & which .ne. 'SM' .and.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	519 & which .ne. 'LA' .and.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	520 & which .ne. 'SA' .and.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	521 & which .ne. 'BE') ierr = -5
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	522 if (bmat .ne. 'I' .and. bmat .ne. 'G') ierr = -6
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	523 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	524 if (lworkl .lt. ncv*2 + 8ncv) ierr = -7
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	525 if (mode .lt. 1 .or. mode .gt. 5) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	526 ierr = -10
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	527 else if (mode .eq. 1 .and. bmat .eq. 'G') then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	528 ierr = -11
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	529 else if (ishift .lt. 0 .or. ishift .gt. 1) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	530 ierr = -12
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	531 else if (nev .eq. 1 .and. which .eq. 'BE') then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	532 ierr = -13
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	533 end if
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	534 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	535 c %------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	536 c \| Error Exit \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	537 c %------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	538 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	539 if (ierr .ne. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	540 info = ierr
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	541 ido = 99
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	542 go to 9000
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	543 end if
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	544 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	545 c %------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	546 c \| Set default parameters \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	547 c %------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	548 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	549 if (nb .le. 0) nb = 1
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	550 if (tol .le. zero) tol = slamch('EpsMach')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	551 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	552 c %----------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	553 c \| NP is the number of additional steps to \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	554 c \| extend the length NEV Lanczos factorization. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	555 c \| NEV0 is the local variable designating the \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	556 c \| size of the invariant subspace desired. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	557 c %----------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	558 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	559 np = ncv - nev
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	560 nev0 = nev
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	561 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	562 c %-----------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	563 c \| Zero out internal workspace \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	564 c %-----------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	565 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	566 do 10 j = 1, ncv*2 + 8ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	567 workl(j) = zero
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	568 10 continue
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	569 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	570 c %-------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	571 c \| Pointer into WORKL for address of H, RITZ, BOUNDS, Q \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	572 c \| etc... and the remaining workspace. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	573 c \| Also update pointer to be used on output. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	574 c \| Memory is laid out as follows: \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	575 c \| workl(1:2*ncv) := generated tridiagonal matrix \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	576 c \| workl(2ncv+1:2ncv+ncv) := ritz values \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	577 c \| workl(3ncv+1:3ncv+ncv) := computed error bounds \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	578 c \| workl(4ncv+1:4ncv+ncv*ncv) := rotation matrix Q \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	579 c \| workl(4ncv+ncvncv+1:7ncv+ncvncv) := workspace \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	580 c %-------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	581 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	582 ldh = ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	583 ldq = ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	584 ih = 1
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	585 ritz = ih + 2*ldh
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	586 bounds = ritz + ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	587 iq = bounds + ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	588 iw = iq + ncv**2
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	589 next = iw + 3*ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	590 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	591 ipntr(4) = next
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	592 ipntr(5) = ih
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	593 ipntr(6) = ritz
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	594 ipntr(7) = bounds
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	595 ipntr(11) = iw
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	596 end if
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	597 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	598 c %-------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	599 c \| Carry out the Implicitly restarted Lanczos Iteration. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	600 c %-------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	601 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	602 call ssaup2
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	603 & ( ido, bmat, n, which, nev0, np, tol, resid, mode, iupd,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	604 & ishift, mxiter, v, ldv, workl(ih), ldh, workl(ritz),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	605 & workl(bounds), workl(iq), ldq, workl(iw), ipntr, workd,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	606 & info )
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	607 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	608 c %--------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	609 c \| ido .ne. 99 implies use of reverse communication \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	610 c \| to compute operations involving OP or shifts. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	611 c %--------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	612 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	613 if (ido .eq. 3) iparam(8) = np
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	614 if (ido .ne. 99) go to 9000
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	615 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	616 iparam(3) = mxiter
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	617 iparam(5) = np
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	618 iparam(9) = nopx
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	619 iparam(10) = nbx
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	620 iparam(11) = nrorth
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	621 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	622 c %------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	623 c \| Exit if there was an informational \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	624 c \| error within ssaup2. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	625 c %------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	626 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	627 if (info .lt. 0) go to 9000
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	628 if (info .eq. 2) info = 3
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	629 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	630 if (msglvl .gt. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	631 call ivout (logfil, 1, mxiter, ndigit,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	632 & '_saupd: number of update iterations taken')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	633 call ivout (logfil, 1, np, ndigit,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	634 & '_saupd: number of "converged" Ritz values')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	635 call svout (logfil, np, workl(Ritz), ndigit,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	636 & '_saupd: final Ritz values')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	637 call svout (logfil, np, workl(Bounds), ndigit,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	638 & '_saupd: corresponding error bounds')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	639 end if
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	640 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	641 call arscnd (t1)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	642 tsaupd = t1 - t0
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	643 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	644 if (msglvl .gt. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	645 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	646 c %--------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	647 c \| Version Number & Version Date are defined in version.h \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	648 c %--------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	649 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	650 write (6,1000)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	651 write (6,1100) mxiter, nopx, nbx, nrorth, nitref, nrstrt,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	652 & tmvopx, tmvbx, tsaupd, tsaup2, tsaitr, titref,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	653 & tgetv0, tseigt, tsgets, tsapps, tsconv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	654 1000 format (//,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	655 & 5x, '==========================================',/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	656 & 5x, '= Symmetric implicit Arnoldi update code =',/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	657 & 5x, '= Version Number:', ' 2.4' , 19x, ' =',/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	658 & 5x, '= Version Date: ', ' 07/31/96' , 14x, ' =',/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	659 & 5x, '==========================================',/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	660 & 5x, '= Summary of timing statistics =',/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	661 & 5x, '==========================================',//)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	662 1100 format (
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	663 & 5x, 'Total number update iterations = ', i5,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	664 & 5x, 'Total number of OP*x operations = ', i5,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	665 & 5x, 'Total number of B*x operations = ', i5,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	666 & 5x, 'Total number of reorthogonalization steps = ', i5,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	667 & 5x, 'Total number of iterative refinement steps = ', i5,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	668 & 5x, 'Total number of restart steps = ', i5,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	669 & 5x, 'Total time in user OP*x operation = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	670 & 5x, 'Total time in user B*x operation = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	671 & 5x, 'Total time in Arnoldi update routine = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	672 & 5x, 'Total time in saup2 routine = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	673 & 5x, 'Total time in basic Arnoldi iteration loop = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	674 & 5x, 'Total time in reorthogonalization phase = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	675 & 5x, 'Total time in (re)start vector generation = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	676 & 5x, 'Total time in trid eigenvalue subproblem = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	677 & 5x, 'Total time in getting the shifts = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	678 & 5x, 'Total time in applying the shifts = ', f12.6,/
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	679 & 5x, 'Total time in convergence testing = ', f12.6)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	680 end if
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	681 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	682 9000 continue
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	683 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	684 return
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	685 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	686 c %---------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	687 c \| End of ssaupd \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	688 c %---------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	689 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	690 end

Mercurial > octave-nkf

annotate libcruft/arpack/src/ssaupd.f @ 12274:9f5d2ef078e8 release-3-4-x