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

« back to all changes in this revision

Viewing changes to PARPACK/SRC/BLACS/pssaup2.f

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2013-11-13 09:56:18 UTC
  • mfrom: (1.3.3)
  • Revision ID: package-import@ubuntu.com-20131113095618-xoq2ckxy41bk7ilb
Tags: 3.1.4-1
New upstream release

Show diffs side-by-side

added added

removed removed

Lines of Context:
153
153
c              last component of the eigenvectors of a symmetric
154
154
c              tridiagonal matrix using the implicit QL or QR method.
155
155
c     pivout   Parallel ARPACK utility routine that prints integers.
156
 
c     second   ARPACK utility routine for timing.
 
156
c     arscnd   ARPACK utility routine for timing.
157
157
c     psvout   Parallel ARPACK utility routine that prints vectors.
158
158
c     pslamch  ScaLAPACK routine that determines machine constants.
159
159
c     scopy    Level 1 BLAS that copies one vector to another.
253
253
c
254
254
      external   scopy, psgetv0, pssaitr, sscal, ssconv, 
255
255
     &           psseigt, pssgets, pssapps,
256
 
     &           ssortr, psvout, pivout, second
 
256
     &           ssortr, psvout, pivout, arscnd
257
257
c
258
258
c     %--------------------%
259
259
c     | External Functions |
280
280
c        | & message level for debugging |
281
281
c        %-------------------------------%
282
282
c
283
 
         call second (t0)
 
283
         call arscnd (t0)
284
284
         msglvl = msaup2
285
285
c
286
286
c        %---------------------------------%
791
791
c        %---------------------------------------------%
792
792
c
793
793
         cnorm = .true.
794
 
         call second (t2)
 
794
         call arscnd (t2)
795
795
         if (bmat .eq. 'G') then
796
796
            nbx = nbx + 1
797
797
            call scopy (n, resid, 1, workd(n+1), 1)
816
816
c        %----------------------------------%
817
817
c
818
818
         if (bmat .eq. 'G') then
819
 
            call second (t3)
 
819
            call arscnd (t3)
820
820
            tmvbx = tmvbx + (t3 - t2)
821
821
         end if
822
822
859
859
c     | Error exit |
860
860
c     %------------%
861
861
c
862
 
      call second (t1)
 
862
      call arscnd (t1)
863
863
      tsaup2 = t1 - t0
864
864
865
865
 9000 continue