← Back to branch summary

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

« back to all changes in this revision

Viewing changes to PARPACK/SRC/BLACS/pdsaitr.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/pdsaitr.f
280
280
c     %----------------------%
281
281
c
282
282
      external   daxpy, dcopy, dscal, dgemv, pdgetv0, pdvout, pdmout,
283
 
     &           dlascl, pivout, second
 
283
     &           dlascl, pivout, arscnd
284
284
c
285
285
c     %--------------------%
286
286
c     | External Functions |
325
325
c        | & message level for debugging |
326
326
c        %-------------------------------%
327
327
c
328
 
         call second (t0)
 
328
         call arscnd (t0)
329
329
         msglvl = msaitr
330
330
c 
331
331
c        %------------------------------%
445
445
c              %------------------------------------------------%
446
446
c
447
447
               info = j - 1
448
 
               call second (t1)
 
448
               call arscnd (t1)
449
449
               tsaitr = tsaitr + (t1 - t0)
450
450
               ido = 99
451
451
               go to 9000
485
485
c
486
486
         step3 = .true.
487
487
         nopx  = nopx + 1
488
 
         call second (t2)
 
488
         call arscnd (t2)
489
489
         call dcopy (n, v(1,j), 1, workd(ivj), 1)
490
490
         ipntr(1) = ivj
491
491
         ipntr(2) = irj
504
504
c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j}.   |
505
505
c        %-----------------------------------%
506
506
c
507
 
         call second (t3)
 
507
         call arscnd (t3)
508
508
         tmvopx = tmvopx + (t3 - t2)
509
509
c 
510
510
         step3 = .false.
525
525
c        %-------------------------------------------%
526
526
c
527
527
         if (mode .eq. 2) go to 65
528
 
         call second (t2)
 
528
         call arscnd (t2)
529
529
         if (bmat .eq. 'G') then
530
530
            nbx = nbx + 1
531
531
            step4 = .true.
549
549
c        %-----------------------------------%
550
550
c
551
551
         if (bmat .eq. 'G') then
552
 
            call second (t3)
 
552
            call arscnd (t3)
553
553
            tmvbx = tmvbx + (t3 - t2)
554
554
         end if 
555
555
c
621
621
         else
622
622
            h(j,1) = rnorm
623
623
         end if
624
 
         call second (t4)
 
624
         call arscnd (t4)
625
625
c 
626
626
         orth1 = .true.
627
627
         iter  = 0
628
628
c 
629
 
         call second (t2)
 
629
         call arscnd (t2)
630
630
         if (bmat .eq. 'G') then
631
631
            nbx = nbx + 1
632
632
            call dcopy (n, resid, 1, workd(irj), 1)
650
650
c        %---------------------------------------------------%
651
651
c
652
652
         if (bmat .eq. 'G') then
653
 
            call second (t3)
 
653
            call arscnd (t3)
654
654
            tmvbx = tmvbx + (t3 - t2)
655
655
         end if
656
656
c 
726
726
         h(j,2) = h(j,2) + workl(j)
727
727
c 
728
728
         orth2 = .true.
729
 
         call second (t2)
 
729
         call arscnd (t2)
730
730
         if (bmat .eq. 'G') then
731
731
            nbx = nbx + 1
732
732
            call dcopy (n, resid, 1, workd(irj), 1)
750
750
c        %---------------------------------------------------%
751
751
c
752
752
         if (bmat .eq. 'G') then
753
 
            call second (t3)
 
753
            call arscnd (t3)
754
754
            tmvbx = tmvbx + (t3 - t2)
755
755
         end if
756
756
c
823
823
         rstart = .false.
824
824
         orth2  = .false.
825
825
c 
826
 
         call second (t5)
 
826
         call arscnd (t5)
827
827
         titref = titref + (t5 - t4)
828
828
c 
829
829
c        %----------------------------------------------------------%
847
847
c
848
848
         j = j + 1
849
849
         if (j .gt. k+np) then
850
 
            call second (t1)
 
850
            call arscnd (t1)
851
851
            tsaitr = tsaitr + (t1 - t0)
852
852
            ido = 99
853
853
c