← Back to branch summary

~ubuntu-branches/debian/sid/arpack/sid

« back to all changes in this revision

Viewing changes to PARPACK/SRC/MPI/pdneigh.f

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2015-10-12 09:33:05 UTC
  • mfrom: (1.3.5)
  • Revision ID: package-import@ubuntu.com-20151012093305-w41sv3g9nzbaodfz
Tags: 3.3.0-1
* New upstream release
* Update of the homepage (moved to github)
* libparpack links against libarpack (instead of static build)
* Update the symbol files. No need to change of SONAME:
  - About the parpack change: https://github.com/opencollab/arpack-ng/pull/21
  - Transfer the call to the library:
    https://github.com/opencollab/arpack-ng/pull/19
* Standards-Version updated to version 3.9.6

Show diffs side-by-side

added added

removed removed

Lines of Context:
PARPACK/SRC/MPI/pdneigh.f
65
65
c     xxxxxx  real
66
66
c
67
67
c\Routines called:
68
 
c     dlaqrb  ARPACK routine to compute the real Schur form of an
 
68
c     dlahqr  ARPACK routine to compute the real Schur form of an
69
69
c             upper Hessenberg matrix and last row of the Schur vectors.
70
70
c     arscnd  ARPACK utility routine for timing.
71
71
c     dmout   ARPACK utility routine that prints matrices
157
157
c     | External Subroutines |
158
158
c     %----------------------%
159
159
c
160
 
      external   dcopy, dlacpy, dlaqrb, dtrevc, pdvout, arscnd
 
160
      external   dcopy, dlacpy, dlahqr, dtrevc, dvout, arscnd
161
161
c
162
162
c     %--------------------%
163
163
c     | External Functions |
195
195
c     | 1. Compute the eigenvalues, the last components of the    |
196
196
c     |    corresponding Schur vectors and the full Schur form T  |
197
197
c     |    of the current upper Hessenberg matrix H.              |
198
 
c     | dlaqrb returns the full Schur form of H in WORKL(1:N**2)  |
 
198
c     | dlahqr returns the full Schur form of H in WORKL(1:N**2)  |
199
199
c     | and the last components of the Schur vectors in BOUNDS.   |
200
200
c     %-----------------------------------------------------------%
201
201
c
202
202
      call dlacpy ('All', n, n, h, ldh, workl, n)
203
 
      call dlaqrb (.true., n, 1, n, workl, n, ritzr, ritzi, bounds,
204
 
     &             ierr)
 
203
      do 5 j = 1, n-1
 
204
          bounds(j) = zero
 
205
   5  continue
 
206
      bounds(n) = 1
 
207
      call dlahqr(.true., .true., n, 1, n, workl, n, ritzr, ritzi, 1, 1,
 
208
     &            bounds, 1, ierr)
205
209
      if (ierr .ne. 0) go to 9000
206
210
c
207
211
      if (msglvl .gt. 1) then