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

annotate libcruft/arpack/src/dneupd.f @ 12277:3d38b4916cbf release-3-4-x

avoid memory overrun in ARPACK

author	David Bateman <dbateman@free.fr>
date	Fri, 28 Jan 2011 14:05:34 -0500
parents	9f5d2ef078e8
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: dneupd
9f5d2ef078e8 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
9f5d2ef078e8 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 This subroutine returns the converged approximations to 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	8 c of Az = lambdaB*z and (optionally):
9f5d2ef078e8 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
9f5d2ef078e8 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 (1) The corresponding approximate eigenvectors;
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	11 c
9f5d2ef078e8 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 (2) An orthonormal basis for the associated approximate
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	13 c invariant subspace;
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	14 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	15 c (3) Both.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	16 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	17 c There is negligible additional cost to obtain eigenvectors. An orthonormal
9f5d2ef078e8 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 basis is always computed. There is an additional storage cost of n*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	19 c if both are requested (in this case a separate array Z must be supplied).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	20 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	21 c The approximate eigenvalues and eigenvectors of Az = lambdaB*z
9f5d2ef078e8 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 are derived from approximate eigenvalues and eigenvectors 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	23 c of the linear operator OP prescribed by the MODE selection 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	24 c call to DNAUPD . DNAUPD must be called before this routine is 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	25 c These approximate eigenvalues and vectors are commonly called 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	26 c values and Ritz vectors respectively. They are referred to as such
9f5d2ef078e8 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 in the comments that follow. The computed orthonormal basis for the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	28 c invariant subspace corresponding to these Ritz values is referred to as a
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	29 c Schur basis.
9f5d2ef078e8 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
9f5d2ef078e8 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 See documentation in the header of the subroutine DNAUPD 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	32 c definition of OP as well as other terms and the relation of computed
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	33 c Ritz values and Ritz vectors of OP with respect to the given problem
9f5d2ef078e8 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 Az = lambdaB*z. For a brief description, see definitions 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	35 c IPARAM(7), MODE and WHICH in the documentation of DNAUPD .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	36 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	37 c\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	38 c call dneupd
9f5d2ef078e8 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 ( RVEC, HOWMNY, SELECT, DR, DI, Z, LDZ, SIGMAR, SIGMAI, WORKEV, BMAT,
9f5d2ef078e8 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 N, WHICH, NEV, TOL, RESID, NCV, V, LDV, IPARAM, IPNTR, WORKD, WORKL,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	41 c LWORKL, INFO )
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	42 c
9f5d2ef078e8 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\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	44 c RVEC LOGICAL (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 Specifies whether a basis for the invariant subspace 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	46 c to the converged Ritz value approximations for the eigenproblem
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	47 c Az = lambdaB*z is computed.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	48 c
9f5d2ef078e8 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 RVEC = .FALSE. Compute Ritz values only.
9f5d2ef078e8 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 RVEC = .TRUE. Compute the Ritz vectors or Schur 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	52 c See Remarks below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	53 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	54 c HOWMNY Character*1 (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	55 c Specifies the form of the basis for the invariant subspace
9f5d2ef078e8 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 corresponding to the converged Ritz values that is to be computed.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 = 'A': Compute NEV Ritz vectors;
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	59 c = 'P': Compute NEV Schur 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	60 c = 'S': compute some of the Ritz vectors, specified
9f5d2ef078e8 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 by the logical array SELECT.
9f5d2ef078e8 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
9f5d2ef078e8 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 SELECT Logical array of dimension NCV. (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	64 c If HOWMNY = 'S', SELECT specifies the Ritz vectors to be
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	65 c computed. To select the Ritz vector corresponding to a
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	66 c Ritz value (DR(j), DI(j)), SELECT(j) must be set to .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	67 c If HOWMNY = 'A' or 'P', SELECT is used as internal workspace.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	68 c
9f5d2ef078e8 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 DR Double precision array of dimension NEV+1. (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	70 c If IPARAM(7) = 1,2 or 3 and SIGMAI=0.0 then on exit: DR contains
9f5d2ef078e8 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 the real part of the Ritz approximations to the eigenvalues of
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	72 c Az = lambdaB*z.
9f5d2ef078e8 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 If IPARAM(7) = 3, 4 and SIGMAI is not equal to zero, then on 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	74 c DR contains the real part of the Ritz values of OP computed 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	75 c DNAUPD . A further computation must be performed by the user
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	76 c to transform the Ritz values computed for OP by DNAUPD to those
9f5d2ef078e8 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 of the original system Az = lambdaB*z. See remark 3 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	78 c
9f5d2ef078e8 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 DI Double precision array of dimension NEV+1. (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	80 c On exit, DI contains the imaginary part of the Ritz value
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	81 c approximations to the eigenvalues of Az = lambdaB*z associated
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	82 c with DR.
9f5d2ef078e8 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
9f5d2ef078e8 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 NOTE: When Ritz values are complex, they will come in 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	85 c conjugate pairs. If eigenvectors are requested, 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	86 c corresponding Ritz vectors will also come in conjugate
9f5d2ef078e8 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 pairs and the real and imaginary parts of these are
9f5d2ef078e8 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 represented in two consecutive columns of the array Z
9f5d2ef078e8 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 (see below).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	90 c
9f5d2ef078e8 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 Z Double precision N by NEV+1 array if RVEC = .TRUE. and HOWMNY = 'A'. (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	92 c On exit, if RVEC = .TRUE. and HOWMNY = 'A', then the columns 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	93 c Z represent approximate eigenvectors (Ritz vectors) 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	94 c to the NCONV=IPARAM(5) Ritz values for eigensystem
9f5d2ef078e8 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 Az = lambdaB*z.
9f5d2ef078e8 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 The complex Ritz vector associated with the Ritz value
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	98 c with positive imaginary part is stored in two consecutive
9f5d2ef078e8 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 columns. The first column holds the real part of the 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	100 c vector and the second column holds the imaginary part. The
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	101 c Ritz vector associated with the Ritz value with negative
9f5d2ef078e8 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 imaginary part is simply the complex conjugate of the Ritz 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	103 c associated with the positive imaginary part.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	104 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	105 c If RVEC = .FALSE. or HOWMNY = 'P', then Z is not referenced.
9f5d2ef078e8 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
9f5d2ef078e8 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 NOTE: If if RVEC = .TRUE. and a Schur basis is not required,
9f5d2ef078e8 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 the array Z may be set equal to first NEV+1 columns of 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	109 c basis array V computed by DNAUPD . In this case the Arnoldi basis
9f5d2ef078e8 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 will be destroyed and overwritten with the eigenvector basis.
9f5d2ef078e8 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
9f5d2ef078e8 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 LDZ 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	113 c The leading dimension of the array Z. If Ritz vectors are
9f5d2ef078e8 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 desired, then LDZ >= max( 1, N ). In any case, LDZ >= 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	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 SIGMAR Double precision (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	117 c If IPARAM(7) = 3 or 4, represents the real part of the shift.
9f5d2ef078e8 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 Not referenced if IPARAM(7) = 1 or 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	119 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	120 c SIGMAI Double precision (INPUT)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	121 c If IPARAM(7) = 3 or 4, represents the imaginary part of the shift.
9f5d2ef078e8 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 Not referenced if IPARAM(7) = 1 or 2. See remark 3 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	123 c
9f5d2ef078e8 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 WORKEV Double precision work array of dimension 3*NCV. (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	125 c
9f5d2ef078e8 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 ** The remaining arguments MUST be the same as for the **
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	127 c ** call to DNAUPD that was just completed. **
9f5d2ef078e8 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
9f5d2ef078e8 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 NOTE: The remaining 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	130 c
9f5d2ef078e8 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 BMAT, N, WHICH, NEV, TOL, RESID, NCV, V, LDV, IPARAM, IPNTR,
9f5d2ef078e8 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 WORKD, WORKL, LWORKL, INFO
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	133 c
9f5d2ef078e8 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 must be passed directly to DNEUPD following the 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	135 c to DNAUPD . These arguments MUST NOT BE MODIFIED between
9f5d2ef078e8 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 the the last call to DNAUPD and the call to DNEUPD .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	137 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	138 c Three of these parameters (V, WORKL, INFO) are also output 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	139 c
9f5d2ef078e8 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 V Double precision N by NCV 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	141 c
9f5d2ef078e8 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 Upon INPUT: the NCV columns of V contain the Arnoldi basis
9f5d2ef078e8 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 vectors for OP as constructed by DNAUPD .
9f5d2ef078e8 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
9f5d2ef078e8 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 Upon OUTPUT: If RVEC = .TRUE. the first NCONV=IPARAM(5) columns
9f5d2ef078e8 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 contain approximate Schur vectors that span 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	147 c desired invariant subspace. See Remark 2 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 NOTE: If the array Z has been set equal to first NEV+1 columns
9f5d2ef078e8 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 of the array V and RVEC=.TRUE. and HOWMNY= 'A', then 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	151 c Arnoldi basis held by V has been overwritten by the desired
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	152 c Ritz vectors. If a separate array Z has been passed 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	153 c the first NCONV=IPARAM(5) columns of V will contain approximate
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	154 c Schur vectors that span the desired invariant subspace.
9f5d2ef078e8 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
9f5d2ef078e8 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 WORKL Double precision work array of length LWORKL. (OUTPUT/WORKSPACE)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	157 c WORKL(1:ncvncv+3ncv) contains information obtained 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	158 c dnaupd . They are not changed by dneupd .
9f5d2ef078e8 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 WORKL(ncvncv+3ncv+1:3ncvncv+6*ncv) holds 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	160 c real and imaginary part of the untransformed 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	161 c the upper quasi-triangular matrix for H, and 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	162 c associated matrix representation of the invariant subspace for 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	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 Note: IPNTR(9:13) contains the pointer into WORKL for addresses
9f5d2ef078e8 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 of the above information computed by dneupd .
9f5d2ef078e8 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 IPNTR(9): pointer to the real part of the NCV RITZ values of the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	168 c original system.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	169 c IPNTR(10): pointer to the imaginary part of the NCV RITZ values 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	170 c the original system.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	171 c IPNTR(11): pointer to the NCV corresponding error bounds.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	172 c IPNTR(12): pointer to the NCV by NCV upper quasi-triangular
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	173 c Schur matrix for 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	174 c IPNTR(13): pointer to the NCV by NCV matrix of eigenvectors
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	175 c of the upper Hessenberg matrix H. Only referenced by
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	176 c dneupd if RVEC = .TRUE. See Remark 2 below.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	177 c -------------------------------------------------------------
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	178 c
9f5d2ef078e8 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 INFO Integer. (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	180 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	181 c
9f5d2ef078e8 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 = 0: Normal exit.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 c = 1: The Schur form computed by LAPACK routine dlahqr
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	185 c could not be reordered by LAPACK routine dtrsen .
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	186 c Re-enter subroutine dneupd with IPARAM(5)=NCV 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	187 c increase the size of the arrays DR and DI to have
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	188 c dimension at least dimension NCV and allocate at least NCV
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	189 c columns for Z. NOTE: Not necessary if Z and V share
9f5d2ef078e8 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 the same space. Please notify the authors if this error
9f5d2ef078e8 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 occurs.
9f5d2ef078e8 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 = -1: N must be positive.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	194 c = -2: NEV must be positive.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	195 c = -3: NCV-NEV >= 2 and less than or equal to N.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	196 c = -5: WHICH must be one of 'LM', 'SM', 'LR', 'SR', 'LI', 'SI'
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	197 c = -6: BMAT must be one of 'I' or 'G'.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	198 c = -7: Length of private work WORKL array is not sufficient.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	199 c = -8: Error return from calculation of a real Schur form.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	200 c Informational error from LAPACK routine dlahqr .
9f5d2ef078e8 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 = -9: Error return from calculation of eigenvectors.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	202 c Informational error from LAPACK routine dtrevc .
9f5d2ef078e8 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 = -10: IPARAM(7) must be 1,2,3,4.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	204 c = -11: IPARAM(7) = 1 and BMAT = 'G' are incompatible.
9f5d2ef078e8 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 = -12: HOWMNY = 'S' not yet implemented
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	206 c = -13: HOWMNY must be one of 'A' or 'P' if RVEC = .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	207 c = -14: DNAUPD did not find any eigenvalues to sufficient
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	208 c accuracy.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	209 c = -15: DNEUPD got a different count of the number of converged
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	210 c Ritz values than DNAUPD got. This indicates the user
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	211 c probably made an error in passing data from DNAUPD 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	212 c DNEUPD or that the data was modified before entering
9f5d2ef078e8 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 DNEUPD
9f5d2ef078e8 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
9f5d2ef078e8 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\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	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 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	218 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	219 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	220 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	221 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	222 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	223 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	224 c 3. B.N. Parlett & Y. Saad, "Complex Shift and Invert Strategies 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	225 c Real Matrices", Linear Algebra and its Applications, vol 88/89,
9f5d2ef078e8 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 pp 575-595, (1987).
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	227 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	228 c\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	229 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	230 c dmout 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	231 c dvout 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	232 c dgeqr2 LAPACK routine that computes the QR factorization 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	233 c a 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	234 c dlacpy LAPACK matrix copy 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	235 c dlahqr LAPACK routine to compute the real Schur form of an
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	236 c 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	237 c dlamch 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	238 c dlapy2 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	239 c dlaset LAPACK matrix initialization 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	240 c dorm2r LAPACK routine that applies an orthogonal matrix 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	241 c factored form.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	242 c dtrevc LAPACK routine to compute the eigenvectors of a 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	243 c in upper quasi-triangular form.
9f5d2ef078e8 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 dtrsen LAPACK routine that re-orders the Schur form.
9f5d2ef078e8 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 dtrmm Level 3 BLAS matrix times an upper triangular 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	246 c dger Level 2 BLAS rank one update to a 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	247 c dcopy 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	248 c ddot 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	249 c dnrm2 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	250 c dscal Level 1 BLAS that scales 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	251 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	252 c\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	253 c
9f5d2ef078e8 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 1. Currently only HOWMNY = 'A' and 'P' are implemented.
9f5d2ef078e8 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 Let trans(X) denote the transpose of 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	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 2. Schur vectors are an orthogonal representation for the basis 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	259 c Ritz vectors. Thus, their numerical properties are often superior.
9f5d2ef078e8 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 If RVEC = .TRUE. then the relationship
9f5d2ef078e8 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 A * V(:,1:IPARAM(5)) = V(:,1:IPARAM(5)) * T, 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	262 c trans(V(:,1:IPARAM(5))) * V(:,1:IPARAM(5)) = I are approximately
9f5d2ef078e8 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 satisfied. Here T is the leading submatrix of order IPARAM(5) of the
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	264 c real upper quasi-triangular matrix stored workl(ipntr(12)). That is,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	265 c T is block upper triangular with 1-by-1 and 2-by-2 diagonal blocks;
9f5d2ef078e8 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 each 2-by-2 diagonal block has its diagonal elements equal and its
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	267 c off-diagonal elements of opposite sign. Corresponding to each 2-by-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	268 c diagonal block is a complex conjugate pair of Ritz values. The 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	269 c Ritz values are stored on the diagonal of T.
9f5d2ef078e8 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 c 3. If IPARAM(7) = 3 or 4 and SIGMAI is not equal zero, then the user must
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	272 c form the IPARAM(5) Rayleigh quotients in order to transform the 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	273 c values computed by DNAUPD for OP to those of Az = lambdaB*z.
9f5d2ef078e8 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 Set RVEC = .true. and HOWMNY = 'A', 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	275 c compute
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	276 c trans(Z(:,I)) * A * Z(:,I) if DI(I) = 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	277 c If DI(I) is not equal to zero and DI(I+1) = - D(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	278 c then the desired real and imaginary parts of the Ritz value are
9f5d2ef078e8 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 trans(Z(:,I)) * A * Z(:,I) + trans(Z(:,I+1)) * A * Z(:,I+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	280 c trans(Z(:,I)) * A * Z(:,I+1) - trans(Z(:,I+1)) * A * Z(:,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	281 c 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	282 c Another possibility is to set RVEC = .true. and HOWMNY = 'P' 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	283 c compute trans(V(:,1:IPARAM(5))) * A * V(:,1:IPARAM(5)) and then an upper
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	284 c quasi-triangular matrix of order IPARAM(5) is computed. See remark
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	285 c 2 above.
9f5d2ef078e8 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\Authors
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	288 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	289 c Richard Lehoucq CRPC / Rice University
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	290 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	291 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	292 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	293 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	294 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	295 c
9f5d2ef078e8 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\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	297 c FILE: neupd.F SID: 2.7 DATE OF SID: 09/20/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	298 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	299 c\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	300 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	301 c-----------------------------------------------------------------------
9f5d2ef078e8 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 subroutine dneupd (rvec , howmny, select, dr , di,
9f5d2ef078e8 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 & z , ldz , sigmar, sigmai, workev,
9f5d2ef078e8 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 & bmat , n , which , nev , 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	305 & resid, ncv , v , ldv , iparam,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	306 & ipntr, workd , workl , lworkl, info)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	307 c
9f5d2ef078e8 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 \| 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	310 c %----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	311 c
9f5d2ef078e8 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 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	313 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	314 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	315 c %------------------%
9f5d2ef078e8 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 \| 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	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
9f5d2ef078e8 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 character bmat, howmny, which*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	320 logical rvec
9f5d2ef078e8 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 integer info, ldz, ldv, lworkl, n, ncv, nev
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	322 Double precision
9f5d2ef078e8 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 & sigmar, sigmai, 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	324 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	325 c %-----------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	326 c \| 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	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 c
9f5d2ef078e8 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 integer iparam(11), ipntr(14)
9f5d2ef078e8 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 logical select(ncv)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	331 Double precision
9f5d2ef078e8 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 & dr(nev+1) , di(nev+1), resid(n) ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	333 & v(ldv,ncv) , z(ldz,) , workd(3n),
9f5d2ef078e8 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 & workl(lworkl), workev(3*ncv)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 c %------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	337 c \| 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	338 c %------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	339 c
9f5d2ef078e8 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 Double precision
9f5d2ef078e8 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 & 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	342 parameter (one = 1.0D+0 , zero = 0.0D+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	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 %---------------%
9f5d2ef078e8 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 \| 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	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 character type*6
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	349 integer bounds, ierr , ih , ihbds ,
9f5d2ef078e8 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 & iheigr, iheigi, iconj , 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	351 & invsub, iuptri, iwev , iwork(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	352 & j , k , ldh , ldq ,
9f5d2ef078e8 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 & mode , msglvl, outncv, ritzr ,
9f5d2ef078e8 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 & ritzi , wri , wrr , irr ,
9f5d2ef078e8 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 & iri , ibd , ishift, numcnv ,
12277 3d38b4916cbf avoid memory overrun in ARPACK David Bateman <dbateman@free.fr> parents: 12274 diff changeset	356 & np , jj , nconv2
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	357 logical reord
9f5d2ef078e8 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 Double precision
9f5d2ef078e8 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 & conds , rnorm, sep , temp,
9f5d2ef078e8 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 & vl(1,1), temp1, 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	361 c
9f5d2ef078e8 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 \| 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	364 c %----------------------%
9f5d2ef078e8 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 external dcopy , dger , dgeqr2 , dlacpy ,
9f5d2ef078e8 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 & dlahqr , dlaset , dmout , dorm2r ,
9f5d2ef078e8 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 & dtrevc , dtrmm , dtrsen , dscal ,
9f5d2ef078e8 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 & dvout , ivout
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	370 c
9f5d2ef078e8 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 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	373 c %--------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	374 c
9f5d2ef078e8 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 Double precision
9f5d2ef078e8 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 & dlapy2 , dnrm2 , dlamch , ddot
9f5d2ef078e8 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 external dlapy2 , dnrm2 , dlamch , ddot
9f5d2ef078e8 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 \| 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	381 c %---------------------%
9f5d2ef078e8 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
9f5d2ef078e8 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 intrinsic abs, min, sqrt
9f5d2ef078e8 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 \| 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	387 c %-----------------------%
9f5d2ef078e8 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 c %------------------------%
9f5d2ef078e8 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 \| Set default parameters \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	391 c %------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	392 c
9f5d2ef078e8 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 msglvl = mneupd
9f5d2ef078e8 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 mode = iparam(7)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	395 nconv = iparam(5)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	396 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	397 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	398 c %---------------------------------%
9f5d2ef078e8 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 \| 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	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 eps23 = dlamch ('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	403 eps23 = eps23**(2.0D+0 / 3.0D+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	404 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	405 c %--------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	406 c \| Quick 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	407 c %--------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	408 c
9f5d2ef078e8 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 ierr = 0
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 if (nconv .le. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	412 ierr = -14
9f5d2ef078e8 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 else if (n .le. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	414 ierr = -1
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	415 else if (nev .le. 0) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	416 ierr = -2
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	417 else if (ncv .le. nev+1 .or. ncv .gt. n) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	418 ierr = -3
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	419 else if (which .ne. 'LM' .and.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	420 & which .ne. 'SM' .and.
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	421 & which .ne. 'LR' .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	422 & which .ne. 'SR' .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	423 & which .ne. 'LI' .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	424 & which .ne. 'SI') 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 ierr = -5
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	426 else if (bmat .ne. 'I' .and. bmat .ne. 'G') then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	427 ierr = -6
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	428 else if (lworkl .lt. 3ncv2 + 6ncv) 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	429 ierr = -7
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	430 else if ( (howmny .ne. 'A' .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	431 & howmny .ne. 'P' .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	432 & howmny .ne. 'S') .and. rvec ) 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 ierr = -13
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	434 else if (howmny .eq. 'S' ) 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	435 ierr = -12
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	436 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	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 if (mode .eq. 1 .or. mode .eq. 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	439 type = 'REGULR'
9f5d2ef078e8 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 else if (mode .eq. 3 .and. sigmai .eq. zero) 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	441 type = 'SHIFTI'
9f5d2ef078e8 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 else if (mode .eq. 3 ) 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	443 type = 'REALPT'
9f5d2ef078e8 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 else if (mode .eq. 4 ) 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	445 type = 'IMAGPT'
9f5d2ef078e8 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 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	447 ierr = -10
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 if (mode .eq. 1 .and. bmat .eq. 'G') ierr = -11
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 %------------%
9f5d2ef078e8 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 \| 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	453 c %------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	454 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	455 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	456 info = ierr
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	457 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	458 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	459 c
9f5d2ef078e8 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 \| Pointer into WORKL for address of H, RITZ, BOUNDS, Q \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	462 c \| etc... and the remaining workspace. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	463 c \| Also update pointer to be used on output. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	464 c \| Memory is laid out as follows: \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	465 c \| workl(1:ncv*ncv) := generated 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	466 c \| workl(ncvncv+1:ncvncv+2*ncv) := real and imaginary \|
9f5d2ef078e8 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 c \| parts of 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	468 c \| workl(ncvncv+2ncv+1:ncvncv+3ncv) := 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	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
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	471 c %-----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	472 c \| The following is used and set by DNEUPD . \|
9f5d2ef078e8 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 \| workl(ncvncv+3ncv+1:ncvncv+4ncv) := The untransformed \|
9f5d2ef078e8 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 \| real part of the 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	475 c \| workl(ncvncv+4ncv+1:ncvncv+5ncv) := The untransformed \|
9f5d2ef078e8 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 \| imaginary part of the 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	477 c \| workl(ncvncv+5ncv+1:ncvncv+6ncv) := The untransformed \|
9f5d2ef078e8 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 \| error bounds of the 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	479 c \| workl(ncvncv+6ncv+1:2ncvncv+6*ncv) := Holds the upper \|
9f5d2ef078e8 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 \| quasi-triangular matrix for 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	481 c \| workl(2ncvncv+6ncv+1: 3ncvncv+6ncv) := Holds 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	482 c \| associated matrix representation of the invariant \|
9f5d2ef078e8 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 \| subspace for 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	484 c \| GRAND total of NCV * ( 3 * NCV + 6 ) locations. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	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 ih = ipntr(5)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	488 ritzr = ipntr(6)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	489 ritzi = ipntr(7)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	490 bounds = ipntr(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	491 ldh = ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	492 ldq = ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	493 iheigr = bounds + ldh
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	494 iheigi = iheigr + ldh
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	495 ihbds = iheigi + ldh
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	496 iuptri = ihbds + ldh
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	497 invsub = iuptri + ldh*ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	498 ipntr(9) = iheigr
9f5d2ef078e8 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 ipntr(10) = iheigi
9f5d2ef078e8 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 ipntr(11) = ihbds
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	501 ipntr(12) = iuptri
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	502 ipntr(13) = invsub
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	503 wrr = 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	504 wri = ncv + 1
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	505 iwev = wri + ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	506 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	507 c %-----------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	508 c \| irr points to the REAL part of the 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	509 c \| values computed by _neigh before \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	510 c \| exiting _naup2. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	511 c \| iri points to the IMAGINARY part of the \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	512 c \| Ritz values computed by _neigh \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	513 c \| before exiting _naup2. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	514 c \| ibd points to the Ritz estimates \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	515 c \| computed by _neigh before exiting \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	516 c \| _naup2. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	517 c %-----------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	518 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	519 irr = ipntr(14)+ncv*ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	520 iri = irr+ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	521 ibd = iri+ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	522 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	523 c %------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	524 c \| RNORM is B-norm of the RESID(1:N). \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	525 c %------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	526 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	527 rnorm = workl(ih+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	528 workl(ih+2) = 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	529 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	530 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	531 call dvout (logfil, ncv, workl(irr), 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	532 & '_neupd: Real part of Ritz values passed in from _NAUPD.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	533 call dvout (logfil, ncv, workl(iri), 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	534 & '_neupd: Imag part of Ritz values passed in from _NAUPD.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	535 call dvout (logfil, ncv, workl(ibd), 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	536 & '_neupd: Ritz estimates passed in from _NAUPD.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	537 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	538 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	539 if (rvec) then
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	540 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	541 reord = .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	542 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	543 c %---------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	544 c \| Use the temporary bounds array to store indices \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	545 c \| These will be used to mark the select array later \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	546 c %---------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	547 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	548 do 10 j = 1,ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	549 workl(bounds+j-1) = 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	550 select(j) = .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	551 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	552 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	553 c %-------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	554 c \| Select the wanted Ritz values. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	555 c \| Sort the Ritz values so that the \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	556 c \| wanted ones appear at the tailing \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	557 c \| NEV positions of workl(irr) 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	558 c \| workl(iri). Move 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	559 c \| error estimates in workl(bound) \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	560 c \| accordingly. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	561 c %-------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	562 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	563 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	564 ishift = 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	565 call dngets (ishift , which , nev ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	566 & np , workl(irr), workl(iri),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	567 & workl(bounds), workl , workl(np+1))
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	568 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	569 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	570 call dvout (logfil, ncv, workl(irr), 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	571 & '_neupd: Real part of Ritz values after calling _NGETS.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	572 call dvout (logfil, ncv, workl(iri), 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	573 & '_neupd: Imag part of Ritz values after calling _NGETS.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	574 call dvout (logfil, ncv, workl(bounds), ndigit,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	575 & '_neupd: Ritz value indices after calling _NGETS.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	576 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	577 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	578 c %-----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	579 c \| Record indices of the 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	580 c \| Mark the select array for possible reordering \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	581 c %-----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	582 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	583 numcnv = 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	584 do 11 j = 1,ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	585 temp1 = max(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	586 & dlapy2 ( workl(irr+ncv-j), workl(iri+ncv-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	587 jj = workl(bounds + ncv - 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	588 if (numcnv .lt. nconv .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	589 & workl(ibd+jj-1) .le. tol*temp1) 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	590 select(jj) = .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	591 numcnv = numcnv + 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	592 if (jj .gt. nev) reord = .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	593 endif
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	594 11 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	595 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	596 c %-----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	597 c \| Check the count (numcnv) of converged Ritz values with \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	598 c \| the number (nconv) reported by dnaupd. If these two \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	599 c \| are different then there has probably been an error \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	600 c \| caused by incorrect passing of the dnaupd data. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	601 c %-----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	602 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	603 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	604 call ivout(logfil, 1, numcnv, 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	605 & '_neupd: Number of specified 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	606 call ivout(logfil, 1, nconv, 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	607 & '_neupd: Number of "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	608 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	609 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	610 if (numcnv .ne. nconv) 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	611 info = -15
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	612 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	613 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	614 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	615 c %-----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	616 c \| Call LAPACK routine dlahqr to compute the real Schur form \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	617 c \| of the upper Hessenberg matrix returned by DNAUPD . \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	618 c \| Make a copy of the 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	619 c \| Initialize the Schur vector matrix Q to the identity. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	620 c %-----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	621 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	622 call dcopy (ldh*ncv, workl(ih), 1, workl(iuptri), 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	623 call dlaset ('All', ncv, ncv,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	624 & zero , one, workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	625 & ldq)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	626 call dlahqr (.true., .true. , ncv,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	627 & 1 , ncv , workl(iuptri),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	628 & ldh , workl(iheigr), workl(iheigi),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	629 & 1 , ncv , workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	630 & ldq , 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	631 call dcopy (ncv , workl(invsub+ncv-1), ldq,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	632 & workl(ihbds), 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	633 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	634 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	635 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	636 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	637 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	638 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	639 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	640 call dvout (logfil, ncv, workl(iheigr), 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	641 & '_neupd: Real part of the eigenvalues of 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	642 call dvout (logfil, ncv, workl(iheigi), 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	643 & '_neupd: Imaginary part of the Eigenvalues of 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	644 call dvout (logfil, ncv, workl(ihbds), 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	645 & '_neupd: Last row of the Schur vector 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	646 if (msglvl .gt. 3) 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	647 call dmout (logfil , ncv, ncv ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	648 & workl(iuptri), ldh, 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	649 & '_neupd: The upper quasi-triangular 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	650 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	651 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	652 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	653 if (reord) 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	654 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	655 c %-----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	656 c \| Reorder the computed upper quasi-triangular 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	657 c %-----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	658 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	659 call dtrsen ('None' , '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	660 & select , ncv ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	661 & workl(iuptri), ldh ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	662 & workl(invsub), ldq ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	663 & workl(iheigr), workl(iheigi),
12277 3d38b4916cbf avoid memory overrun in ARPACK David Bateman <dbateman@free.fr> parents: 12274 diff changeset	664 & nconv2 , conds ,
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	665 & sep , workl(ihbds) ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	666 & ncv , iwork ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	667 & 1 , 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	668 c
12277 3d38b4916cbf avoid memory overrun in ARPACK David Bateman <dbateman@free.fr> parents: 12274 diff changeset	669 if (nconv2 .lt. nconv) then
3d38b4916cbf avoid memory overrun in ARPACK David Bateman <dbateman@free.fr> parents: 12274 diff changeset	670 nconv = nconv2
3d38b4916cbf avoid memory overrun in ARPACK David Bateman <dbateman@free.fr> parents: 12274 diff changeset	671 end if
3d38b4916cbf avoid memory overrun in ARPACK David Bateman <dbateman@free.fr> parents: 12274 diff changeset	672
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	673 if (ierr .eq. 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	674 info = 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	675 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	676 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	677 c
12277 3d38b4916cbf avoid memory overrun in ARPACK David Bateman <dbateman@free.fr> parents: 12274 diff changeset	678
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	679 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	680 call dvout (logfil, ncv, workl(iheigr), 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	681 & '_neupd: Real part of the eigenvalues of H--reordered')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	682 call dvout (logfil, ncv, workl(iheigi), 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	683 & '_neupd: Imag part of the eigenvalues of H--reordered')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	684 if (msglvl .gt. 3) 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	685 call dmout (logfil , ncv, ncv ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	686 & workl(iuptri), ldq, 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	687 & '_neupd: Quasi-triangular matrix after re-ordering')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	688 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	689 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	690 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	691 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	692 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	693 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	694 c \| Copy the last row of the Schur 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	695 c \| into workl(ihbds). This will be 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	696 c \| to compute the Ritz estimates 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	697 c \| 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	698 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	699 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	700 call dcopy (ncv, workl(invsub+ncv-1), ldq, workl(ihbds), 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	701 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	702 c %----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	703 c \| Place the computed eigenvalues of H into DR and DI \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	704 c \| if a spectral transformation was not 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	705 c %----------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	706 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	707 if (type .eq. 'REGULR') 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	708 call dcopy (nconv, workl(iheigr), 1, dr, 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	709 call dcopy (nconv, workl(iheigi), 1, di, 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	710 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	711 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	712 c %----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	713 c \| Compute the QR factorization of the matrix representing \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	714 c \| the wanted invariant subspace located in the first 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	715 c \| columns of workl(invsub,ldq). \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	716 c %----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	717 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	718 call dgeqr2 (ncv, nconv , workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	719 & ldq, workev, workev(ncv+1),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	720 & 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	721 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	722 c %---------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	723 c \| * Postmultiply V by Q using dorm2r . \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	724 c \| * Copy the first NCONV columns of VQ into Z. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	725 c \| * Postmultiply Z by R. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	726 c \| The N by NCONV matrix Z is now a matrix representation \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	727 c \| of the approximate invariant subspace associated with \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	728 c \| the Ritz values in workl(iheigr) and workl(iheigi) \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	729 c \| The first NCONV columns of V are now approximate Schur \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	730 c \| vectors associated with the real upper quasi-triangular \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	731 c \| matrix of order NCONV in workl(iuptri) \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	732 c %---------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	733 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	734 call dorm2r ('Right', 'Notranspose', n ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	735 & ncv , nconv , workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	736 & ldq , workev , 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	737 & ldv , workd(n+1) , 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	738 call dlacpy ('All', n, nconv, v, ldv, z, ldz)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	739 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	740 do 20 j=1, 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	741 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	742 c %---------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	743 c \| Perform both a column and row scaling if 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	744 c \| diagonal element of workl(invsub,ldq) is negative \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	745 c \| I'm lazy and don't take advantage of the upper \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	746 c \| quasi-triangular form of workl(iuptri,ldq) \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	747 c \| Note that since Q is orthogonal, R is a diagonal \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	748 c \| matrix consisting of plus or minus ones \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	749 c %---------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	750 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	751 if (workl(invsub+(j-1)*ldq+j-1) .lt. zero) 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	752 call dscal (nconv, -one, workl(iuptri+j-1), ldq)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	753 call dscal (nconv, -one, workl(iuptri+(j-1)*ldq), 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	754 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	755 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	756 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	757 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	758 if (howmny .eq. 'A') 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	759 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	760 c %--------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	761 c \| Compute the NCONV wanted eigenvectors of T \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	762 c \| located in workl(iuptri,ldq). \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	763 c %--------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	764 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	765 do 30 j=1, ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	766 if (j .le. nconv) 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	767 select(j) = .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	768 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	769 select(j) = .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	770 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	771 30 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	772 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	773 call dtrevc ('Right', 'Select' , select ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	774 & ncv , workl(iuptri), ldq ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	775 & vl , 1 , workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	776 & ldq , ncv , outncv ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	777 & workev , 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	778 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	779 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	780 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	781 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	782 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	783 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	784 c %------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	785 c \| Scale the returning eigenvectors so that their \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	786 c \| Euclidean norms are all one. LAPACK 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	787 c \| dtrevc returns each eigenvector normalized so \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	788 c \| that the element of largest magnitude has \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	789 c \| magnitude 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	790 c %------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	791 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	792 iconj = 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	793 do 40 j=1, 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	794 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	795 if ( workl(iheigi+j-1) .eq. zero ) 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	796 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	797 c %----------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	798 c \| real eigenvalue case \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	799 c %----------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	800 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	801 temp = dnrm2 ( ncv, workl(invsub+(j-1)*ldq), 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	802 call dscal ( ncv, one / temp,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	803 & workl(invsub+(j-1)*ldq), 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	804 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	805 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	806 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	807 c %-------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	808 c \| Complex conjugate pair case. Note that \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	809 c \| since the real and imaginary part 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	810 c \| the eigenvector are stored in consecutive \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	811 c \| columns, we further normalize by 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	812 c \| square root of two. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	813 c %-------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	814 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	815 if (iconj .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	816 temp = dlapy2 (dnrm2 (ncv,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	817 & workl(invsub+(j-1)*ldq),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	818 & 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	819 & dnrm2 (ncv,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	820 & workl(invsub+j*ldq),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	821 & 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	822 call dscal (ncv, one/temp,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	823 & workl(invsub+(j-1)*ldq), 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	824 call dscal (ncv, one/temp,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	825 & workl(invsub+j*ldq), 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	826 iconj = 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	827 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	828 iconj = 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	829 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	830 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	831 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	832 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	833 40 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	834 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	835 call dgemv ('T', ncv, nconv, one, workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	836 & ldq, workl(ihbds), 1, zero, workev, 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	837 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	838 iconj = 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	839 do 45 j=1, 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	840 if (workl(iheigi+j-1) .ne. zero) 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	841 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	842 c %-------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	843 c \| Complex conjugate pair case. Note that \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	844 c \| since the real and imaginary part 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	845 c \| the eigenvector are stored in consecutive \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	846 c %-------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	847 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	848 if (iconj .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	849 workev(j) = dlapy2 (workev(j), workev(j+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	850 workev(j+1) = workev(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	851 iconj = 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	852 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	853 iconj = 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	854 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	855 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	856 45 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	857 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	858 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	859 call dcopy (ncv, workl(invsub+ncv-1), ldq,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	860 & workl(ihbds), 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	861 call dvout (logfil, ncv, workl(ihbds), 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	862 & '_neupd: Last row of the eigenvector matrix for T')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	863 if (msglvl .gt. 3) 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	864 call dmout (logfil, ncv, ncv, workl(invsub), ldq,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	865 & ndigit, '_neupd: The eigenvector matrix for T')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	866 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	867 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	868 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	869 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	870 c \| Copy Ritz estimates into workl(ihbds) \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	871 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	872 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	873 call dcopy (nconv, workev, 1, workl(ihbds), 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	874 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	875 c %---------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	876 c \| Compute the QR factorization of the eigenvector 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	877 c \| associated with leading portion of T in the first 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	878 c \| columns of workl(invsub,ldq). \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	879 c %---------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	880 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	881 call dgeqr2 (ncv, nconv , workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	882 & ldq, workev, workev(ncv+1),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	883 & 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	884 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	885 c %----------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	886 c \| * Postmultiply Z by Q. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	887 c \| * Postmultiply Z by R. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	888 c \| The N by NCONV matrix Z is now contains 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	889 c \| Ritz vectors associated with the 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	890 c \| in workl(iheigr) and workl(iheigi). \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	891 c %----------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	892 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	893 call dorm2r ('Right', 'Notranspose', n ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	894 & ncv , nconv , workl(invsub),
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	895 & ldq , workev , z ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	896 & ldz , workd(n+1) , 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	897 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	898 call dtrmm ('Right' , 'Upper' , 'No transpose',
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	899 & 'Non-unit', n , 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	900 & one , workl(invsub), ldq ,
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	901 & z , ldz)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	902 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	903 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	904 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	905 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	906 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	907 c %------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	908 c \| An approximate invariant subspace is not needed. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	909 c \| Place the Ritz values computed DNAUPD into DR and DI \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	910 c %------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	911 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	912 call dcopy (nconv, workl(ritzr), 1, dr, 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	913 call dcopy (nconv, workl(ritzi), 1, di, 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	914 call dcopy (nconv, workl(ritzr), 1, workl(iheigr), 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	915 call dcopy (nconv, workl(ritzi), 1, workl(iheigi), 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	916 call dcopy (nconv, workl(bounds), 1, workl(ihbds), 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	917 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	918 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	919 c %------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	920 c \| Transform the Ritz values and possibly 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	921 c \| and corresponding error bounds of OP to those \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	922 c \| of Ax = lambdaB*x. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	923 c %------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	924 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	925 if (type .eq. 'REGULR') 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	926 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	927 if (rvec)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	928 & call dscal (ncv, rnorm, workl(ihbds), 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	929 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	930 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	931 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	932 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	933 c \| A spectral transformation was 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	934 c \| * Determine the Ritz estimates of the \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	935 c \| Ritz values in the original system. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	936 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	937 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	938 if (type .eq. 'SHIFTI') 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	939 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	940 if (rvec)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	941 & call dscal (ncv, rnorm, workl(ihbds), 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	942 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	943 do 50 k=1, ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	944 temp = dlapy2 ( workl(iheigr+k-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	945 & workl(iheigi+k-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	946 workl(ihbds+k-1) = abs( workl(ihbds+k-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	947 & / temp / temp
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	948 50 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	949 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	950 else if (type .eq. 'REALPT') 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	951 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	952 do 60 k=1, ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	953 60 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	954 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	955 else if (type .eq. 'IMAGPT') 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	956 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	957 do 70 k=1, ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	958 70 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	959 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	960 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	961 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	962 c %-----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	963 c \| * Transform the Ritz values back to the original system. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	964 c \| For TYPE = 'SHIFTI' the transformation is \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	965 c \| lambda = 1/theta + sigma \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	966 c \| For TYPE = 'REALPT' or 'IMAGPT' the user must from \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	967 c \| Rayleigh quotients or a projection. See remark 3 above.\|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	968 c \| NOTES: \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	969 c \| *The Ritz vectors are not affected by the transformation. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	970 c %-----------------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	971 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	972 if (type .eq. 'SHIFTI') 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	973 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	974 do 80 k=1, ncv
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	975 temp = dlapy2 ( workl(iheigr+k-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	976 & workl(iheigi+k-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	977 workl(iheigr+k-1) = workl(iheigr+k-1)/temp/temp
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	978 & + sigmar
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	979 workl(iheigi+k-1) = -workl(iheigi+k-1)/temp/temp
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	980 & + sigmai
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	981 80 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	982 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	983 call dcopy (nconv, workl(iheigr), 1, dr, 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	984 call dcopy (nconv, workl(iheigi), 1, di, 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	985 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	986 else if (type .eq. 'REALPT' .or. type .eq. 'IMAGPT') 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	987 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	988 call dcopy (nconv, workl(iheigr), 1, dr, 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	989 call dcopy (nconv, workl(iheigi), 1, di, 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	990 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	991 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	992 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	993 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	994 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	995 if (type .eq. 'SHIFTI' .and. 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	996 call dvout (logfil, nconv, dr, 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	997 & '_neupd: Untransformed real part of the Ritz valuess.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	998 call dvout (logfil, nconv, di, 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	999 & '_neupd: Untransformed imag part of the Ritz valuess.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1000 call dvout (logfil, nconv, workl(ihbds), 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	1001 & '_neupd: Ritz estimates of untransformed 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	1002 else if (type .eq. 'REGULR' .and. 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	1003 call dvout (logfil, nconv, dr, 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	1004 & '_neupd: Real parts 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	1005 call dvout (logfil, nconv, di, 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	1006 & '_neupd: Imag parts 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	1007 call dvout (logfil, nconv, workl(ihbds), 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	1008 & '_neupd: Associated Ritz estimates.')
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1009 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	1010 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1011 c %-------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1012 c \| Eigenvector Purification step. Formally perform \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1013 c \| one of inverse subspace iteration. Only 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	1014 c \| for MODE = 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	1015 c %-------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1016 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1017 if (rvec .and. howmny .eq. 'A' .and. type .eq. 'SHIFTI') 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	1018 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1019 c %------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1020 c \| Purify the computed Ritz vectors by adding a \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1021 c \| little bit of the 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	1022 c \| T \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1023 c \| resid(:)*( e s ) / theta \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1024 c \| NCV \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1025 c \| where H s = s theta. Remember that when theta \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1026 c \| has nonzero imaginary part, 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	1027 c \| Ritz vector is stored across two columns of Z. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1028 c %------------------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1029 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1030 iconj = 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	1031 do 110 j=1, 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	1032 if (workl(iheigi+j-1) .eq. zero) 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	1033 workev(j) = workl(invsub+(j-1)*ldq+ncv-1) /
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1034 & workl(iheigr+j-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	1035 else if (iconj .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	1036 temp = dlapy2 ( workl(iheigr+j-1), workl(iheigi+j-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	1037 workev(j) = ( workl(invsub+(j-1)ldq+ncv-1)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1038 & workl(iheigr+j-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	1039 & workl(invsub+jldq+ncv-1)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1040 & workl(iheigi+j-1) ) / temp / temp
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1041 workev(j+1) = ( workl(invsub+jldq+ncv-1)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1042 & workl(iheigr+j-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	1043 & workl(invsub+(j-1)ldq+ncv-1)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1044 & workl(iheigi+j-1) ) / temp / temp
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1045 iconj = 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	1046 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	1047 iconj = 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	1048 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	1049 110 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	1050 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1051 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1052 c \| Perform a rank one update to Z 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	1053 c \| purify all the Ritz vectors together. \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1054 c %---------------------------------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1055 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1056 call dger (n, nconv, one, resid, 1, workev, 1, z, ldz)
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1057 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1058 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	1059 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1060 9000 continue
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1061 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1062 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	1063 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1064 c %---------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1065 c \| End of DNEUPD \|
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1066 c %---------------%
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1067 c
9f5d2ef078e8 import ARPACK sources to libcruft from Debian package libarpack2 2.1+parpack96.dfsg-3+b1 John W. Eaton <jwe@octave.org> parents: diff changeset	1068 end

Mercurial > octave-nkf

annotate libcruft/arpack/src/dneupd.f @ 12277:3d38b4916cbf release-3-4-x