1 diff -urN original/ARPACK/SRC/cgetv0.f ARPACK/SRC/cgetv0.f
2 --- original/ARPACK/SRC/cgetv0.f 1999-04-08 00:46:44.000000000 +0200
3 +++ ARPACK/SRC/cgetv0.f 2005-06-30 23:20:38.850796569 +0200
8 -c second ARPACK utility routine for timing.
9 +c arscnd ARPACK utility routine for timing.
10 c cvout ARPACK utility routine that prints vectors.
11 c clarnv LAPACK routine for generating a random vector.
12 c cgemv Level 2 BLAS routine for matrix vector multiplication.
14 c | External Subroutines |
15 c %----------------------%
17 - external ccopy, cgemv, clarnv, cvout, second
18 + external ccopy, cgemv, clarnv, cvout, arscnd
20 c %--------------------%
21 c | External Functions |
23 c | & message level for debugging |
24 c %-------------------------------%
32 c | the generalized problem when B is possibly (singular). |
33 c %----------------------------------------------------------%
37 if (bmat .eq. 'G') then
46 tmvopx = tmvopx + (t3 - t2)
48 c %------------------------------------------------------%
50 c | Compute B-norm of starting vector. |
51 c %------------------------------------------------------%
56 if (bmat .eq. 'G') then
61 if (bmat .eq. 'G') then
64 tmvbx = tmvbx + (t3 - t2)
68 c | Compute the B-norm of the orthogonalized starting vector |
69 c %----------------------------------------------------------%
73 if (bmat .eq. 'G') then
75 call ccopy (n, resid, 1, workd(n+1), 1)
79 if (bmat .eq. 'G') then
82 tmvbx = tmvbx + (t3 - t2)
91 tgetv0 = tgetv0 + (t1 - t0)
94 diff -urN original/ARPACK/SRC/cnaitr.f ARPACK/SRC/cnaitr.f
95 --- original/ARPACK/SRC/cnaitr.f 1996-08-27 18:12:56.000000000 +0200
96 +++ ARPACK/SRC/cnaitr.f 2005-06-30 23:20:38.856795869 +0200
99 c cgetv0 ARPACK routine to generate the initial vector.
100 c ivout ARPACK utility routine that prints integers.
101 -c second ARPACK utility routine for timing.
102 +c arscnd ARPACK utility routine for timing.
103 c cmout ARPACK utility routine that prints matrices
104 c cvout ARPACK utility routine that prints vectors.
105 c clanhs LAPACK routine that computes various norms of a matrix.
107 c %----------------------%
109 external caxpy, ccopy, cscal, csscal, cgemv, cgetv0,
110 - & slabad, cvout, cmout, ivout, second
111 + & slabad, cvout, cmout, ivout, arscnd
113 c %--------------------%
114 c | External Functions |
116 c | & message level for debugging |
117 c %-------------------------------%
123 c %------------------------------%
125 c %------------------------------------------------%
130 tcaitr = tcaitr + (t1 - t0)
139 call ccopy (n, v(1,j), 1, workd(ivj), 1)
143 c | if step3 = .true. |
144 c %----------------------------------%
148 tmvopx = tmvopx + (t3 - t2)
152 c | factorization to length j. |
153 c %---------------------------------------%
157 if (bmat .eq. 'G') then
161 c %----------------------------------%
163 if (bmat .eq. 'G') then
166 tmvbx = tmvbx + (t3 - t2)
169 @@ -583,11 +583,11 @@
171 if (j .gt. 1) h(j,j-1) = cmplx(betaj, rzero)
180 if (bmat .eq. 'G') then
182 call ccopy (n, resid, 1, workd(irj), 1)
184 c %---------------------------------------------------%
186 if (bmat .eq. 'G') then
189 tmvbx = tmvbx + (t3 - t2)
193 call caxpy (j, one, workd(irj), 1, h(1,j), 1)
198 if (bmat .eq. 'G') then
200 call ccopy (n, resid, 1, workd(irj), 1)
202 c %---------------------------------------------------%
204 if (bmat .eq. 'G') then
207 tmvbx = tmvbx + (t3 - t2)
216 titref = titref + (t5 - t4)
218 c %------------------------------------%
222 if (j .gt. k+np) then
225 tcaitr = tcaitr + (t1 - t0)
227 do 110 i = max(1,k), k+np-1
228 diff -urN original/ARPACK/SRC/cnapps.f ARPACK/SRC/cnapps.f
229 --- original/ARPACK/SRC/cnapps.f 1998-05-20 16:45:03.000000000 +0200
230 +++ ARPACK/SRC/cnapps.f 2005-06-30 23:20:38.859795518 +0200
234 c ivout ARPACK utility routine that prints integers.
235 -c second ARPACK utility routine for timing.
236 +c arscnd ARPACK utility routine for timing.
237 c cmout ARPACK utility routine that prints matrices
238 c cvout ARPACK utility routine that prints vectors.
239 c clacpy LAPACK matrix copy routine.
241 c %----------------------%
243 external caxpy, ccopy, cgemv, cscal, clacpy, clartg,
244 - & cvout, claset, slabad, cmout, second, ivout
245 + & cvout, claset, slabad, cmout, arscnd, ivout
247 c %--------------------%
248 c | External Functions |
250 c | & message level for debugging |
251 c %-------------------------------%
264 tcapps = tcapps + (t1 - t0)
267 diff -urN original/ARPACK/SRC/cnaup2.f ARPACK/SRC/cnaup2.f
268 --- original/ARPACK/SRC/cnaup2.f 2000-06-01 22:28:53.000000000 +0200
269 +++ ARPACK/SRC/cnaup2.f 2005-06-30 23:20:38.864794935 +0200
271 c cngets ARPACK reorder Ritz values and error bounds routine.
272 c csortc ARPACK sorting routine.
273 c ivout ARPACK utility routine that prints integers.
274 -c second ARPACK utility routine for timing.
275 +c arscnd ARPACK utility routine for timing.
276 c cmout ARPACK utility routine that prints matrices
277 c cvout ARPACK utility routine that prints vectors.
278 c svout ARPACK utility routine that prints vectors.
280 c %----------------------%
282 external ccopy, cgetv0, cnaitr, cneigh, cngets, cnapps,
283 - & csortc, cswap, cmout, cvout, ivout, second
284 + & csortc, cswap, cmout, cvout, ivout, arscnd
286 c %--------------------%
287 c | External functions |
298 c %---------------------------------------------%
303 if (bmat .eq. 'G') then
305 call ccopy (n, resid, 1, workd(n+1), 1)
307 c %----------------------------------%
309 if (bmat .eq. 'G') then
312 tmvbx = tmvbx + (t3 - t2)
324 diff -urN original/ARPACK/SRC/cnaupd.f ARPACK/SRC/cnaupd.f
325 --- original/ARPACK/SRC/cnaupd.f 2002-08-15 07:50:57.000000000 +0200
326 +++ ARPACK/SRC/cnaupd.f 2005-06-30 23:20:38.869794351 +0200
328 c cstatn ARPACK routine that initializes the timing variables.
329 c ivout ARPACK utility routine that prints integers.
330 c cvout ARPACK utility routine that prints vectors.
331 -c second ARPACK utility routine for timing.
332 +c arscnd ARPACK utility routine for timing.
333 c slamch LAPACK routine that determines machine constants.
337 c | External Subroutines |
338 c %----------------------%
340 - external cnaup2, cvout, ivout, second, cstatn
341 + external cnaup2, cvout, ivout, arscnd, cstatn
343 c %--------------------%
344 c | External Functions |
346 c %-------------------------------%
355 & '_naupd: Associated Ritz estimates')
362 if (msglvl .gt. 0) then
363 diff -urN original/ARPACK/SRC/cneigh.f ARPACK/SRC/cneigh.f
364 --- original/ARPACK/SRC/cneigh.f 1996-04-20 19:15:53.000000000 +0200
365 +++ ARPACK/SRC/cneigh.f 2005-06-30 23:20:38.870794234 +0200
369 c ivout ARPACK utility routine that prints integers.
370 -c second ARPACK utility routine for timing.
371 +c arscnd ARPACK utility routine for timing.
372 c cmout ARPACK utility routine that prints matrices
373 c cvout ARPACK utility routine that prints vectors.
374 c svout ARPACK utility routine that prints vectors.
376 c %----------------------%
378 external clacpy, clahqr, ctrevc, ccopy,
379 - & csscal, cmout, cvout, second
380 + & csscal, cmout, cvout, arscnd
382 c %--------------------%
383 c | External Functions |
385 c | & message level for debugging |
386 c %-------------------------------%
392 if (msglvl .gt. 2) then
394 & '_neigh: Ritz estimates for the eigenvalues of H')
399 tceigh = tceigh + (t1 - t0)
402 diff -urN original/ARPACK/SRC/cngets.f ARPACK/SRC/cngets.f
403 --- original/ARPACK/SRC/cngets.f 1996-04-20 19:15:56.000000000 +0200
404 +++ ARPACK/SRC/cngets.f 2005-06-30 23:20:38.877793417 +0200
407 c csortc ARPACK sorting routine.
408 c ivout ARPACK utility routine that prints integers.
409 -c second ARPACK utility routine for timing.
410 +c arscnd ARPACK utility routine for timing.
411 c cvout ARPACK utility routine that prints vectors.
415 c | External Subroutines |
416 c %----------------------%
418 - external cvout, csortc, second
419 + external cvout, csortc, arscnd
421 c %-----------------------%
422 c | Executable Statements |
424 c | & message level for debugging |
425 c %-------------------------------%
431 call csortc (which, .true., kev+np, ritz, bounds)
438 tcgets = tcgets + (t1 - t0)
440 if (msglvl .gt. 0) then
441 diff -urN original/ARPACK/SRC/dgetv0.f ARPACK/SRC/dgetv0.f
442 --- original/ARPACK/SRC/dgetv0.f 1999-04-08 00:42:47.000000000 +0200
443 +++ ARPACK/SRC/dgetv0.f 2005-06-30 23:20:38.882792833 +0200
445 c TR95-13, Department of Computational and Applied Mathematics.
448 -c second ARPACK utility routine for timing.
449 +c arscnd ARPACK utility routine for timing.
450 c dvout ARPACK utility routine for vector output.
451 c dlarnv LAPACK routine for generating a random vector.
452 c dgemv Level 2 BLAS routine for matrix vector multiplication.
454 c | External Subroutines |
455 c %----------------------%
457 - external dlarnv, dvout, dcopy, dgemv, second
458 + external dlarnv, dvout, dcopy, dgemv, arscnd
460 c %--------------------%
461 c | External Functions |
463 c | & message level for debugging |
464 c %-------------------------------%
472 c | the generalized problem when B is possibly (singular). |
473 c %----------------------------------------------------------%
477 if (bmat .eq. 'G') then
483 if (bmat .eq. 'G') then
486 tmvopx = tmvopx + (t3 - t2)
490 c | Compute B-norm of starting vector. |
491 c %------------------------------------------------------%
496 if (bmat .eq. 'G') then
501 if (bmat .eq. 'G') then
504 tmvbx = tmvbx + (t3 - t2)
508 c | Compute the B-norm of the orthogonalized starting vector |
509 c %----------------------------------------------------------%
513 if (bmat .eq. 'G') then
515 call dcopy (n, resid, 1, workd(n+1), 1)
519 if (bmat .eq. 'G') then
522 tmvbx = tmvbx + (t3 - t2)
531 tgetv0 = tgetv0 + (t1 - t0)
534 diff -urN original/ARPACK/SRC/dnaitr.f ARPACK/SRC/dnaitr.f
535 --- original/ARPACK/SRC/dnaitr.f 1996-08-27 18:09:00.000000000 +0200
536 +++ ARPACK/SRC/dnaitr.f 2005-06-30 23:20:38.891791783 +0200
539 c dgetv0 ARPACK routine to generate the initial vector.
540 c ivout ARPACK utility routine that prints integers.
541 -c second ARPACK utility routine for timing.
542 +c arscnd ARPACK utility routine for timing.
543 c dmout ARPACK utility routine that prints matrices
544 c dvout ARPACK utility routine that prints vectors.
545 c dlabad LAPACK routine that computes machine constants.
547 c %----------------------%
549 external daxpy, dcopy, dscal, dgemv, dgetv0, dlabad,
550 - & dvout, dmout, ivout, second
551 + & dvout, dmout, ivout, arscnd
553 c %--------------------%
554 c | External Functions |
556 c | & message level for debugging |
557 c %-------------------------------%
563 c %------------------------------%
565 c %------------------------------------------------%
570 tnaitr = tnaitr + (t1 - t0)
579 call dcopy (n, v(1,j), 1, workd(ivj), 1)
583 c | if step3 = .true. |
584 c %----------------------------------%
588 tmvopx = tmvopx + (t3 - t2)
592 c | factorization to length j. |
593 c %---------------------------------------%
597 if (bmat .eq. 'G') then
601 c %----------------------------------%
603 if (bmat .eq. 'G') then
606 tmvbx = tmvbx + (t3 - t2)
609 @@ -576,11 +576,11 @@
611 if (j .gt. 1) h(j,j-1) = betaj
620 if (bmat .eq. 'G') then
622 call dcopy (n, resid, 1, workd(irj), 1)
624 c %---------------------------------------------------%
626 if (bmat .eq. 'G') then
629 tmvbx = tmvbx + (t3 - t2)
633 call daxpy (j, one, workd(irj), 1, h(1,j), 1)
638 if (bmat .eq. 'G') then
640 call dcopy (n, resid, 1, workd(irj), 1)
642 c %---------------------------------------------------%
644 if (bmat .eq. 'G') then
647 tmvbx = tmvbx + (t3 - t2)
656 titref = titref + (t5 - t4)
658 c %------------------------------------%
662 if (j .gt. k+np) then
665 tnaitr = tnaitr + (t1 - t0)
667 do 110 i = max(1,k), k+np-1
668 diff -urN original/ARPACK/SRC/dnapps.f ARPACK/SRC/dnapps.f
669 --- original/ARPACK/SRC/dnapps.f 1998-05-20 16:58:56.000000000 +0200
670 +++ ARPACK/SRC/dnapps.f 2005-06-30 23:20:38.896791199 +0200
674 c ivout ARPACK utility routine that prints integers.
675 -c second ARPACK utility routine for timing.
676 +c arscnd ARPACK utility routine for timing.
677 c dmout ARPACK utility routine that prints matrices.
678 c dvout ARPACK utility routine that prints vectors.
679 c dlabad LAPACK routine that computes machine constants.
681 c %----------------------%
683 external daxpy, dcopy, dscal, dlacpy, dlarfg, dlarf,
684 - & dlaset, dlabad, second, dlartg
685 + & dlaset, dlabad, arscnd, dlartg
687 c %--------------------%
688 c | External Functions |
690 c | & message level for debugging |
691 c %-------------------------------%
704 tnapps = tnapps + (t1 - t0)
707 --- ARPACK/SRC/dnaup2.f.orig 1996-08-16 22:28:09.000000000 +0200
708 +++ ARPACK/SRC/dnaup2.f 2022-12-15 21:12:42.177193805 +0100
709 @@ -143,7 +143,7 @@ c dneigh ARPACK compute Ritz values
710 c dngets ARPACK reorder Ritz values and error bounds routine.
711 c dsortc ARPACK sorting routine.
712 c ivout ARPACK utility routine that prints integers.
713 -c second ARPACK utility routine for timing.
714 +c arscnd ARPACK utility routine for timing.
715 c dmout ARPACK utility routine that prints matrices
716 c dvout ARPACK utility routine that prints vectors.
717 c dlamch LAPACK routine that determines machine constants.
718 @@ -234,7 +234,7 @@ c | External Subroutines |
719 c %----------------------%
721 external dcopy, dgetv0, dnaitr, dnconv, dneigh, dngets, dnapps,
722 - & dvout, ivout, second
723 + & dvout, ivout, arscnd
725 c %--------------------%
726 c | External Functions |
727 @@ -256,7 +256,7 @@ c %-----------------------%
736 @@ -755,7 +755,7 @@ c | the first step of the next ca
737 c %---------------------------------------------%
742 if (bmat .eq. 'G') then
744 call dcopy (n, resid, 1, workd(n+1), 1)
745 @@ -780,7 +780,7 @@ c | WORKD(1:N) := B*RESID
746 c %----------------------------------%
748 if (bmat .eq. 'G') then
751 tmvbx = tmvbx + (t3 - t2)
754 @@ -819,7 +819,7 @@ c %------------%
763 diff -urN original/ARPACK/SRC/dnaupd.f ARPACK/SRC/dnaupd.f
764 --- original/ARPACK/SRC/dnaupd.f 2002-08-23 07:29:34.000000000 +0200
765 +++ ARPACK/SRC/dnaupd.f 2005-06-30 23:20:38.906790031 +0200
767 c dnaup2 ARPACK routine that implements the Implicitly Restarted
769 c ivout ARPACK utility routine that prints integers.
770 -c second ARPACK utility routine for timing.
771 +c arscnd ARPACK utility routine for timing.
772 c dvout ARPACK utility routine that prints vectors.
773 c dlamch LAPACK routine that determines machine constants.
776 c | External Subroutines |
777 c %----------------------%
779 - external dnaup2, dvout, ivout, second, dstatn
780 + external dnaup2, dvout, ivout, arscnd, dstatn
782 c %--------------------%
783 c | External Functions |
785 c %-------------------------------%
794 & '_naupd: Associated Ritz estimates')
801 if (msglvl .gt. 0) then
802 diff -urN original/ARPACK/SRC/dnconv.f ARPACK/SRC/dnconv.f
803 --- original/ARPACK/SRC/dnconv.f 1996-04-20 18:27:05.000000000 +0200
804 +++ ARPACK/SRC/dnconv.f 2005-06-30 23:20:38.907789915 +0200
809 -c second ARPACK utility routine for timing.
810 +c arscnd ARPACK utility routine for timing.
811 c dlamch LAPACK routine that determines machine constants.
812 c dlapy2 LAPACK routine to compute sqrt(x**2+y**2) carefully.
815 c | for some appropriate choice of norm. |
816 c %-------------------------------------------------------------%
821 c %---------------------------------%
822 c | Get machine dependent constant. |
824 if (bounds(i) .le. tol*temp) nconv = nconv + 1
829 tnconv = tnconv + (t1 - t0)
832 diff -urN original/ARPACK/SRC/dneigh.f ARPACK/SRC/dneigh.f
833 --- original/ARPACK/SRC/dneigh.f 1996-04-20 18:27:06.000000000 +0200
834 +++ ARPACK/SRC/dneigh.f 2005-06-30 23:20:38.909789681 +0200
837 c dlaqrb ARPACK routine to compute the real Schur form of an
838 c upper Hessenberg matrix and last row of the Schur vectors.
839 -c second ARPACK utility routine for timing.
840 +c arscnd ARPACK utility routine for timing.
841 c dmout ARPACK utility routine that prints matrices
842 c dvout ARPACK utility routine that prints vectors.
843 c dlacpy LAPACK matrix copy routine.
845 c | External Subroutines |
846 c %----------------------%
848 - external dcopy, dlacpy, dlaqrb, dtrevc, dvout, second
849 + external dcopy, dlacpy, dlaqrb, dtrevc, dvout, arscnd
851 c %--------------------%
852 c | External Functions |
854 c | & message level for debugging |
855 c %-------------------------------%
861 if (msglvl .gt. 2) then
863 & '_neigh: Ritz estimates for the eigenvalues of H')
868 tneigh = tneigh + (t1 - t0)
871 diff -urN original/ARPACK/SRC/dngets.f ARPACK/SRC/dngets.f
872 --- original/ARPACK/SRC/dngets.f 1996-04-20 18:27:07.000000000 +0200
873 +++ ARPACK/SRC/dngets.f 2005-06-30 23:20:38.918788630 +0200
875 c | External Subroutines |
876 c %----------------------%
878 - external dcopy, dsortc, second
879 + external dcopy, dsortc, arscnd
881 c %----------------------%
882 c | Intrinsics Functions |
884 c | & message level for debugging |
885 c %-------------------------------%
891 c %----------------------------------------------------%
893 call dsortc ( 'SR', .true., np, bounds, ritzr, ritzi )
898 tngets = tngets + (t1 - t0)
900 if (msglvl .gt. 0) then
901 diff -urN original/ARPACK/SRC/dsaitr.f ARPACK/SRC/dsaitr.f
902 --- original/ARPACK/SRC/dsaitr.f 1996-08-28 16:21:43.000000000 +0200
903 +++ ARPACK/SRC/dsaitr.f 2005-06-30 23:20:38.923788047 +0200
905 c %----------------------%
907 external daxpy, dcopy, dscal, dgemv, dgetv0, dvout, dmout,
908 - & dlascl, ivout, second
909 + & dlascl, ivout, arscnd
911 c %--------------------%
912 c | External Functions |
914 c | & message level for debugging |
915 c %-------------------------------%
921 c %------------------------------%
923 c %------------------------------------------------%
928 tsaitr = tsaitr + (t1 - t0)
937 call dcopy (n, v(1,j), 1, workd(ivj), 1)
941 c | WORKD(IRJ:IRJ+N-1) := OP*v_{j}. |
942 c %-----------------------------------%
946 tmvopx = tmvopx + (t3 - t2)
950 c %-------------------------------------------%
952 if (mode .eq. 2) go to 65
955 if (bmat .eq. 'G') then
959 c %-----------------------------------%
961 if (bmat .eq. 'G') then
964 tmvbx = tmvbx + (t3 - t2)
967 @@ -592,12 +592,12 @@
979 if (bmat .eq. 'G') then
981 call dcopy (n, resid, 1, workd(irj), 1)
983 c %---------------------------------------------------%
985 if (bmat .eq. 'G') then
988 tmvbx = tmvbx + (t3 - t2)
992 h(j,2) = h(j,2) + workd(irj + j - 1)
997 if (bmat .eq. 'G') then
999 call dcopy (n, resid, 1, workd(irj), 1)
1001 c %---------------------------------------------------%
1003 if (bmat .eq. 'G') then
1006 tmvbx = tmvbx + (t3 - t2)
1015 titref = titref + (t5 - t4)
1017 c %----------------------------------------------------------%
1021 if (j .gt. k+np) then
1024 tsaitr = tsaitr + (t1 - t0)
1027 diff -urN original/ARPACK/SRC/dsapps.f ARPACK/SRC/dsapps.f
1028 --- original/ARPACK/SRC/dsapps.f 1998-05-20 16:58:59.000000000 +0200
1029 +++ ARPACK/SRC/dsapps.f 2005-06-30 23:20:38.926787696 +0200
1033 c ivout ARPACK utility routine that prints integers.
1034 -c second ARPACK utility routine for timing.
1035 +c arscnd ARPACK utility routine for timing.
1036 c dvout ARPACK utility routine that prints vectors.
1037 c dlamch LAPACK routine that determines machine constants.
1038 c dlartg LAPACK Givens rotation construction routine.
1040 c %----------------------%
1042 external daxpy, dcopy, dscal, dlacpy, dlartg, dlaset, dvout,
1043 - & ivout, second, dgemv
1044 + & ivout, arscnd, dgemv
1046 c %--------------------%
1047 c | External Functions |
1049 c | & message level for debugging |
1050 c %-------------------------------%
1063 tsapps = tsapps + (t1 - t0)
1066 diff -urN original/ARPACK/SRC/dsaup2.f ARPACK/SRC/dsaup2.f
1067 --- original/ARPACK/SRC/dsaup2.f 1998-05-20 16:59:00.000000000 +0200
1068 +++ ARPACK/SRC/dsaup2.f 2005-06-30 23:20:38.933786879 +0200
1070 c dsgets ARPACK reorder Ritz values and error bounds routine.
1071 c dsortr ARPACK sorting routine.
1072 c ivout ARPACK utility routine that prints integers.
1073 -c second ARPACK utility routine for timing.
1074 +c arscnd ARPACK utility routine for timing.
1075 c dvout ARPACK utility routine that prints vectors.
1076 c dlamch LAPACK routine that determines machine constants.
1077 c dcopy Level 1 BLAS that copies one vector to another.
1079 c %----------------------%
1081 external dcopy, dgetv0, dsaitr, dscal, dsconv, dseigt, dsgets,
1082 - & dsapps, dsortr, dvout, ivout, second, dswap
1083 + & dsapps, dsortr, dvout, ivout, arscnd, dswap
1085 c %--------------------%
1086 c | External Functions |
1088 c | & message level for debugging |
1089 c %-------------------------------%
1095 c %---------------------------------%
1097 c %---------------------------------------------%
1102 if (bmat .eq. 'G') then
1104 call dcopy (n, resid, 1, workd(n+1), 1)
1106 c %----------------------------------%
1108 if (bmat .eq. 'G') then
1111 tmvbx = tmvbx + (t3 - t2)
1123 --- ARPACK/SRC/dsaupd.f.orig 1996-08-27 07:28:56.000000000 +0200
1124 +++ ARPACK/SRC/dsaupd.f 2022-12-15 21:14:09.256722054 +0100
1125 @@ -380,7 +380,7 @@ c Arnoldi Iteration.
1126 c dstats ARPACK routine that initialize timing and other statistics
1128 c ivout ARPACK utility routine that prints integers.
1129 -c second ARPACK utility routine for timing.
1130 +c arscnd ARPACK utility routine for timing.
1131 c dvout ARPACK utility routine that prints vectors.
1132 c dlamch LAPACK routine that determines machine constants.
1134 @@ -456,7 +456,7 @@ c %----------------------%
1135 c | External Subroutines |
1136 c %----------------------%
1138 - external dsaup2, dvout, ivout, second, dstats
1139 + external dsaup2, dvout, ivout, arscnd, dstats
1141 c %--------------------%
1142 c | External Functions |
1143 @@ -478,7 +478,7 @@ c | & message level for debugging
1144 c %-------------------------------%
1152 @@ -637,7 +637,7 @@ c
1153 & '_saupd: corresponding error bounds')
1160 if (msglvl .gt. 0) then
1161 diff -urN original/ARPACK/SRC/dsconv.f ARPACK/SRC/dsconv.f
1162 --- original/ARPACK/SRC/dsconv.f 1996-04-20 18:27:11.000000000 +0200
1163 +++ ARPACK/SRC/dsconv.f 2005-06-30 23:20:38.939786179 +0200
1168 -c second ARPACK utility routine for timing.
1169 +c arscnd ARPACK utility routine for timing.
1170 c dlamch LAPACK routine that determines machine constants.
1174 c | Executable Statements |
1175 c %-----------------------%
1180 eps23 = dlamch('Epsilon-Machine')
1181 eps23 = eps23**(2.0D+0 / 3.0D+0)
1188 tsconv = tsconv + (t1 - t0)
1191 diff -urN original/ARPACK/SRC/dseigt.f ARPACK/SRC/dseigt.f
1192 --- original/ARPACK/SRC/dseigt.f 1996-08-27 07:29:04.000000000 +0200
1193 +++ ARPACK/SRC/dseigt.f 2005-06-30 23:20:38.940786062 +0200
1195 c dstqrb ARPACK routine that computes the eigenvalues and the
1196 c last components of the eigenvectors of a symmetric
1197 c and tridiagonal matrix.
1198 -c second ARPACK utility routine for timing.
1199 +c arscnd ARPACK utility routine for timing.
1200 c dvout ARPACK utility routine that prints vectors.
1201 c dcopy Level 1 BLAS that copies one vector to another.
1204 c | External Subroutines |
1205 c %----------------------%
1207 - external dcopy, dstqrb, dvout, second
1208 + external dcopy, dstqrb, dvout, arscnd
1210 c %-----------------------%
1211 c | Executable Statements |
1213 c | & message level for debugging |
1214 c %-------------------------------%
1220 if (msglvl .gt. 0) then
1222 bounds(k) = rnorm*abs(bounds(k))
1227 tseigt = tseigt + (t1 - t0)
1230 diff -urN original/ARPACK/SRC/dsgets.f ARPACK/SRC/dsgets.f
1231 --- original/ARPACK/SRC/dsgets.f 1996-04-20 18:27:14.000000000 +0200
1232 +++ ARPACK/SRC/dsgets.f 2005-06-30 23:20:38.949785011 +0200
1235 c dsortr ARPACK utility sorting routine.
1236 c ivout ARPACK utility routine that prints integers.
1237 -c second ARPACK utility routine for timing.
1238 +c arscnd ARPACK utility routine for timing.
1239 c dvout ARPACK utility routine that prints vectors.
1240 c dcopy Level 1 BLAS that copies one vector to another.
1241 c dswap Level 1 BLAS that swaps the contents of two vectors.
1243 c | External Subroutines |
1244 c %----------------------%
1246 - external dswap, dcopy, dsortr, second
1247 + external dswap, dcopy, dsortr, arscnd
1249 c %---------------------%
1250 c | Intrinsic Functions |
1252 c | & message level for debugging |
1253 c %-------------------------------%
1259 if (which .eq. 'BE') then
1261 call dcopy (np, ritz, 1, shifts, 1)
1266 tsgets = tsgets + (t1 - t0)
1268 if (msglvl .gt. 0) then
1269 diff -urN original/ARPACK/SRC/sgetv0.f ARPACK/SRC/sgetv0.f
1270 --- original/ARPACK/SRC/sgetv0.f 1999-04-08 00:42:46.000000000 +0200
1271 +++ ARPACK/SRC/sgetv0.f 2005-06-30 23:20:38.960783727 +0200
1273 c TR95-13, Department of Computational and Applied Mathematics.
1276 -c second ARPACK utility routine for timing.
1277 +c arscnd ARPACK utility routine for timing.
1278 c svout ARPACK utility routine for vector output.
1279 c slarnv LAPACK routine for generating a random vector.
1280 c sgemv Level 2 BLAS routine for matrix vector multiplication.
1282 c | External Subroutines |
1283 c %----------------------%
1285 - external slarnv, svout, scopy, sgemv, second
1286 + external slarnv, svout, scopy, sgemv, arscnd
1288 c %--------------------%
1289 c | External Functions |
1291 c | & message level for debugging |
1292 c %-------------------------------%
1300 c | the generalized problem when B is possibly (singular). |
1301 c %----------------------------------------------------------%
1305 if (bmat .eq. 'G') then
1311 if (bmat .eq. 'G') then
1314 tmvopx = tmvopx + (t3 - t2)
1318 c | Compute B-norm of starting vector. |
1319 c %------------------------------------------------------%
1324 if (bmat .eq. 'G') then
1329 if (bmat .eq. 'G') then
1332 tmvbx = tmvbx + (t3 - t2)
1336 c | Compute the B-norm of the orthogonalized starting vector |
1337 c %----------------------------------------------------------%
1341 if (bmat .eq. 'G') then
1343 call scopy (n, resid, 1, workd(n+1), 1)
1347 if (bmat .eq. 'G') then
1350 tmvbx = tmvbx + (t3 - t2)
1359 tgetv0 = tgetv0 + (t1 - t0)
1362 diff -urN original/ARPACK/SRC/snaitr.f ARPACK/SRC/snaitr.f
1363 --- original/ARPACK/SRC/snaitr.f 1996-08-27 18:09:01.000000000 +0200
1364 +++ ARPACK/SRC/snaitr.f 2005-06-30 23:20:38.969782676 +0200
1367 c sgetv0 ARPACK routine to generate the initial vector.
1368 c ivout ARPACK utility routine that prints integers.
1369 -c second ARPACK utility routine for timing.
1370 +c arscnd ARPACK utility routine for timing.
1371 c smout ARPACK utility routine that prints matrices
1372 c svout ARPACK utility routine that prints vectors.
1373 c slabad LAPACK routine that computes machine constants.
1375 c %----------------------%
1377 external saxpy, scopy, sscal, sgemv, sgetv0, slabad,
1378 - & svout, smout, ivout, second
1379 + & svout, smout, ivout, arscnd
1381 c %--------------------%
1382 c | External Functions |
1384 c | & message level for debugging |
1385 c %-------------------------------%
1391 c %------------------------------%
1393 c %------------------------------------------------%
1398 tnaitr = tnaitr + (t1 - t0)
1407 call scopy (n, v(1,j), 1, workd(ivj), 1)
1411 c | if step3 = .true. |
1412 c %----------------------------------%
1416 tmvopx = tmvopx + (t3 - t2)
1420 c | factorization to length j. |
1421 c %---------------------------------------%
1425 if (bmat .eq. 'G') then
1429 c %----------------------------------%
1431 if (bmat .eq. 'G') then
1434 tmvbx = tmvbx + (t3 - t2)
1437 @@ -576,11 +576,11 @@
1439 if (j .gt. 1) h(j,j-1) = betaj
1448 if (bmat .eq. 'G') then
1450 call scopy (n, resid, 1, workd(irj), 1)
1452 c %---------------------------------------------------%
1454 if (bmat .eq. 'G') then
1457 tmvbx = tmvbx + (t3 - t2)
1461 call saxpy (j, one, workd(irj), 1, h(1,j), 1)
1466 if (bmat .eq. 'G') then
1468 call scopy (n, resid, 1, workd(irj), 1)
1470 c %---------------------------------------------------%
1472 if (bmat .eq. 'G') then
1475 tmvbx = tmvbx + (t3 - t2)
1484 titref = titref + (t5 - t4)
1486 c %------------------------------------%
1490 if (j .gt. k+np) then
1493 tnaitr = tnaitr + (t1 - t0)
1495 do 110 i = max(1,k), k+np-1
1496 diff -urN original/ARPACK/SRC/snapps.f ARPACK/SRC/snapps.f
1497 --- original/ARPACK/SRC/snapps.f 1998-05-20 16:58:55.000000000 +0200
1498 +++ ARPACK/SRC/snapps.f 2005-06-30 23:20:38.973782209 +0200
1502 c ivout ARPACK utility routine that prints integers.
1503 -c second ARPACK utility routine for timing.
1504 +c arscnd ARPACK utility routine for timing.
1505 c smout ARPACK utility routine that prints matrices.
1506 c svout ARPACK utility routine that prints vectors.
1507 c slabad LAPACK routine that computes machine constants.
1509 c %----------------------%
1511 external saxpy, scopy, sscal, slacpy, slarfg, slarf,
1512 - & slaset, slabad, second, slartg
1513 + & slaset, slabad, arscnd, slartg
1515 c %--------------------%
1516 c | External Functions |
1518 c | & message level for debugging |
1519 c %-------------------------------%
1532 tnapps = tnapps + (t1 - t0)
1535 --- ARPACK/SRC/snaup2.f.orig 1996-08-16 22:28:10.000000000 +0200
1536 +++ ARPACK/SRC/snaup2.f 2022-12-15 21:16:41.849228723 +0100
1537 @@ -143,7 +143,7 @@ c sneigh ARPACK compute Ritz values
1538 c sngets ARPACK reorder Ritz values and error bounds routine.
1539 c ssortc ARPACK sorting routine.
1540 c ivout ARPACK utility routine that prints integers.
1541 -c second ARPACK utility routine for timing.
1542 +c arscnd ARPACK utility routine for timing.
1543 c smout ARPACK utility routine that prints matrices
1544 c svout ARPACK utility routine that prints vectors.
1545 c slamch LAPACK routine that determines machine constants.
1546 @@ -234,7 +234,7 @@ c | External Subroutines |
1547 c %----------------------%
1549 external scopy, sgetv0, snaitr, snconv, sneigh, sngets, snapps,
1550 - & svout, ivout, second
1551 + & svout, ivout, arscnd
1553 c %--------------------%
1554 c | External Functions |
1555 @@ -256,7 +256,7 @@ c %-----------------------%
1557 if (ido .eq. 0) then
1564 @@ -755,7 +755,7 @@ c | the first step of the next ca
1565 c %---------------------------------------------%
1570 if (bmat .eq. 'G') then
1572 call scopy (n, resid, 1, workd(n+1), 1)
1573 @@ -780,7 +780,7 @@ c | WORKD(1:N) := B*RESID
1574 c %----------------------------------%
1576 if (bmat .eq. 'G') then
1579 tmvbx = tmvbx + (t3 - t2)
1582 @@ -819,7 +819,7 @@ c %------------%
1591 diff -urN original/ARPACK/SRC/snaupd.f ARPACK/SRC/snaupd.f
1592 --- original/ARPACK/SRC/snaupd.f 2002-08-23 07:29:34.000000000 +0200
1593 +++ ARPACK/SRC/snaupd.f 2005-06-30 23:20:38.984780925 +0200
1595 c snaup2 ARPACK routine that implements the Implicitly Restarted
1596 c Arnoldi Iteration.
1597 c ivout ARPACK utility routine that prints integers.
1598 -c second ARPACK utility routine for timing.
1599 +c arscnd ARPACK utility routine for timing.
1600 c svout ARPACK utility routine that prints vectors.
1601 c slamch LAPACK routine that determines machine constants.
1604 c | External Subroutines |
1605 c %----------------------%
1607 - external snaup2, svout, ivout, second, sstatn
1608 + external snaup2, svout, ivout, arscnd, sstatn
1610 c %--------------------%
1611 c | External Functions |
1613 c %-------------------------------%
1620 c %----------------%
1622 & '_naupd: Associated Ritz estimates')
1629 if (msglvl .gt. 0) then
1630 diff -urN original/ARPACK/SRC/snconv.f ARPACK/SRC/snconv.f
1631 --- original/ARPACK/SRC/snconv.f 1996-04-20 18:27:21.000000000 +0200
1632 +++ ARPACK/SRC/snconv.f 2005-06-30 23:20:38.985780808 +0200
1637 -c second ARPACK utility routine for timing.
1638 +c arscnd ARPACK utility routine for timing.
1639 c slamch LAPACK routine that determines machine constants.
1640 c slapy2 LAPACK routine to compute sqrt(x**2+y**2) carefully.
1643 c | for some appropriate choice of norm. |
1644 c %-------------------------------------------------------------%
1649 c %---------------------------------%
1650 c | Get machine dependent constant. |
1652 if (bounds(i) .le. tol*temp) nconv = nconv + 1
1657 tnconv = tnconv + (t1 - t0)
1660 diff -urN original/ARPACK/SRC/sneigh.f ARPACK/SRC/sneigh.f
1661 --- original/ARPACK/SRC/sneigh.f 1996-04-20 18:27:22.000000000 +0200
1662 +++ ARPACK/SRC/sneigh.f 2005-06-30 23:20:38.987780575 +0200
1665 c slaqrb ARPACK routine to compute the real Schur form of an
1666 c upper Hessenberg matrix and last row of the Schur vectors.
1667 -c second ARPACK utility routine for timing.
1668 +c arscnd ARPACK utility routine for timing.
1669 c smout ARPACK utility routine that prints matrices
1670 c svout ARPACK utility routine that prints vectors.
1671 c slacpy LAPACK matrix copy routine.
1673 c | External Subroutines |
1674 c %----------------------%
1676 - external scopy, slacpy, slaqrb, strevc, svout, second
1677 + external scopy, slacpy, slaqrb, strevc, svout, arscnd
1679 c %--------------------%
1680 c | External Functions |
1682 c | & message level for debugging |
1683 c %-------------------------------%
1689 if (msglvl .gt. 2) then
1691 & '_neigh: Ritz estimates for the eigenvalues of H')
1696 tneigh = tneigh + (t1 - t0)
1699 diff -urN original/ARPACK/SRC/sngets.f ARPACK/SRC/sngets.f
1700 --- original/ARPACK/SRC/sngets.f 1996-04-20 18:27:24.000000000 +0200
1701 +++ ARPACK/SRC/sngets.f 2005-06-30 23:20:38.995779641 +0200
1703 c | External Subroutines |
1704 c %----------------------%
1706 - external scopy, ssortc, second
1707 + external scopy, ssortc, arscnd
1709 c %----------------------%
1710 c | Intrinsics Functions |
1712 c | & message level for debugging |
1713 c %-------------------------------%
1719 c %----------------------------------------------------%
1721 call ssortc ( 'SR', .true., np, bounds, ritzr, ritzi )
1726 tngets = tngets + (t1 - t0)
1728 if (msglvl .gt. 0) then
1729 diff -urN original/ARPACK/SRC/ssaitr.f ARPACK/SRC/ssaitr.f
1730 --- original/ARPACK/SRC/ssaitr.f 1996-08-28 16:21:43.000000000 +0200
1731 +++ ARPACK/SRC/ssaitr.f 2005-06-30 23:20:39.001778940 +0200
1733 c %----------------------%
1735 external saxpy, scopy, sscal, sgemv, sgetv0, svout, smout,
1736 - & slascl, ivout, second
1737 + & slascl, ivout, arscnd
1739 c %--------------------%
1740 c | External Functions |
1742 c | & message level for debugging |
1743 c %-------------------------------%
1749 c %------------------------------%
1751 c %------------------------------------------------%
1756 tsaitr = tsaitr + (t1 - t0)
1765 call scopy (n, v(1,j), 1, workd(ivj), 1)
1769 c | WORKD(IRJ:IRJ+N-1) := OP*v_{j}. |
1770 c %-----------------------------------%
1774 tmvopx = tmvopx + (t3 - t2)
1778 c %-------------------------------------------%
1780 if (mode .eq. 2) go to 65
1783 if (bmat .eq. 'G') then
1787 c %-----------------------------------%
1789 if (bmat .eq. 'G') then
1792 tmvbx = tmvbx + (t3 - t2)
1795 @@ -592,12 +592,12 @@
1807 if (bmat .eq. 'G') then
1809 call scopy (n, resid, 1, workd(irj), 1)
1811 c %---------------------------------------------------%
1813 if (bmat .eq. 'G') then
1816 tmvbx = tmvbx + (t3 - t2)
1820 h(j,2) = h(j,2) + workd(irj + j - 1)
1825 if (bmat .eq. 'G') then
1827 call scopy (n, resid, 1, workd(irj), 1)
1829 c %---------------------------------------------------%
1831 if (bmat .eq. 'G') then
1834 tmvbx = tmvbx + (t3 - t2)
1843 titref = titref + (t5 - t4)
1845 c %----------------------------------------------------------%
1849 if (j .gt. k+np) then
1852 tsaitr = tsaitr + (t1 - t0)
1855 diff -urN original/ARPACK/SRC/ssapps.f ARPACK/SRC/ssapps.f
1856 --- original/ARPACK/SRC/ssapps.f 1998-05-20 16:58:58.000000000 +0200
1857 +++ ARPACK/SRC/ssapps.f 2005-06-30 23:20:39.004778590 +0200
1861 c ivout ARPACK utility routine that prints integers.
1862 -c second ARPACK utility routine for timing.
1863 +c arscnd ARPACK utility routine for timing.
1864 c svout ARPACK utility routine that prints vectors.
1865 c slamch LAPACK routine that determines machine constants.
1866 c slartg LAPACK Givens rotation construction routine.
1868 c %----------------------%
1870 external saxpy, scopy, sscal, slacpy, slartg, slaset, svout,
1871 - & ivout, second, sgemv
1872 + & ivout, arscnd, sgemv
1874 c %--------------------%
1875 c | External Functions |
1877 c | & message level for debugging |
1878 c %-------------------------------%
1891 tsapps = tsapps + (t1 - t0)
1894 diff -urN original/ARPACK/SRC/ssaup2.f ARPACK/SRC/ssaup2.f
1895 --- original/ARPACK/SRC/ssaup2.f 1998-05-20 16:58:59.000000000 +0200
1896 +++ ARPACK/SRC/ssaup2.f 2005-06-30 23:20:39.010777890 +0200
1898 c ssgets ARPACK reorder Ritz values and error bounds routine.
1899 c ssortr ARPACK sorting routine.
1900 c ivout ARPACK utility routine that prints integers.
1901 -c second ARPACK utility routine for timing.
1902 +c arscnd ARPACK utility routine for timing.
1903 c svout ARPACK utility routine that prints vectors.
1904 c slamch LAPACK routine that determines machine constants.
1905 c scopy Level 1 BLAS that copies one vector to another.
1907 c %----------------------%
1909 external scopy, sgetv0, ssaitr, sscal, ssconv, sseigt, ssgets,
1910 - & ssapps, ssortr, svout, ivout, second, sswap
1911 + & ssapps, ssortr, svout, ivout, arscnd, sswap
1913 c %--------------------%
1914 c | External Functions |
1916 c | & message level for debugging |
1917 c %-------------------------------%
1923 c %---------------------------------%
1925 c %---------------------------------------------%
1930 if (bmat .eq. 'G') then
1932 call scopy (n, resid, 1, workd(n+1), 1)
1934 c %----------------------------------%
1936 if (bmat .eq. 'G') then
1939 tmvbx = tmvbx + (t3 - t2)
1951 diff -urN original/ARPACK/SRC/ssaupd.f ARPACK/SRC/ssaupd.f
1952 --- original/ARPACK/SRC/ssaupd.f 2001-04-11 00:52:40.000000000 +0200
1953 +++ ARPACK/SRC/ssaupd.f 2005-06-30 23:20:39.015777306 +0200
1955 c sstats ARPACK routine that initialize timing and other statistics
1957 c ivout ARPACK utility routine that prints integers.
1958 -c second ARPACK utility routine for timing.
1959 +c arscnd ARPACK utility routine for timing.
1960 c svout ARPACK utility routine that prints vectors.
1961 c slamch LAPACK routine that determines machine constants.
1964 c | External Subroutines |
1965 c %----------------------%
1967 - external ssaup2, svout, ivout, second, sstats
1968 + external ssaup2, svout, ivout, arscnd, sstats
1970 c %--------------------%
1971 c | External Functions |
1973 c %-------------------------------%
1982 & '_saupd: corresponding error bounds')
1989 if (msglvl .gt. 0) then
1990 diff -urN original/ARPACK/SRC/ssconv.f ARPACK/SRC/ssconv.f
1991 --- original/ARPACK/SRC/ssconv.f 1996-04-20 18:27:27.000000000 +0200
1992 +++ ARPACK/SRC/ssconv.f 2005-06-30 23:20:39.016777190 +0200
1997 -c second ARPACK utility routine for timing.
1998 +c arscnd ARPACK utility routine for timing.
1999 c slamch LAPACK routine that determines machine constants.
2003 c | Executable Statements |
2004 c %-----------------------%
2009 eps23 = slamch('Epsilon-Machine')
2010 eps23 = eps23**(2.0E+0 / 3.0E+0)
2017 tsconv = tsconv + (t1 - t0)
2020 diff -urN original/ARPACK/SRC/sseigt.f ARPACK/SRC/sseigt.f
2021 --- original/ARPACK/SRC/sseigt.f 1996-08-27 07:29:04.000000000 +0200
2022 +++ ARPACK/SRC/sseigt.f 2005-06-30 23:20:39.017777073 +0200
2024 c sstqrb ARPACK routine that computes the eigenvalues and the
2025 c last components of the eigenvectors of a symmetric
2026 c and tridiagonal matrix.
2027 -c second ARPACK utility routine for timing.
2028 +c arscnd ARPACK utility routine for timing.
2029 c svout ARPACK utility routine that prints vectors.
2030 c scopy Level 1 BLAS that copies one vector to another.
2033 c | External Subroutines |
2034 c %----------------------%
2036 - external scopy, sstqrb, svout, second
2037 + external scopy, sstqrb, svout, arscnd
2039 c %-----------------------%
2040 c | Executable Statements |
2042 c | & message level for debugging |
2043 c %-------------------------------%
2049 if (msglvl .gt. 0) then
2051 bounds(k) = rnorm*abs(bounds(k))
2056 tseigt = tseigt + (t1 - t0)
2059 diff -urN original/ARPACK/SRC/ssgets.f ARPACK/SRC/ssgets.f
2060 --- original/ARPACK/SRC/ssgets.f 1996-04-20 18:27:30.000000000 +0200
2061 +++ ARPACK/SRC/ssgets.f 2005-06-30 23:20:39.026776022 +0200
2064 c ssortr ARPACK utility sorting routine.
2065 c ivout ARPACK utility routine that prints integers.
2066 -c second ARPACK utility routine for timing.
2067 +c arscnd ARPACK utility routine for timing.
2068 c svout ARPACK utility routine that prints vectors.
2069 c scopy Level 1 BLAS that copies one vector to another.
2070 c sswap Level 1 BLAS that swaps the contents of two vectors.
2072 c | External Subroutines |
2073 c %----------------------%
2075 - external sswap, scopy, ssortr, second
2076 + external sswap, scopy, ssortr, arscnd
2078 c %---------------------%
2079 c | Intrinsic Functions |
2081 c | & message level for debugging |
2082 c %-------------------------------%
2088 if (which .eq. 'BE') then
2090 call scopy (np, ritz, 1, shifts, 1)
2095 tsgets = tsgets + (t1 - t0)
2097 if (msglvl .gt. 0) then
2098 diff -urN original/ARPACK/SRC/zgetv0.f ARPACK/SRC/zgetv0.f
2099 --- original/ARPACK/SRC/zgetv0.f 1999-04-08 00:46:45.000000000 +0200
2100 +++ ARPACK/SRC/zgetv0.f 2005-06-30 23:20:39.037774738 +0200
2105 -c second ARPACK utility routine for timing.
2106 +c arscnd ARPACK utility routine for timing.
2107 c zvout ARPACK utility routine that prints vectors.
2108 c zlarnv LAPACK routine for generating a random vector.
2109 c zgemv Level 2 BLAS routine for matrix vector multiplication.
2111 c | External Subroutines |
2112 c %----------------------%
2114 - external zcopy, zgemv, zlarnv, zvout, second
2115 + external zcopy, zgemv, zlarnv, zvout, arscnd
2117 c %--------------------%
2118 c | External Functions |
2120 c | & message level for debugging |
2121 c %-------------------------------%
2129 c | the generalized problem when B is possibly (singular). |
2130 c %----------------------------------------------------------%
2134 if (bmat .eq. 'G') then
2143 tmvopx = tmvopx + (t3 - t2)
2145 c %------------------------------------------------------%
2147 c | Compute B-norm of starting vector. |
2148 c %------------------------------------------------------%
2153 if (bmat .eq. 'G') then
2158 if (bmat .eq. 'G') then
2161 tmvbx = tmvbx + (t3 - t2)
2165 c | Compute the B-norm of the orthogonalized starting vector |
2166 c %----------------------------------------------------------%
2170 if (bmat .eq. 'G') then
2172 call zcopy (n, resid, 1, workd(n+1), 1)
2176 if (bmat .eq. 'G') then
2179 tmvbx = tmvbx + (t3 - t2)
2188 tgetv0 = tgetv0 + (t1 - t0)
2191 diff -urN original/ARPACK/SRC/znaitr.f ARPACK/SRC/znaitr.f
2192 --- original/ARPACK/SRC/znaitr.f 1996-08-27 18:12:57.000000000 +0200
2193 +++ ARPACK/SRC/znaitr.f 2005-06-30 23:20:39.042774155 +0200
2196 c zgetv0 ARPACK routine to generate the initial vector.
2197 c ivout ARPACK utility routine that prints integers.
2198 -c second ARPACK utility routine for timing.
2199 +c arscnd ARPACK utility routine for timing.
2200 c zmout ARPACK utility routine that prints matrices
2201 c zvout ARPACK utility routine that prints vectors.
2202 c zlanhs LAPACK routine that computes various norms of a matrix.
2204 c %----------------------%
2206 external zaxpy, zcopy, zscal, zdscal, zgemv, zgetv0,
2207 - & dlabad, zvout, zmout, ivout, second
2208 + & dlabad, zvout, zmout, ivout, arscnd
2210 c %--------------------%
2211 c | External Functions |
2213 c | & message level for debugging |
2214 c %-------------------------------%
2220 c %------------------------------%
2222 c %------------------------------------------------%
2227 tcaitr = tcaitr + (t1 - t0)
2236 call zcopy (n, v(1,j), 1, workd(ivj), 1)
2240 c | if step3 = .true. |
2241 c %----------------------------------%
2245 tmvopx = tmvopx + (t3 - t2)
2249 c | factorization to length j. |
2250 c %---------------------------------------%
2254 if (bmat .eq. 'G') then
2258 c %----------------------------------%
2260 if (bmat .eq. 'G') then
2263 tmvbx = tmvbx + (t3 - t2)
2266 @@ -583,11 +583,11 @@
2268 if (j .gt. 1) h(j,j-1) = dcmplx(betaj, rzero)
2277 if (bmat .eq. 'G') then
2279 call zcopy (n, resid, 1, workd(irj), 1)
2281 c %---------------------------------------------------%
2283 if (bmat .eq. 'G') then
2286 tmvbx = tmvbx + (t3 - t2)
2290 call zaxpy (j, one, workd(irj), 1, h(1,j), 1)
2295 if (bmat .eq. 'G') then
2297 call zcopy (n, resid, 1, workd(irj), 1)
2299 c %---------------------------------------------------%
2301 if (bmat .eq. 'G') then
2304 tmvbx = tmvbx + (t3 - t2)
2313 titref = titref + (t5 - t4)
2315 c %------------------------------------%
2319 if (j .gt. k+np) then
2322 tcaitr = tcaitr + (t1 - t0)
2324 do 110 i = max(1,k), k+np-1
2325 diff -urN original/ARPACK/SRC/znapps.f ARPACK/SRC/znapps.f
2326 --- original/ARPACK/SRC/znapps.f 1998-05-20 16:45:03.000000000 +0200
2327 +++ ARPACK/SRC/znapps.f 2005-06-30 23:20:39.045773805 +0200
2331 c ivout ARPACK utility routine that prints integers.
2332 -c second ARPACK utility routine for timing.
2333 +c arscnd ARPACK utility routine for timing.
2334 c zmout ARPACK utility routine that prints matrices
2335 c zvout ARPACK utility routine that prints vectors.
2336 c zlacpy LAPACK matrix copy routine.
2338 c %----------------------%
2340 external zaxpy, zcopy, zgemv, zscal, zlacpy, zlartg,
2341 - & zvout, zlaset, dlabad, zmout, second, ivout
2342 + & zvout, zlaset, dlabad, zmout, arscnd, ivout
2344 c %--------------------%
2345 c | External Functions |
2347 c | & message level for debugging |
2348 c %-------------------------------%
2361 tcapps = tcapps + (t1 - t0)
2364 --- ARPACK/SRC/znaup2.f.orig 1996-09-24 22:40:01.000000000 +0200
2365 +++ ARPACK/SRC/znaup2.f 2022-12-15 21:20:55.564520897 +0100
2366 @@ -135,7 +135,7 @@ c zneigh ARPACK compute Ritz values
2367 c zngets ARPACK reorder Ritz values and error bounds routine.
2368 c zsortc ARPACK sorting routine.
2369 c ivout ARPACK utility routine that prints integers.
2370 -c second ARPACK utility routine for timing.
2371 +c arscnd ARPACK utility routine for timing.
2372 c zmout ARPACK utility routine that prints matrices
2373 c zvout ARPACK utility routine that prints vectors.
2374 c dvout ARPACK utility routine that prints vectors.
2375 @@ -239,7 +239,7 @@ c | External Subroutines |
2376 c %----------------------%
2378 external zcopy, zgetv0, znaitr, zneigh, zngets, znapps,
2379 - & zsortc, zswap, zmout, zvout, ivout, second
2380 + & zsortc, zswap, zmout, zvout, ivout, arscnd
2382 c %--------------------%
2383 c | External functions |
2384 @@ -263,7 +263,7 @@ c %-----------------------%
2386 if (ido .eq. 0) then
2393 @@ -723,7 +723,7 @@ c | the first step of the next ca
2394 c %---------------------------------------------%
2399 if (bmat .eq. 'G') then
2401 call zcopy (n, resid, 1, workd(n+1), 1)
2402 @@ -748,7 +748,7 @@ c | WORKD(1:N) := B*RESID
2403 c %----------------------------------%
2405 if (bmat .eq. 'G') then
2408 tmvbx = tmvbx + (t3 - t2)
2411 @@ -787,7 +787,7 @@ c %------------%
2420 diff -urN original/ARPACK/SRC/znaupd.f ARPACK/SRC/znaupd.f
2421 --- original/ARPACK/SRC/znaupd.f 2002-08-15 07:50:57.000000000 +0200
2422 +++ ARPACK/SRC/znaupd.f 2005-06-30 23:20:39.055772637 +0200
2424 c zstatn ARPACK routine that initializes the timing variables.
2425 c ivout ARPACK utility routine that prints integers.
2426 c zvout ARPACK utility routine that prints vectors.
2427 -c second ARPACK utility routine for timing.
2428 +c arscnd ARPACK utility routine for timing.
2429 c dlamch LAPACK routine that determines machine constants.
2433 c | External Subroutines |
2434 c %----------------------%
2436 - external znaup2, zvout, ivout, second, zstatn
2437 + external znaup2, zvout, ivout, arscnd, zstatn
2439 c %--------------------%
2440 c | External Functions |
2442 c %-------------------------------%
2449 c %----------------%
2451 & '_naupd: Associated Ritz estimates')
2458 if (msglvl .gt. 0) then
2459 diff -urN original/ARPACK/SRC/zneigh.f ARPACK/SRC/zneigh.f
2460 --- original/ARPACK/SRC/zneigh.f 1996-04-20 19:16:01.000000000 +0200
2461 +++ ARPACK/SRC/zneigh.f 2005-06-30 23:20:39.056772521 +0200
2465 c ivout ARPACK utility routine that prints integers.
2466 -c second ARPACK utility routine for timing.
2467 +c arscnd ARPACK utility routine for timing.
2468 c zmout ARPACK utility routine that prints matrices
2469 c zvout ARPACK utility routine that prints vectors.
2470 c dvout ARPACK utility routine that prints vectors.
2472 c %----------------------%
2474 external zlacpy, zlahqr, ztrevc, zcopy,
2475 - & zdscal, zmout, zvout, second
2476 + & zdscal, zmout, zvout, arscnd
2478 c %--------------------%
2479 c | External Functions |
2481 c | & message level for debugging |
2482 c %-------------------------------%
2488 if (msglvl .gt. 2) then
2490 & '_neigh: Ritz estimates for the eigenvalues of H')
2495 tceigh = tceigh + (t1 - t0)
2498 diff -urN original/ARPACK/SRC/zngets.f ARPACK/SRC/zngets.f
2499 --- original/ARPACK/SRC/zngets.f 1996-04-20 19:16:02.000000000 +0200
2500 +++ ARPACK/SRC/zngets.f 2005-06-30 23:20:39.063771704 +0200
2503 c zsortc ARPACK sorting routine.
2504 c ivout ARPACK utility routine that prints integers.
2505 -c second ARPACK utility routine for timing.
2506 +c arscnd ARPACK utility routine for timing.
2507 c zvout ARPACK utility routine that prints vectors.
2511 c | External Subroutines |
2512 c %----------------------%
2514 - external zvout, zsortc, second
2515 + external zvout, zsortc, arscnd
2517 c %-----------------------%
2518 c | Executable Statements |
2520 c | & message level for debugging |
2521 c %-------------------------------%
2527 call zsortc (which, .true., kev+np, ritz, bounds)
2534 tcgets = tcgets + (t1 - t0)
2536 if (msglvl .gt. 0) then
2537 diff -urN original/ARPACK/UTIL/Makefile ARPACK/UTIL/Makefile
2538 --- original/ARPACK/UTIL/Makefile 1996-09-09 16:30:48.000000000 +0200
2539 +++ ARPACK/UTIL/Makefile 2005-06-30 23:23:46.746224656 +0200
2541 # ../ARmake.inc and is created at the next higher directory level.
2544 -OBJS = icnteq.o icopy.o iset.o iswap.o ivout.o second.o
2545 +OBJS = icnteq.o icopy.o iset.o iswap.o ivout.o arscnd.o
2547 SOBJ = svout.o smout.o
2549 diff -urN original/ARPACK/UTIL/arscnd.f ARPACK/UTIL/arscnd.f
2550 --- original/ARPACK/UTIL/arscnd.f 1970-01-01 01:00:00.000000000 +0100
2551 +++ ARPACK/UTIL/arscnd.f 2005-06-30 23:22:55.982983811 +0200
2553 + SUBROUTINE ARSCND( T )
2557 +* -- LAPACK auxiliary routine (preliminary version) --
2558 +* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
2559 +* Courant Institute, Argonne National Lab, and Rice University
2565 +* ARSCND returns the user time for a process in seconds.
2566 +* This version gets the time from the system function ETIME.
2568 +* .. Local Scalars ..
2571 +* .. Local Arrays ..
2574 +* .. External Functions ..
2578 +* .. Executable Statements ..
2581 + T1 = ETIME( TARRAY )
2589 diff -urN original/ARPACK/UTIL/second.f ARPACK/UTIL/second.f
2590 --- original/ARPACK/UTIL/second.f 1995-10-02 19:56:44.000000000 +0100
2591 +++ ARPACK/UTIL/second.f 1970-01-01 01:00:00.000000000 +0100
2593 - SUBROUTINE SECOND( T )
2597 -* -- LAPACK auxiliary routine (preliminary version) --
2598 -* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
2599 -* Courant Institute, Argonne National Lab, and Rice University
2605 -* SECOND returns the user time for a process in seconds.
2606 -* This version gets the time from the system function ETIME.
2608 -* .. Local Scalars ..
2611 -* .. Local Arrays ..
2614 -* .. External Functions ..
2618 -* .. Executable Statements ..
2621 - T1 = ETIME( TARRAY )