← Back to branch summary

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

« back to all changes in this revision

Viewing changes to PARPACK/SRC/BLACS/pdnaitr.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/pdnaitr.f
138
138
c\Routines called:
139
139
c     pdgetv0  Parallel ARPACK routine to generate the initial vector.
140
140
c     pivout   Parallel ARPACK utility routine that prints integers.
141
 
c     second   ARPACK utility routine for timing.
 
141
c     arscnd   ARPACK utility routine for timing.
142
142
c     pdmout   Parallel ARPACK utility routine that prints matrices
143
143
c     pdvout   Parallel ARPACK utility routine that prints vectors.
144
144
c     dlabad   LAPACK routine that computes machine constants.
287
287
c     %----------------------%
288
288
c
289
289
      external   daxpy, dcopy, dscal, dgemv, pdgetv0, dlabad, 
290
 
     &           pdvout, pdmout, pivout, second
 
290
     &           pdvout, pdmout, pivout, arscnd
291
291
c
292
292
c     %--------------------%
293
293
c     | External Functions |
338
338
c        | & message level for debugging |
339
339
c        %-------------------------------%
340
340
c
341
 
         call second (t0)
 
341
         call arscnd (t0)
342
342
         msglvl = mnaitr
343
343
c 
344
344
c        %------------------------------%
449
449
c              %------------------------------------------------%
450
450
c
451
451
               info = j - 1
452
 
               call second (t1)
 
452
               call arscnd (t1)
453
453
               tnaitr = tnaitr + (t1 - t0)
454
454
               ido = 99
455
455
               go to 9000
489
489
c
490
490
         step3 = .true.
491
491
         nopx  = nopx + 1
492
 
         call second (t2)
 
492
         call arscnd (t2)
493
493
         call dcopy (n, v(1,j), 1, workd(ivj), 1)
494
494
         ipntr(1) = ivj
495
495
         ipntr(2) = irj
509
509
c        | if step3 = .true.                |
510
510
c        %----------------------------------%
511
511
c
512
 
         call second (t3)
 
512
         call arscnd (t3)
513
513
         tmvopx = tmvopx + (t3 - t2)
514
514
 
515
515
         step3 = .false.
525
525
c        |          factorization to length j.   |
526
526
c        %---------------------------------------%
527
527
c
528
 
         call second (t2)
 
528
         call arscnd (t2)
529
529
         if (bmat .eq. 'G') then
530
530
            nbx = nbx + 1
531
531
            step4 = .true.
550
550
c        %----------------------------------%
551
551
c
552
552
         if (bmat .eq. 'G') then
553
 
            call second (t3)
 
553
            call arscnd (t3)
554
554
            tmvbx = tmvbx + (t3 - t2)
555
555
         end if
556
556
c 
597
597
c
598
598
         if (j .gt. 1) h(j,j-1) = betaj
599
599
c
600
 
         call second (t4)
 
600
         call arscnd (t4)
601
601
c 
602
602
         orth1 = .true.
603
603
c
604
 
         call second (t2)
 
604
         call arscnd (t2)
605
605
         if (bmat .eq. 'G') then
606
606
            nbx = nbx + 1
607
607
            call dcopy (n, resid, 1, workd(irj), 1)
625
625
c        %---------------------------------------------------%
626
626
c
627
627
         if (bmat .eq. 'G') then
628
 
            call second (t3)
 
628
            call arscnd (t3)
629
629
            tmvbx = tmvbx + (t3 - t2)
630
630
         end if
631
631
c 
704
704
         call daxpy (j, one, workl(1), 1, h(1,j), 1)
705
705
c 
706
706
         orth2 = .true.
707
 
         call second (t2)
 
707
         call arscnd (t2)
708
708
         if (bmat .eq. 'G') then
709
709
            nbx = nbx + 1
710
710
            call dcopy (n, resid, 1, workd(irj), 1)
728
728
c        %---------------------------------------------------%
729
729
c
730
730
         if (bmat .eq. 'G') then
731
 
            call second (t3)
 
731
            call arscnd (t3)
732
732
            tmvbx = tmvbx + (t3 - t2)
733
733
         end if
734
734
c
807
807
         rstart = .false.
808
808
         orth2  = .false.
809
809
c 
810
 
         call second (t5)
 
810
         call arscnd (t5)
811
811
         titref = titref + (t5 - t4)
812
812
c 
813
813
c        %------------------------------------%
816
816
c
817
817
         j = j + 1
818
818
         if (j .gt. k+np) then
819
 
            call second (t1)
 
819
            call arscnd (t1)
820
820
            tnaitr = tnaitr + (t1 - t0)
821
821
            ido = 99
822
822
            do 110 i = max(1,k), k+np-1