← Back to branch summary

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

« back to all changes in this revision

Viewing changes to PARPACK/SRC/BLACS/pznaitr.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:
PARPACK/SRC/BLACS/pznaitr.f
137
137
c\Routines called:
138
138
c     pzgetv0   Parallel ARPACK routine to generate the initial vector.
139
139
c     pivout    Parallel ARPACK utility routine that prints integers.
140
 
c     second    ARPACK utility routine for timing.
 
140
c     arscnd    ARPACK utility routine for timing.
141
141
c     pzmout    Parallel ARPACK utility routine that prints matrices
142
142
c     pzvout    Parallel ARPACK utility routine that prints vectors.
143
143
c     zlanhs    LAPACK routine that computes various norms of a matrix.
296
296
c     %----------------------%
297
297
c
298
298
      external   zaxpy, zcopy, zscal, zgemv, pzgetv0, dlabad, 
299
 
     &           zdscal, pzvout, pzmout, pivout, second
 
299
     &           zdscal, pzvout, pzmout, pivout, arscnd
300
300
c
301
301
c     %--------------------%
302
302
c     | External Functions |
349
349
c        | & message level for debugging |
350
350
c        %-------------------------------%
351
351
c
352
 
         call second (t0)
 
352
         call arscnd (t0)
353
353
         msglvl = mcaitr
354
354
c 
355
355
c        %------------------------------%
460
460
c              %------------------------------------------------%
461
461
c
462
462
               info = j - 1
463
 
               call second (t1)
 
463
               call arscnd (t1)
464
464
               tcaitr = tcaitr + (t1 - t0)
465
465
               ido = 99
466
466
               go to 9000
500
500
c
501
501
         step3 = .true.
502
502
         nopx  = nopx + 1
503
 
         call second (t2)
 
503
         call arscnd (t2)
504
504
         call zcopy (n, v(1,j), 1, workd(ivj), 1)
505
505
         ipntr(1) = ivj
506
506
         ipntr(2) = irj
520
520
c        | if step3 = .true.                |
521
521
c        %----------------------------------%
522
522
c
523
 
         call second (t3)
 
523
         call arscnd (t3)
524
524
         tmvopx = tmvopx + (t3 - t2)
525
525
 
526
526
         step3 = .false.
536
536
c        |          factorization to length j.   |
537
537
c        %---------------------------------------%
538
538
c
539
 
         call second (t2)
 
539
         call arscnd (t2)
540
540
         if (bmat .eq. 'G') then
541
541
            nbx = nbx + 1
542
542
            step4 = .true.
561
561
c        %----------------------------------%
562
562
c
563
563
         if (bmat .eq. 'G') then
564
 
            call second (t3)
 
564
            call arscnd (t3)
565
565
            tmvbx = tmvbx + (t3 - t2)
566
566
         end if
567
567
c 
608
608
c
609
609
         if (j .gt. 1) h(j,j-1) = dcmplx(betaj, rzero)
610
610
c
611
 
         call second (t4)
 
611
         call arscnd (t4)
612
612
c 
613
613
         orth1 = .true.
614
614
c 
615
 
         call second (t2)
 
615
         call arscnd (t2)
616
616
         if (bmat .eq. 'G') then
617
617
            nbx = nbx + 1
618
618
            call zcopy (n, resid, 1, workd(irj), 1)
636
636
c        %---------------------------------------------------%
637
637
c
638
638
         if (bmat .eq. 'G') then
639
 
            call second (t3)
 
639
            call arscnd (t3)
640
640
            tmvbx = tmvbx + (t3 - t2)
641
641
         end if
642
642
c 
716
716
         call zaxpy (j, one, workl(1), 1, h(1,j), 1)
717
717
c 
718
718
         orth2 = .true.
719
 
         call second (t2)
 
719
         call arscnd (t2)
720
720
         if (bmat .eq. 'G') then
721
721
            nbx = nbx + 1
722
722
            call zcopy (n, resid, 1, workd(irj), 1)
740
740
c        %---------------------------------------------------%
741
741
c
742
742
         if (bmat .eq. 'G') then
743
 
            call second (t3)
 
743
            call arscnd (t3)
744
744
            tmvbx = tmvbx + (t3 - t2)
745
745
         end if 
746
746
c
819
819
         rstart = .false.
820
820
         orth2  = .false.
821
821
c 
822
 
         call second (t5)
 
822
         call arscnd (t5)
823
823
         titref = titref + (t5 - t4)
824
824
c 
825
825
c        %------------------------------------%
828
828
c
829
829
         j = j + 1
830
830
         if (j .gt. k+np) then
831
 
            call second (t1)
 
831
            call arscnd (t1)
832
832
            tcaitr = tcaitr + (t1 - t0)
833
833
            ido = 99
834
834
            do 110 i = max(1,k), k+np-1