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

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\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	2 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	3 c\Name: cnaup2
9f5d2ef078e8 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
9f5d2ef078e8 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\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	6 c Intermediate level interface called by cnaupd.
9f5d2ef078e8 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\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	9 c call cnaup2
9f5d2ef078e8 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 ( IDO, BMAT, N, WHICH, NEV, 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	11 c ISHIFT, MXITER, V, LDV, H, LDH, RITZ, 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	12 c Q, LDQ, WORKL, IPNTR, WORKD, RWORK, 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	13 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	14 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	15 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	16 c IDO, BMAT, N, WHICH, NEV, TOL, RESID: same as defined in cnaupd.
9f5d2ef078e8 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 MODE, ISHIFT, MXITER: see the definition of IPARAM in cnaupd.
9f5d2ef078e8 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 NP 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	20 c Contains the number of implicit shifts to apply during
9f5d2ef078e8 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 each 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	22 c If ISHIFT=1, NP is adjusted dynamically at each 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	23 c to accelerate convergence and prevent stagnation.
9f5d2ef078e8 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 This is also roughly equal to the number of matrix-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	25 c products (involving the operator OP) per 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	26 c The logic for adjusting is contained within the current
9f5d2ef078e8 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 subroutine.
9f5d2ef078e8 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 If ISHIFT=0, NP is the number of shifts the user needs
9f5d2ef078e8 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 to provide via reverse comunication. 0 < 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	30 c NP may be less than NCV-NEV since a leading block of the current
9f5d2ef078e8 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 upper Hessenberg matrix has split off and contains "unwanted"
9f5d2ef078e8 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 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	33 c Upon termination of the IRA iteration, NP contains 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	34 c of "converged" 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	35 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	36 c IUPD 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	37 c IUPD .EQ. 0: use explicit restart instead implicit update.
9f5d2ef078e8 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 IUPD .NE. 0: use implicit update.
9f5d2ef078e8 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
9f5d2ef078e8 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 V Complex N by (NEV+NP) array. (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	41 c The Arnoldi basis vectors are returned in the first 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	42 c columns of 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	43 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	44 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	45 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	46 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	47 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	48 c H Complex (NEV+NP) by (NEV+NP) 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	49 c H is used to store the generated upper Hessenberg 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	50 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	51 c LDH 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	52 c Leading dimension of H 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	53 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	54 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	55 c RITZ Complex array of length NEV+NP. (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	56 c RITZ(1:NEV) contains the computed Ritz values 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	57 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	58 c BOUNDS Complex array of length NEV+NP. (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	59 c BOUNDS(1:NEV) contain the error bounds corresponding 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	60 c the computed 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	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 Q Complex (NEV+NP) by (NEV+NP) array. (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	63 c Private (replicated) work array used to accumulate 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	64 c rotation in the shift application step.
9f5d2ef078e8 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
9f5d2ef078e8 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 LDQ 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	67 c Leading dimension of Q 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	68 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	69 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	70 c WORKL Complex work array of length at least
9f5d2ef078e8 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 (NEV+NP)*2 + 3(NEV+NP). (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	72 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	73 c the front end. It is used in shifts calculation, 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	74 c application and convergence 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	75 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	76 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	77 c IPNTR Integer array of length 3. (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	78 c Pointer to mark the starting locations in the WORKD 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	79 c vectors used by the 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	80 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	81 c IPNTR(1): pointer to the current operand vector 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	82 c IPNTR(2): pointer to the current result vector 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	83 c IPNTR(3): pointer to the vector B * X when used 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	84 c shift-and-invert mode. X is the current operand.
9f5d2ef078e8 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 -------------------------------------------------------------
9f5d2ef078e8 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
9f5d2ef078e8 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 WORKD Complex work array of length 3*N. (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	88 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	89 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	90 c as temporary workspace during the 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	91 c See Data Distribution Note in CNAUPD.
9f5d2ef078e8 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
9f5d2ef078e8 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 RWORK Real work array of length NEV+NP ( 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	94 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	95 c the front 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	96 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	97 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	98 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	99 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	100 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	101 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	102 c = 0: Normal 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	103 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	104 c All possible eigenvalues of OP has been found.
9f5d2ef078e8 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 NP returns the 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	106 c = 2: No shifts could be applied.
9f5d2ef078e8 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 = -8: Error return from LAPACK 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	108 c This should never happen.
9f5d2ef078e8 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 = -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	110 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	111 c Size that was built in returned in 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	112 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	113 c\EndDoc
9f5d2ef078e8 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
9f5d2ef078e8 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-----------------------------------------------------------------------
9f5d2ef078e8 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
9f5d2ef078e8 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\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	118 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	119 c\Local 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	120 c xxxxxx Complex
9f5d2ef078e8 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
9f5d2ef078e8 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\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	123 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	124 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	125 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	126 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	127 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	128 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	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\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	131 c cgetv0 ARPACK initial vector generation routine.
9f5d2ef078e8 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 cnaitr ARPACK Arnoldi factorization routine.
9f5d2ef078e8 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 cnapps ARPACK application of implicit shifts routine.
9f5d2ef078e8 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 cneigh ARPACK compute Ritz values and error bounds routine.
9f5d2ef078e8 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 cngets ARPACK reorder Ritz values and error bounds routine.
9f5d2ef078e8 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 csortc ARPACK sorting routine.
9f5d2ef078e8 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 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	138 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	139 c cmout ARPACK utility routine that prints matrices
9f5d2ef078e8 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 cvout 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	141 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	142 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	143 c slapy2 LAPACK routine to compute sqrt(x2+y2) carefully.
9f5d2ef078e8 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 ccopy Level 1 BLAS that copies one vector to another .
9f5d2ef078e8 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 cdotc Level 1 BLAS that computes the scalar product of two 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	146 c cswap Level 1 BLAS that swaps two 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	147 c scnrm2 Level 1 BLAS that computes the norm of a 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	148 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	149 c\Author
9f5d2ef078e8 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 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	151 c Richard Lehoucq CRPC / Rice Universitya
9f5d2ef078e8 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 Chao Yang 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	153 c Dept. of Computational &
9f5d2ef078e8 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 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	155 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	156 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	157 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	158 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	159 c FILE: naup2.F SID: 2.6 DATE OF SID: 06/01/00 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	160 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	161 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	162 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	163 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	164 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	165 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	166 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	167 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	168 subroutine cnaup2
9f5d2ef078e8 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 & ( ido, bmat, n, which, nev, 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	170 & ishift, mxiter, v, ldv, h, ldh, ritz, 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	171 & q, ldq, workl, ipntr, workd, rwork, 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	172 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	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 \| 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	175 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	176 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	177 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	178 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	179 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	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 \| 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	182 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	183 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	184 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	185 integer ido, info, ishift, iupd, mode, ldh, ldq, ldv, 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	186 & n, nev, 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	187 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	188 & 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	189 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	190 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	191 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	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
9f5d2ef078e8 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 integer ipntr(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	195 Complex
9f5d2ef078e8 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 & bounds(nev+np), h(ldh,nev+np), q(ldq,nev+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	197 & resid(n), ritz(nev+np), v(ldv,nev+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	198 & workd(3n), workl( (nev+np)(nev+np+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	199 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	200 & rwork(nev+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	201 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	202 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	203 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	204 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	205 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	206 Complex
9f5d2ef078e8 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 & 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	208 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	209 & rzero
9f5d2ef078e8 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 parameter (one = (1.0E+0, 0.0E+0) , zero = (0.0E+0, 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	211 & rzero = 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	212 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	213 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	214 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	215 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	216 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	217 logical cnorm , getv0, initv , update, ushift
9f5d2ef078e8 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 integer ierr , iter , kplusp, msglvl, nconv,
9f5d2ef078e8 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 & nevbef, nev0 , np0 , nptemp, 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	220 & 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	221 Complex
9f5d2ef078e8 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 & cmpnorm
9f5d2ef078e8 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 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	224 & rnorm , eps23, rtemp
9f5d2ef078e8 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 character wprime*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	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 save cnorm, getv0, initv , update, ushift,
9f5d2ef078e8 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 & rnorm, iter , kplusp, msglvl, nconv ,
9f5d2ef078e8 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 & nevbef, nev0 , np0 , eps23
9f5d2ef078e8 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
9f5d2ef078e8 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
9f5d2ef078e8 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 %-----------------------%
9f5d2ef078e8 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 \| Local 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	234 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	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 integer kp(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	237 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	238 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	239 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	240 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	241 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	242 external ccopy, cgetv0, cnaitr, cneigh, cngets, cnapps,
9f5d2ef078e8 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 & csortc, cswap, cmout, cvout, ivout, arscnd
9f5d2ef078e8 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
9f5d2ef078e8 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 %--------------------%
9f5d2ef078e8 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 \| 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	247 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	248 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	249 Complex
9f5d2ef078e8 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 & cdotc
9f5d2ef078e8 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 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	252 & scnrm2, slamch, slapy2
9f5d2ef078e8 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 external cdotc, scnrm2, slamch, slapy2
9f5d2ef078e8 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
9f5d2ef078e8 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 %---------------------%
9f5d2ef078e8 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 \| Intrinsic 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	257 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	258 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	259 intrinsic aimag, real , min, max
9f5d2ef078e8 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
9f5d2ef078e8 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 %-----------------------%
9f5d2ef078e8 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 \| 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	263 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	264 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	265 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	266 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	267 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	268 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	269 msglvl = mcaup2
9f5d2ef078e8 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
9f5d2ef078e8 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 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	272 np0 = 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	273 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	274 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	275 c \| kplusp is the bound on the largest \|
9f5d2ef078e8 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 \| Lanczos factorization built. \|
9f5d2ef078e8 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 \| nconv is the current number 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	278 c \| "converged" 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	279 c \| iter is the counter on the current \|
9f5d2ef078e8 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 \| iteration step. \|
9f5d2ef078e8 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 %-------------------------------------%
9f5d2ef078e8 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
9f5d2ef078e8 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 kplusp = nev + 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	284 nconv = 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	285 iter = 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	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 %---------------------------------%
9f5d2ef078e8 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 \| Get machine dependent constant. \|
9f5d2ef078e8 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 %---------------------------------%
9f5d2ef078e8 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 eps23 = slamch('Epsilon-Machine')
9f5d2ef078e8 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 eps23 = eps23**(2.0E+0 / 3.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	293 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	294 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	295 c \| Set flags for computing the first 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	296 c \| steps of the 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	297 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	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 getv0 = .true.
9f5d2ef078e8 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 update = .false.
9f5d2ef078e8 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 ushift = .false.
9f5d2ef078e8 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 cnorm = .false.
9f5d2ef078e8 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
9f5d2ef078e8 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 if (info .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	305 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	306 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	307 c \| User provides 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	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
9f5d2ef078e8 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 initv = .true.
9f5d2ef078e8 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 info = 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	312 else
9f5d2ef078e8 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 initv = .false.
9f5d2ef078e8 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 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	315 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	316 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	317 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	318 c \| Get a possibly random starting vector 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	319 c \| force it into the range of 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	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 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	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 if (getv0) 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	325 call cgetv0 (ido, bmat, 1, initv, n, 1, v, ldv, resid, rnorm,
9f5d2ef078e8 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 & ipntr, workd, 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	327 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	328 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	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 if (rnorm .eq. rzero) 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	331 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	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 \| The initial vector is zero. 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	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
9f5d2ef078e8 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 info = -9
9f5d2ef078e8 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 go to 1100
9f5d2ef078e8 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 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	339 getv0 = .false.
9f5d2ef078e8 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 ido = 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	341 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	342 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	343 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	344 c \| Back from 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	345 c \| continue with update step \|
9f5d2ef078e8 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
9f5d2ef078e8 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 if (update) go to 20
9f5d2ef078e8 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
9f5d2ef078e8 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 \| Back from computing user specified 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	352 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	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 if (ushift) go to 50
9f5d2ef078e8 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
9f5d2ef078e8 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 %-------------------------------------%
9f5d2ef078e8 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 \| Back from computing residual norm \|
9f5d2ef078e8 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 \| at the end of the current 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	359 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	360 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	361 if (cnorm) go to 100
9f5d2ef078e8 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
9f5d2ef078e8 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 %----------------------------------------------------------%
9f5d2ef078e8 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 \| Compute the first NEV steps of the 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	365 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	366 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	367 call cnaitr (ido, bmat, n, 0, nev, mode, resid, rnorm, v, ldv,
9f5d2ef078e8 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 & h, ldh, ipntr, workd, 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	369 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	370 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	371 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	372 if (info .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	373 np = 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	374 mxiter = iter
9f5d2ef078e8 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 info = -9999
9f5d2ef078e8 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 go to 1200
9f5d2ef078e8 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 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	378 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	379 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	380 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	381 c \| M A I N ARNOLDI I T E R A T I O N L O O P \|
9f5d2ef078e8 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 \| Each iteration implicitly restarts the 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	383 c \| factorization in place. \|
9f5d2ef078e8 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 \| \|
9f5d2ef078e8 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 %--------------------------------------------------------------%
9f5d2ef078e8 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 1000 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	388 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	389 iter = iter + 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	390 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	391 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	392 call ivout (logfil, 1, iter, 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	393 & '_naup2: ** Start of major iteration 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	394 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	395 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	396 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	397 c \| Compute NP additional steps of the 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	398 c \| Adjust NP since NEV might have been updated by last 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	399 c \| to the shift application routine cnapps. \|
9f5d2ef078e8 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
9f5d2ef078e8 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 np = kplusp - 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	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 if (msglvl .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	405 call ivout (logfil, 1, nev, 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	406 & '_naup2: The length of the current 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	407 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	408 & '_naup2: Extend the Arnoldi factorization 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	409 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	410 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	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 \| Compute NP additional steps of the 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	413 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	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 ido = 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	416 20 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	417 update = .true.
9f5d2ef078e8 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 call cnaitr(ido, bmat, n, nev, np, mode, resid, rnorm,
9f5d2ef078e8 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 & v , ldv , h, ldh, ipntr, workd, 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	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 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	423 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	424 if (info .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	425 np = 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	426 mxiter = iter
9f5d2ef078e8 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 info = -9999
9f5d2ef078e8 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 go to 1200
9f5d2ef078e8 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 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	430 update = .false.
9f5d2ef078e8 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 if (msglvl .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	433 call svout (logfil, 1, rnorm, 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	434 & '_naup2: Corresponding B-norm of the residual')
9f5d2ef078e8 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 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	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 %--------------------------------------------------------%
9f5d2ef078e8 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 \| Compute the eigenvalues and 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	439 c \| of the current upper Hessenberg 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	440 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	441 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	442 call cneigh (rnorm, kplusp, h, ldh, ritz, 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	443 & q, ldq, workl, rwork, 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	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 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	446 info = -8
9f5d2ef078e8 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 go to 1200
9f5d2ef078e8 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 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	449 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	450 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	451 c \| Select the wanted Ritz values and their 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	452 c \| to be used in the convergence test. \|
9f5d2ef078e8 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 c \| The wanted part of the spectrum and corresponding \|
9f5d2ef078e8 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 \| error bounds are in the last NEV loc. of 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	455 c \| and BOUNDS respectively. \|
9f5d2ef078e8 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 %---------------------------------------------------%
9f5d2ef078e8 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 nev = nev0
9f5d2ef078e8 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 np = np0
9f5d2ef078e8 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 \| Make a copy of Ritz values and the corresponding \|
9f5d2ef078e8 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 \| Ritz estimates obtained from cneigh. \|
9f5d2ef078e8 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 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	466 call ccopy(kplusp,ritz,1,workl(kplusp**2+1),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	467 call ccopy(kplusp,bounds,1,workl(kplusp**2+kplusp+1),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	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 \| Select the wanted Ritz values and their 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	471 c \| to be used in the convergence test. \|
9f5d2ef078e8 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 \| The wanted part of the spectrum and corresponding \|
9f5d2ef078e8 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 c \| bounds are in the last NEV loc. of 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	474 c \| BOUNDS respectively. \|
9f5d2ef078e8 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
9f5d2ef078e8 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 call cngets (ishift, which, nev, np, ritz, 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	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 c \| Convergence test: currently we use the following criteria. \|
9f5d2ef078e8 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 c \| The relative accuracy of a Ritz value is considered \|
9f5d2ef078e8 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 c \| acceptable 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	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 c \| error_bounds(i) .le. tol*max(eps23, magnitude_of_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	485 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	486 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	487 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	488 nconv = 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	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 do 25 i = 1, 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	491 rtemp = max( eps23, slapy2( real (ritz(np+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	492 & aimag(ritz(np+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	493 if ( slapy2(real (bounds(np+i)),aimag(bounds(np+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	494 & .le. tol*rtemp ) 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	495 nconv = nconv + 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	496 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	497 25 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	498 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	499 if (msglvl .gt. 2) 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	500 kp(1) = nev

12274

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