diff -urN original/ARPACK/SRC/cgetv0.f ARPACK/SRC/cgetv0.f
--- original/ARPACK/SRC/cgetv0.f	1999-04-08 00:46:44.000000000 +0200
+++ ARPACK/SRC/cgetv0.f	2005-06-30 23:20:38.850796569 +0200
@@ -89,7 +89,7 @@
 c     pp 357-385.
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     cvout   ARPACK utility routine that prints vectors.
 c     clarnv  LAPACK routine for generating a random vector. 
 c     cgemv   Level 2 BLAS routine for matrix vector multiplication.
@@ -168,7 +168,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   ccopy, cgemv, clarnv, cvout, second
+      external   ccopy, cgemv, clarnv, cvout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -211,7 +211,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mgetv0
 c 
          ierr   = 0
@@ -238,7 +238,7 @@
 c        | the generalized problem when B is possibly (singular).   |
 c        %----------------------------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nopx = nopx + 1
             ipntr(1) = 1
@@ -261,7 +261,7 @@
 c
       if (orth)  go to 40
 c 
-      call second (t3)
+      call arscnd (t3)
       tmvopx = tmvopx + (t3 - t2)
 c 
 c     %------------------------------------------------------%
@@ -269,7 +269,7 @@
 c     | Compute B-norm of starting vector.                   |
 c     %------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       first = .TRUE.
       if (bmat .eq. 'G') then
          nbx = nbx + 1
@@ -285,7 +285,7 @@
    20 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -328,7 +328,7 @@
 c     | Compute the B-norm of the orthogonalized starting vector |
 c     %----------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       if (bmat .eq. 'G') then
          nbx = nbx + 1
          call ccopy (n, resid, 1, workd(n+1), 1)
@@ -343,7 +343,7 @@
    40 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -401,7 +401,7 @@
       end if
       ido = 99
 c 
-      call second (t1)
+      call arscnd (t1)
       tgetv0 = tgetv0 + (t1 - t0)
 c 
  9000 continue
diff -urN original/ARPACK/SRC/cnaitr.f ARPACK/SRC/cnaitr.f
--- original/ARPACK/SRC/cnaitr.f	1996-08-27 18:12:56.000000000 +0200
+++ ARPACK/SRC/cnaitr.f	2005-06-30 23:20:38.856795869 +0200
@@ -131,7 +131,7 @@
 c\Routines called:
 c     cgetv0  ARPACK routine to generate the initial vector.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     cmout   ARPACK utility routine that prints matrices
 c     cvout   ARPACK utility routine that prints vectors.
 c     clanhs  LAPACK routine that computes various norms of a matrix.
@@ -273,7 +273,7 @@
 c     %----------------------%
 c
       external   caxpy, ccopy, cscal, csscal, cgemv, cgetv0, 
-     &           slabad, cvout, cmout, ivout, second
+     &           slabad, cvout, cmout, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -326,7 +326,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mcaitr
 c 
 c        %------------------------------%
@@ -437,7 +437,7 @@
 c              %------------------------------------------------%
 c
                info = j - 1
-               call second (t1)
+               call arscnd (t1)
                tcaitr = tcaitr + (t1 - t0)
                ido = 99
                go to 9000
@@ -477,7 +477,7 @@
 c
          step3 = .true.
          nopx  = nopx + 1
-         call second (t2)
+         call arscnd (t2)
          call ccopy (n, v(1,j), 1, workd(ivj), 1)
          ipntr(1) = ivj
          ipntr(2) = irj
@@ -497,7 +497,7 @@
 c        | if step3 = .true.                |
 c        %----------------------------------%
 c
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
  
          step3 = .false.
@@ -513,7 +513,7 @@
 c        |          factorization to length j.   |
 c        %---------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             step4 = .true.
@@ -538,7 +538,7 @@
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -583,11 +583,11 @@
 c
          if (j .gt. 1) h(j,j-1) = cmplx(betaj, rzero)
 c
-         call second (t4)
+         call arscnd (t4)
 c 
          orth1 = .true.
 c 
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call ccopy (n, resid, 1, workd(irj), 1)
@@ -611,7 +611,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -689,7 +689,7 @@
          call caxpy (j, one, workd(irj), 1, h(1,j), 1)
 c 
          orth2 = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call ccopy (n, resid, 1, workd(irj), 1)
@@ -713,7 +713,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if 
 c
@@ -791,7 +791,7 @@
          rstart = .false.
          orth2  = .false.
 c 
-         call second (t5)
+         call arscnd (t5)
          titref = titref + (t5 - t4)
 c 
 c        %------------------------------------%
@@ -800,7 +800,7 @@
 c
          j = j + 1
          if (j .gt. k+np) then
-            call second (t1)
+            call arscnd (t1)
             tcaitr = tcaitr + (t1 - t0)
             ido = 99
             do 110 i = max(1,k), k+np-1
diff -urN original/ARPACK/SRC/cnapps.f ARPACK/SRC/cnapps.f
--- original/ARPACK/SRC/cnapps.f	1998-05-20 16:45:03.000000000 +0200
+++ ARPACK/SRC/cnapps.f	2005-06-30 23:20:38.859795518 +0200
@@ -92,7 +92,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     cmout   ARPACK utility routine that prints matrices
 c     cvout   ARPACK utility routine that prints vectors.
 c     clacpy  LAPACK matrix copy routine.
@@ -184,7 +184,7 @@
 c     %----------------------%
 c
       external   caxpy, ccopy, cgemv, cscal, clacpy, clartg, 
-     &           cvout, claset, slabad, cmout, second, ivout
+     &           cvout, claset, slabad, cmout, arscnd, ivout
 c
 c     %--------------------%
 c     | External Functions |
@@ -240,7 +240,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mcapps
 c 
       kplusp = kev + np 
@@ -495,7 +495,7 @@
       end if
 c
  9000 continue
-      call second (t1)
+      call arscnd (t1)
       tcapps = tcapps + (t1 - t0)
 c 
       return
diff -urN original/ARPACK/SRC/cnaup2.f ARPACK/SRC/cnaup2.f
--- original/ARPACK/SRC/cnaup2.f	2000-06-01 22:28:53.000000000 +0200
+++ ARPACK/SRC/cnaup2.f	2005-06-30 23:20:38.864794935 +0200
@@ -135,7 +135,7 @@
 c     cngets  ARPACK reorder Ritz values and error bounds routine.
 c     csortc  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     cmout   ARPACK utility routine that prints matrices
 c     cvout   ARPACK utility routine that prints vectors.
 c     svout   ARPACK utility routine that prints vectors.
@@ -240,7 +240,7 @@
 c     %----------------------%
 c
       external   ccopy, cgetv0, cnaitr, cneigh, cngets, cnapps,
-     &           csortc, cswap, cmout, cvout, ivout, second
+     &           csortc, cswap, cmout, cvout, ivout, arscnd
 c
 c     %--------------------%
 c     | External functions |
@@ -264,7 +264,7 @@
 c
       if (ido .eq. 0) then
 c 
-         call second (t0)
+         call arscnd (t0)
 c 
          msglvl = mcaup2
 c 
@@ -724,7 +724,7 @@
 c        %---------------------------------------------%
 c
          cnorm = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call ccopy (n, resid, 1, workd(n+1), 1)
@@ -749,7 +749,7 @@
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -788,7 +788,7 @@
 c     | Error Exit |
 c     %------------%
 c
-      call second (t1)
+      call arscnd (t1)
       tcaup2 = t1 - t0
 c     
  9000 continue
diff -urN original/ARPACK/SRC/cnaupd.f ARPACK/SRC/cnaupd.f
--- original/ARPACK/SRC/cnaupd.f	2002-08-15 07:50:57.000000000 +0200
+++ ARPACK/SRC/cnaupd.f	2005-06-30 23:20:38.869794351 +0200
@@ -356,7 +356,7 @@
 c     cstatn  ARPACK routine that initializes the timing variables.
 c     ivout   ARPACK utility routine that prints integers.
 c     cvout   ARPACK utility routine that prints vectors.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     slamch  LAPACK routine that determines machine constants.
 c
 c\Author
@@ -429,7 +429,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   cnaup2, cvout, ivout, second, cstatn
+      external   cnaup2, cvout, ivout, arscnd, cstatn
 c
 c     %--------------------%
 c     | External Functions |
@@ -451,7 +451,7 @@
 c        %-------------------------------%
 c
          call cstatn
-         call second (t0)
+         call arscnd (t0)
          msglvl = mcaupd
 c
 c        %----------------%
@@ -611,7 +611,7 @@
      &               '_naupd: Associated Ritz estimates')
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tcaupd = t1 - t0
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/cneigh.f ARPACK/SRC/cneigh.f
--- original/ARPACK/SRC/cneigh.f	1996-04-20 19:15:53.000000000 +0200
+++ ARPACK/SRC/cneigh.f	2005-06-30 23:20:38.870794234 +0200
@@ -64,7 +64,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     cmout   ARPACK utility routine that prints matrices
 c     cvout   ARPACK utility routine that prints vectors.
 c     svout   ARPACK utility routine that prints vectors.
@@ -152,7 +152,7 @@
 c     %----------------------%
 c
       external   clacpy, clahqr, ctrevc, ccopy, 
-     &           csscal, cmout, cvout, second
+     &           csscal, cmout, cvout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -171,7 +171,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mceigh
 c 
       if (msglvl .gt. 2) then
@@ -244,7 +244,7 @@
      &              '_neigh: Ritz estimates for the eigenvalues of H')
       end if
 c
-      call second(t1)
+      call arscnd(t1)
       tceigh = tceigh + (t1 - t0)
 c
  9000 continue
diff -urN original/ARPACK/SRC/cngets.f ARPACK/SRC/cngets.f
--- original/ARPACK/SRC/cngets.f	1996-04-20 19:15:56.000000000 +0200
+++ ARPACK/SRC/cngets.f	2005-06-30 23:20:38.877793417 +0200
@@ -63,7 +63,7 @@
 c\Routines called:
 c     csortc  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     cvout   ARPACK utility routine that prints vectors.
 c
 c\Author
@@ -126,7 +126,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   cvout,  csortc, second
+      external   cvout,  csortc, arscnd
 c
 c     %-----------------------%
 c     | Executable Statements |
@@ -137,7 +137,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c 
-      call second (t0)
+      call arscnd (t0)
       msglvl = mcgets
 c 
       call csortc (which, .true., kev+np, ritz, bounds)
@@ -157,7 +157,7 @@
 c
       end if
 c     
-      call second (t1)
+      call arscnd (t1)
       tcgets = tcgets + (t1 - t0)
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/dgetv0.f ARPACK/SRC/dgetv0.f
--- original/ARPACK/SRC/dgetv0.f	1999-04-08 00:42:47.000000000 +0200
+++ ARPACK/SRC/dgetv0.f	2005-06-30 23:20:38.882792833 +0200
@@ -93,7 +93,7 @@
 c     TR95-13, Department of Computational and Applied Mathematics.
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dvout   ARPACK utility routine for vector output.
 c     dlarnv  LAPACK routine for generating a random vector.
 c     dgemv   Level 2 BLAS routine for matrix vector multiplication.
@@ -167,7 +167,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   dlarnv, dvout, dcopy, dgemv, second
+      external   dlarnv, dvout, dcopy, dgemv, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -214,7 +214,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mgetv0
 c 
          ierr   = 0
@@ -241,7 +241,7 @@
 c        | the generalized problem when B is possibly (singular).   |
 c        %----------------------------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nopx = nopx + 1
             ipntr(1) = 1
@@ -265,7 +265,7 @@
       if (orth)  go to 40
 c 
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
       end if
 c 
@@ -274,7 +274,7 @@
 c     | Compute B-norm of starting vector.                   |
 c     %------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       first = .TRUE.
       if (bmat .eq. 'G') then
          nbx = nbx + 1
@@ -290,7 +290,7 @@
    20 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -333,7 +333,7 @@
 c     | Compute the B-norm of the orthogonalized starting vector |
 c     %----------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       if (bmat .eq. 'G') then
          nbx = nbx + 1
          call dcopy (n, resid, 1, workd(n+1), 1)
@@ -348,7 +348,7 @@
    40 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -406,7 +406,7 @@
       end if
       ido = 99
 c 
-      call second (t1)
+      call arscnd (t1)
       tgetv0 = tgetv0 + (t1 - t0)
 c 
  9000 continue
diff -urN original/ARPACK/SRC/dnaitr.f ARPACK/SRC/dnaitr.f
--- original/ARPACK/SRC/dnaitr.f	1996-08-27 18:09:00.000000000 +0200
+++ ARPACK/SRC/dnaitr.f	2005-06-30 23:20:38.891791783 +0200
@@ -132,7 +132,7 @@
 c\Routines called:
 c     dgetv0  ARPACK routine to generate the initial vector.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dmout   ARPACK utility routine that prints matrices
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlabad  LAPACK routine that computes machine constants.
@@ -268,7 +268,7 @@
 c     %----------------------%
 c
       external   daxpy, dcopy, dscal, dgemv, dgetv0, dlabad, 
-     &           dvout, dmout, ivout, second
+     &           dvout, dmout, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -319,7 +319,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mnaitr
 c 
 c        %------------------------------%
@@ -430,7 +430,7 @@
 c              %------------------------------------------------%
 c
                info = j - 1
-               call second (t1)
+               call arscnd (t1)
                tnaitr = tnaitr + (t1 - t0)
                ido = 99
                go to 9000
@@ -470,7 +470,7 @@
 c
          step3 = .true.
          nopx  = nopx + 1
-         call second (t2)
+         call arscnd (t2)
          call dcopy (n, v(1,j), 1, workd(ivj), 1)
          ipntr(1) = ivj
          ipntr(2) = irj
@@ -490,7 +490,7 @@
 c        | if step3 = .true.                |
 c        %----------------------------------%
 c
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
  
          step3 = .false.
@@ -506,7 +506,7 @@
 c        |          factorization to length j.   |
 c        %---------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             step4 = .true.
@@ -531,7 +531,7 @@
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -576,11 +576,11 @@
 c
          if (j .gt. 1) h(j,j-1) = betaj
 c
-         call second (t4)
+         call arscnd (t4)
 c 
          orth1 = .true.
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call dcopy (n, resid, 1, workd(irj), 1)
@@ -604,7 +604,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -681,7 +681,7 @@
          call daxpy (j, one, workd(irj), 1, h(1,j), 1)
 c 
          orth2 = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call dcopy (n, resid, 1, workd(irj), 1)
@@ -705,7 +705,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c
@@ -783,7 +783,7 @@
          rstart = .false.
          orth2  = .false.
 c 
-         call second (t5)
+         call arscnd (t5)
          titref = titref + (t5 - t4)
 c 
 c        %------------------------------------%
@@ -792,7 +792,7 @@
 c
          j = j + 1
          if (j .gt. k+np) then
-            call second (t1)
+            call arscnd (t1)
             tnaitr = tnaitr + (t1 - t0)
             ido = 99
             do 110 i = max(1,k), k+np-1
diff -urN original/ARPACK/SRC/dnapps.f ARPACK/SRC/dnapps.f
--- original/ARPACK/SRC/dnapps.f	1998-05-20 16:58:56.000000000 +0200
+++ ARPACK/SRC/dnapps.f	2005-06-30 23:20:38.896791199 +0200
@@ -97,7 +97,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dmout   ARPACK utility routine that prints matrices.
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlabad  LAPACK routine that computes machine constants.
@@ -189,7 +189,7 @@
 c     %----------------------%
 c
       external   daxpy, dcopy, dscal, dlacpy, dlarfg, dlarf,
-     &           dlaset, dlabad, second, dlartg
+     &           dlaset, dlabad, arscnd, dlartg
 c
 c     %--------------------%
 c     | External Functions |
@@ -237,7 +237,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mnapps
       kplusp = kev + np 
 c 
@@ -635,7 +635,7 @@
       end if
 c 
  9000 continue
-      call second (t1)
+      call arscnd (t1)
       tnapps = tnapps + (t1 - t0)
 c 
       return
--- ARPACK/SRC/dnaup2.f.orig	1996-08-16 22:28:09.000000000 +0200
+++ ARPACK/SRC/dnaup2.f	2022-12-15 21:12:42.177193805 +0100
@@ -143,7 +143,7 @@ c     dneigh  ARPACK compute Ritz values
 c     dngets  ARPACK reorder Ritz values and error bounds routine.
 c     dsortc  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dmout   ARPACK utility routine that prints matrices
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlamch  LAPACK routine that determines machine constants.
@@ -234,7 +234,7 @@ c     | External Subroutines |
 c     %----------------------%
 c
       external   dcopy, dgetv0, dnaitr, dnconv, dneigh, dngets, dnapps,
-     &           dvout, ivout, second
+     &           dvout, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -256,7 +256,7 @@ c     %-----------------------%
 c
       if (ido .eq. 0) then
 c 
-         call second (t0)
+         call arscnd (t0)
 c 
          msglvl = mnaup2
 c 
@@ -755,7 +755,7 @@ c        | the first step of the next ca
 c        %---------------------------------------------%
 c
          cnorm = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call dcopy (n, resid, 1, workd(n+1), 1)
@@ -780,7 +780,7 @@ c        | WORKD(1:N) := B*RESID
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -819,7 +819,7 @@ c     %------------%
 c     | Error Exit |
 c     %------------%
 c
-      call second (t1)
+      call arscnd (t1)
       tnaup2 = t1 - t0
 c     
  9000 continue
diff -urN original/ARPACK/SRC/dnaupd.f ARPACK/SRC/dnaupd.f
--- original/ARPACK/SRC/dnaupd.f	2002-08-23 07:29:34.000000000 +0200
+++ ARPACK/SRC/dnaupd.f	2005-06-30 23:20:38.906790031 +0200
@@ -379,7 +379,7 @@
 c     dnaup2  ARPACK routine that implements the Implicitly Restarted
 c             Arnoldi Iteration.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlamch  LAPACK routine that determines machine constants.
 c
@@ -454,7 +454,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   dnaup2, dvout, ivout, second, dstatn
+      external   dnaup2, dvout, ivout, arscnd, dstatn
 c
 c     %--------------------%
 c     | External Functions |
@@ -476,7 +476,7 @@
 c        %-------------------------------%
 c
          call dstatn
-         call second (t0)
+         call arscnd (t0)
          msglvl = mnaupd
 c
 c        %----------------%
@@ -640,7 +640,7 @@
      &               '_naupd: Associated Ritz estimates')
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tnaupd = t1 - t0
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/dnconv.f ARPACK/SRC/dnconv.f
--- original/ARPACK/SRC/dnconv.f	1996-04-20 18:27:05.000000000 +0200
+++ ARPACK/SRC/dnconv.f	2005-06-30 23:20:38.907789915 +0200
@@ -38,7 +38,7 @@
 c     xxxxxx  real
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dlamch  LAPACK routine that determines machine constants.
 c     dlapy2  LAPACK routine to compute sqrt(x**2+y**2) carefully.
 c
@@ -119,7 +119,7 @@
 c     | for some appropriate choice of norm.                        |
 c     %-------------------------------------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
 c
 c     %---------------------------------%
 c     | Get machine dependent constant. |
@@ -134,7 +134,7 @@
          if (bounds(i) .le. tol*temp)   nconv = nconv + 1
    20 continue
 c 
-      call second (t1)
+      call arscnd (t1)
       tnconv = tnconv + (t1 - t0)
 c 
       return
diff -urN original/ARPACK/SRC/dneigh.f ARPACK/SRC/dneigh.f
--- original/ARPACK/SRC/dneigh.f	1996-04-20 18:27:06.000000000 +0200
+++ ARPACK/SRC/dneigh.f	2005-06-30 23:20:38.909789681 +0200
@@ -63,7 +63,7 @@
 c\Routines called:
 c     dlaqrb  ARPACK routine to compute the real Schur form of an
 c             upper Hessenberg matrix and last row of the Schur vectors.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dmout   ARPACK utility routine that prints matrices
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlacpy  LAPACK matrix copy routine.
@@ -144,7 +144,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   dcopy, dlacpy, dlaqrb, dtrevc, dvout, second
+      external   dcopy, dlacpy, dlaqrb, dtrevc, dvout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -170,7 +170,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mneigh
 c 
       if (msglvl .gt. 2) then
@@ -301,7 +301,7 @@
      &              '_neigh: Ritz estimates for the eigenvalues of H')
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tneigh = tneigh + (t1 - t0)
 c
  9000 continue
diff -urN original/ARPACK/SRC/dngets.f ARPACK/SRC/dngets.f
--- original/ARPACK/SRC/dngets.f	1996-04-20 18:27:07.000000000 +0200
+++ ARPACK/SRC/dngets.f	2005-06-30 23:20:38.918788630 +0200
@@ -135,7 +135,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   dcopy, dsortc, second
+      external   dcopy, dsortc, arscnd
 c
 c     %----------------------%
 c     | Intrinsics Functions |
@@ -152,7 +152,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c 
-      call second (t0)
+      call arscnd (t0)
       msglvl = mngets
 c 
 c     %----------------------------------------------------%
@@ -208,7 +208,7 @@
          call dsortc ( 'SR', .true., np, bounds, ritzr, ritzi )
       end if
 c     
-      call second (t1)
+      call arscnd (t1)
       tngets = tngets + (t1 - t0)
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/dsaitr.f ARPACK/SRC/dsaitr.f
--- original/ARPACK/SRC/dsaitr.f	1996-08-28 16:21:43.000000000 +0200
+++ ARPACK/SRC/dsaitr.f	2005-06-30 23:20:38.923788047 +0200
@@ -262,7 +262,7 @@
 c     %----------------------%
 c
       external   daxpy, dcopy, dscal, dgemv, dgetv0, dvout, dmout,
-     &           dlascl, ivout, second
+     &           dlascl, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -300,7 +300,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = msaitr
 c 
 c        %------------------------------%
@@ -420,7 +420,7 @@
 c              %------------------------------------------------%
 c
                info = j - 1
-               call second (t1)
+               call arscnd (t1)
                tsaitr = tsaitr + (t1 - t0)
                ido = 99
                go to 9000
@@ -460,7 +460,7 @@
 c
          step3 = .true.
          nopx  = nopx + 1
-         call second (t2)
+         call arscnd (t2)
          call dcopy (n, v(1,j), 1, workd(ivj), 1)
          ipntr(1) = ivj
          ipntr(2) = irj
@@ -479,7 +479,7 @@
 c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j}.   |
 c        %-----------------------------------%
 c
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
 c 
          step3 = .false.
@@ -500,7 +500,7 @@
 c        %-------------------------------------------%
 c
          if (mode .eq. 2) go to 65
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             step4 = .true.
@@ -524,7 +524,7 @@
 c        %-----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if 
 c
@@ -592,12 +592,12 @@
          else
             h(j,1) = rnorm
          end if
-         call second (t4)
+         call arscnd (t4)
 c 
          orth1 = .true.
          iter  = 0
 c 
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call dcopy (n, resid, 1, workd(irj), 1)
@@ -621,7 +621,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -695,7 +695,7 @@
          h(j,2) = h(j,2) + workd(irj + j - 1)
 c 
          orth2 = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call dcopy (n, resid, 1, workd(irj), 1)
@@ -719,7 +719,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c
@@ -791,7 +791,7 @@
          rstart = .false.
          orth2  = .false.
 c 
-         call second (t5)
+         call arscnd (t5)
          titref = titref + (t5 - t4)
 c 
 c        %----------------------------------------------------------%
@@ -815,7 +815,7 @@
 c
          j = j + 1
          if (j .gt. k+np) then
-            call second (t1)
+            call arscnd (t1)
             tsaitr = tsaitr + (t1 - t0)
             ido = 99
 c
diff -urN original/ARPACK/SRC/dsapps.f ARPACK/SRC/dsapps.f
--- original/ARPACK/SRC/dsapps.f	1998-05-20 16:58:59.000000000 +0200
+++ ARPACK/SRC/dsapps.f	2005-06-30 23:20:38.926787696 +0200
@@ -91,7 +91,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers. 
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlamch  LAPACK routine that determines machine constants.
 c     dlartg  LAPACK Givens rotation construction routine.
@@ -176,7 +176,7 @@
 c     %----------------------%
 c
       external   daxpy, dcopy, dscal, dlacpy, dlartg, dlaset, dvout, 
-     &           ivout, second, dgemv
+     &           ivout, arscnd, dgemv
 c
 c     %--------------------%
 c     | External Functions |
@@ -213,7 +213,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = msapps
 c 
       kplusp = kev + np 
@@ -503,7 +503,7 @@
          end if
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tsapps = tsapps + (t1 - t0)
 c 
  9000 continue 
diff -urN original/ARPACK/SRC/dsaup2.f ARPACK/SRC/dsaup2.f
--- original/ARPACK/SRC/dsaup2.f	1998-05-20 16:59:00.000000000 +0200
+++ ARPACK/SRC/dsaup2.f	2005-06-30 23:20:38.933786879 +0200
@@ -148,7 +148,7 @@
 c     dsgets  ARPACK reorder Ritz values and error bounds routine.
 c     dsortr  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlamch  LAPACK routine that determines machine constants.
 c     dcopy   Level 1 BLAS that copies one vector to another.
@@ -235,7 +235,7 @@
 c     %----------------------%
 c
       external   dcopy, dgetv0, dsaitr, dscal, dsconv, dseigt, dsgets, 
-     &           dsapps, dsortr, dvout, ivout, second, dswap
+     &           dsapps, dsortr, dvout, ivout, arscnd, dswap
 c
 c     %--------------------%
 c     | External Functions |
@@ -262,7 +262,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = msaup2
 c
 c        %---------------------------------%
@@ -770,7 +770,7 @@
 c        %---------------------------------------------%
 c
          cnorm = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call dcopy (n, resid, 1, workd(n+1), 1)
@@ -795,7 +795,7 @@
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -837,7 +837,7 @@
 c     | Error exit |
 c     %------------%
 c
-      call second (t1)
+      call arscnd (t1)
       tsaup2 = t1 - t0
 c 
  9000 continue
--- ARPACK/SRC/dsaupd.f.orig	1996-08-27 07:28:56.000000000 +0200
+++ ARPACK/SRC/dsaupd.f	2022-12-15 21:14:09.256722054 +0100
@@ -380,7 +380,7 @@ c             Arnoldi Iteration.
 c     dstats  ARPACK routine that initialize timing and other statistics
 c             variables.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dvout   ARPACK utility routine that prints vectors.
 c     dlamch  LAPACK routine that determines machine constants.
 c
@@ -456,7 +456,7 @@ c     %----------------------%
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   dsaup2,  dvout, ivout, second, dstats
+      external   dsaup2,  dvout, ivout, arscnd, dstats
 c
 c     %--------------------%
 c     | External Functions |
@@ -478,7 +478,7 @@ c        | & message level for debugging
 c        %-------------------------------%
 c
          call dstats
-         call second (t0)
+         call arscnd (t0)
          msglvl = msaupd
 c
          ierr   = 0
@@ -637,7 +637,7 @@ c
      &               '_saupd: corresponding error bounds')
       end if 
 c
-      call second (t1)
+      call arscnd (t1)
       tsaupd = t1 - t0
 c 
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/dsconv.f ARPACK/SRC/dsconv.f
--- original/ARPACK/SRC/dsconv.f	1996-04-20 18:27:11.000000000 +0200
+++ ARPACK/SRC/dsconv.f	2005-06-30 23:20:38.939786179 +0200
@@ -34,7 +34,7 @@
 c\BeginLib
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dlamch  LAPACK routine that determines machine constants. 
 c
 c\Author
@@ -106,7 +106,7 @@
 c     | Executable Statements |
 c     %-----------------------%
 c
-      call second (t0)
+      call arscnd (t0)
 c
       eps23 = dlamch('Epsilon-Machine') 
       eps23 = eps23**(2.0D+0 / 3.0D+0)
@@ -126,7 +126,7 @@
 c
    10 continue
 c 
-      call second (t1)
+      call arscnd (t1)
       tsconv = tsconv + (t1 - t0)
 c 
       return
diff -urN original/ARPACK/SRC/dseigt.f ARPACK/SRC/dseigt.f
--- original/ARPACK/SRC/dseigt.f	1996-08-27 07:29:04.000000000 +0200
+++ ARPACK/SRC/dseigt.f	2005-06-30 23:20:38.940786062 +0200
@@ -59,7 +59,7 @@
 c     dstqrb  ARPACK routine that computes the eigenvalues and the
 c             last components of the eigenvectors of a symmetric
 c             and tridiagonal matrix.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dvout   ARPACK utility routine that prints vectors.
 c     dcopy   Level 1 BLAS that copies one vector to another.
 c
@@ -127,7 +127,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   dcopy, dstqrb, dvout, second
+      external   dcopy, dstqrb, dvout, arscnd
 c
 c     %-----------------------%
 c     | Executable Statements |
@@ -138,7 +138,7 @@
 c     | & message level for debugging |
 c     %-------------------------------% 
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mseigt
 c
       if (msglvl .gt. 0) then
@@ -168,7 +168,7 @@
          bounds(k) = rnorm*abs(bounds(k))
    30 continue
 c 
-      call second (t1)
+      call arscnd (t1)
       tseigt = tseigt + (t1 - t0)
 c
  9000 continue
diff -urN original/ARPACK/SRC/dsgets.f ARPACK/SRC/dsgets.f
--- original/ARPACK/SRC/dsgets.f	1996-04-20 18:27:14.000000000 +0200
+++ ARPACK/SRC/dsgets.f	2005-06-30 23:20:38.949785011 +0200
@@ -65,7 +65,7 @@
 c\Routines called:
 c     dsortr  ARPACK utility sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dvout   ARPACK utility routine that prints vectors.
 c     dcopy   Level 1 BLAS that copies one vector to another.
 c     dswap   Level 1 BLAS that swaps the contents of two vectors.
@@ -131,7 +131,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   dswap, dcopy, dsortr, second
+      external   dswap, dcopy, dsortr, arscnd
 c
 c     %---------------------%
 c     | Intrinsic Functions |
@@ -148,7 +148,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = msgets
 c 
       if (which .eq. 'BE') then
@@ -198,7 +198,7 @@
          call dcopy (np, ritz, 1, shifts, 1)
       end if
 c 
-      call second (t1)
+      call arscnd (t1)
       tsgets = tsgets + (t1 - t0)
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/sgetv0.f ARPACK/SRC/sgetv0.f
--- original/ARPACK/SRC/sgetv0.f	1999-04-08 00:42:46.000000000 +0200
+++ ARPACK/SRC/sgetv0.f	2005-06-30 23:20:38.960783727 +0200
@@ -93,7 +93,7 @@
 c     TR95-13, Department of Computational and Applied Mathematics.
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     svout   ARPACK utility routine for vector output.
 c     slarnv  LAPACK routine for generating a random vector.
 c     sgemv   Level 2 BLAS routine for matrix vector multiplication.
@@ -167,7 +167,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   slarnv, svout, scopy, sgemv, second
+      external   slarnv, svout, scopy, sgemv, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -214,7 +214,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mgetv0
 c 
          ierr   = 0
@@ -241,7 +241,7 @@
 c        | the generalized problem when B is possibly (singular).   |
 c        %----------------------------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nopx = nopx + 1
             ipntr(1) = 1
@@ -265,7 +265,7 @@
       if (orth)  go to 40
 c 
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
       end if
 c 
@@ -274,7 +274,7 @@
 c     | Compute B-norm of starting vector.                   |
 c     %------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       first = .TRUE.
       if (bmat .eq. 'G') then
          nbx = nbx + 1
@@ -290,7 +290,7 @@
    20 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -333,7 +333,7 @@
 c     | Compute the B-norm of the orthogonalized starting vector |
 c     %----------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       if (bmat .eq. 'G') then
          nbx = nbx + 1
          call scopy (n, resid, 1, workd(n+1), 1)
@@ -348,7 +348,7 @@
    40 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -406,7 +406,7 @@
       end if
       ido = 99
 c 
-      call second (t1)
+      call arscnd (t1)
       tgetv0 = tgetv0 + (t1 - t0)
 c 
  9000 continue
diff -urN original/ARPACK/SRC/snaitr.f ARPACK/SRC/snaitr.f
--- original/ARPACK/SRC/snaitr.f	1996-08-27 18:09:01.000000000 +0200
+++ ARPACK/SRC/snaitr.f	2005-06-30 23:20:38.969782676 +0200
@@ -132,7 +132,7 @@
 c\Routines called:
 c     sgetv0  ARPACK routine to generate the initial vector.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     smout   ARPACK utility routine that prints matrices
 c     svout   ARPACK utility routine that prints vectors.
 c     slabad  LAPACK routine that computes machine constants.
@@ -268,7 +268,7 @@
 c     %----------------------%
 c
       external   saxpy, scopy, sscal, sgemv, sgetv0, slabad, 
-     &           svout, smout, ivout, second
+     &           svout, smout, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -319,7 +319,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mnaitr
 c 
 c        %------------------------------%
@@ -430,7 +430,7 @@
 c              %------------------------------------------------%
 c
                info = j - 1
-               call second (t1)
+               call arscnd (t1)
                tnaitr = tnaitr + (t1 - t0)
                ido = 99
                go to 9000
@@ -470,7 +470,7 @@
 c
          step3 = .true.
          nopx  = nopx + 1
-         call second (t2)
+         call arscnd (t2)
          call scopy (n, v(1,j), 1, workd(ivj), 1)
          ipntr(1) = ivj
          ipntr(2) = irj
@@ -490,7 +490,7 @@
 c        | if step3 = .true.                |
 c        %----------------------------------%
 c
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
  
          step3 = .false.
@@ -506,7 +506,7 @@
 c        |          factorization to length j.   |
 c        %---------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             step4 = .true.
@@ -531,7 +531,7 @@
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -576,11 +576,11 @@
 c
          if (j .gt. 1) h(j,j-1) = betaj
 c
-         call second (t4)
+         call arscnd (t4)
 c 
          orth1 = .true.
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call scopy (n, resid, 1, workd(irj), 1)
@@ -604,7 +604,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -681,7 +681,7 @@
          call saxpy (j, one, workd(irj), 1, h(1,j), 1)
 c 
          orth2 = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call scopy (n, resid, 1, workd(irj), 1)
@@ -705,7 +705,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c
@@ -783,7 +783,7 @@
          rstart = .false.
          orth2  = .false.
 c 
-         call second (t5)
+         call arscnd (t5)
          titref = titref + (t5 - t4)
 c 
 c        %------------------------------------%
@@ -792,7 +792,7 @@
 c
          j = j + 1
          if (j .gt. k+np) then
-            call second (t1)
+            call arscnd (t1)
             tnaitr = tnaitr + (t1 - t0)
             ido = 99
             do 110 i = max(1,k), k+np-1
diff -urN original/ARPACK/SRC/snapps.f ARPACK/SRC/snapps.f
--- original/ARPACK/SRC/snapps.f	1998-05-20 16:58:55.000000000 +0200
+++ ARPACK/SRC/snapps.f	2005-06-30 23:20:38.973782209 +0200
@@ -97,7 +97,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     smout   ARPACK utility routine that prints matrices.
 c     svout   ARPACK utility routine that prints vectors.
 c     slabad  LAPACK routine that computes machine constants.
@@ -189,7 +189,7 @@
 c     %----------------------%
 c
       external   saxpy, scopy, sscal, slacpy, slarfg, slarf,
-     &           slaset, slabad, second, slartg
+     &           slaset, slabad, arscnd, slartg
 c
 c     %--------------------%
 c     | External Functions |
@@ -237,7 +237,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mnapps
       kplusp = kev + np 
 c 
@@ -635,7 +635,7 @@
       end if
 c 
  9000 continue
-      call second (t1)
+      call arscnd (t1)
       tnapps = tnapps + (t1 - t0)
 c 
       return
--- ARPACK/SRC/snaup2.f.orig	1996-08-16 22:28:10.000000000 +0200
+++ ARPACK/SRC/snaup2.f	2022-12-15 21:16:41.849228723 +0100
@@ -143,7 +143,7 @@ c     sneigh  ARPACK compute Ritz values
 c     sngets  ARPACK reorder Ritz values and error bounds routine.
 c     ssortc  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     smout   ARPACK utility routine that prints matrices
 c     svout   ARPACK utility routine that prints vectors.
 c     slamch  LAPACK routine that determines machine constants.
@@ -234,7 +234,7 @@ c     | External Subroutines |
 c     %----------------------%
 c
       external   scopy, sgetv0, snaitr, snconv, sneigh, sngets, snapps,
-     &           svout, ivout, second
+     &           svout, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -256,7 +256,7 @@ c     %-----------------------%
 c
       if (ido .eq. 0) then
 c 
-         call second (t0)
+         call arscnd (t0)
 c 
          msglvl = mnaup2
 c 
@@ -755,7 +755,7 @@ c        | the first step of the next ca
 c        %---------------------------------------------%
 c
          cnorm = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call scopy (n, resid, 1, workd(n+1), 1)
@@ -780,7 +780,7 @@ c        | WORKD(1:N) := B*RESID
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -819,7 +819,7 @@ c     %------------%
 c     | Error Exit |
 c     %------------%
 c
-      call second (t1)
+      call arscnd (t1)
       tnaup2 = t1 - t0
 c     
  9000 continue
diff -urN original/ARPACK/SRC/snaupd.f ARPACK/SRC/snaupd.f
--- original/ARPACK/SRC/snaupd.f	2002-08-23 07:29:34.000000000 +0200
+++ ARPACK/SRC/snaupd.f	2005-06-30 23:20:38.984780925 +0200
@@ -379,7 +379,7 @@
 c     snaup2  ARPACK routine that implements the Implicitly Restarted
 c             Arnoldi Iteration.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     svout   ARPACK utility routine that prints vectors.
 c     slamch  LAPACK routine that determines machine constants.
 c
@@ -454,7 +454,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   snaup2, svout, ivout, second, sstatn
+      external   snaup2, svout, ivout, arscnd, sstatn
 c
 c     %--------------------%
 c     | External Functions |
@@ -476,7 +476,7 @@
 c        %-------------------------------%
 c
          call sstatn
-         call second (t0)
+         call arscnd (t0)
          msglvl = mnaupd
 c
 c        %----------------%
@@ -640,7 +640,7 @@
      &               '_naupd: Associated Ritz estimates')
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tnaupd = t1 - t0
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/snconv.f ARPACK/SRC/snconv.f
--- original/ARPACK/SRC/snconv.f	1996-04-20 18:27:21.000000000 +0200
+++ ARPACK/SRC/snconv.f	2005-06-30 23:20:38.985780808 +0200
@@ -38,7 +38,7 @@
 c     xxxxxx  real
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     slamch  LAPACK routine that determines machine constants.
 c     slapy2  LAPACK routine to compute sqrt(x**2+y**2) carefully.
 c
@@ -119,7 +119,7 @@
 c     | for some appropriate choice of norm.                        |
 c     %-------------------------------------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
 c
 c     %---------------------------------%
 c     | Get machine dependent constant. |
@@ -134,7 +134,7 @@
          if (bounds(i) .le. tol*temp)   nconv = nconv + 1
    20 continue
 c 
-      call second (t1)
+      call arscnd (t1)
       tnconv = tnconv + (t1 - t0)
 c 
       return
diff -urN original/ARPACK/SRC/sneigh.f ARPACK/SRC/sneigh.f
--- original/ARPACK/SRC/sneigh.f	1996-04-20 18:27:22.000000000 +0200
+++ ARPACK/SRC/sneigh.f	2005-06-30 23:20:38.987780575 +0200
@@ -63,7 +63,7 @@
 c\Routines called:
 c     slaqrb  ARPACK routine to compute the real Schur form of an
 c             upper Hessenberg matrix and last row of the Schur vectors.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     smout   ARPACK utility routine that prints matrices
 c     svout   ARPACK utility routine that prints vectors.
 c     slacpy  LAPACK matrix copy routine.
@@ -144,7 +144,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   scopy, slacpy, slaqrb, strevc, svout, second
+      external   scopy, slacpy, slaqrb, strevc, svout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -170,7 +170,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mneigh
 c 
       if (msglvl .gt. 2) then
@@ -301,7 +301,7 @@
      &              '_neigh: Ritz estimates for the eigenvalues of H')
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tneigh = tneigh + (t1 - t0)
 c
  9000 continue
diff -urN original/ARPACK/SRC/sngets.f ARPACK/SRC/sngets.f
--- original/ARPACK/SRC/sngets.f	1996-04-20 18:27:24.000000000 +0200
+++ ARPACK/SRC/sngets.f	2005-06-30 23:20:38.995779641 +0200
@@ -135,7 +135,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   scopy, ssortc, second
+      external   scopy, ssortc, arscnd
 c
 c     %----------------------%
 c     | Intrinsics Functions |
@@ -152,7 +152,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c 
-      call second (t0)
+      call arscnd (t0)
       msglvl = mngets
 c 
 c     %----------------------------------------------------%
@@ -208,7 +208,7 @@
          call ssortc ( 'SR', .true., np, bounds, ritzr, ritzi )
       end if
 c     
-      call second (t1)
+      call arscnd (t1)
       tngets = tngets + (t1 - t0)
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/ssaitr.f ARPACK/SRC/ssaitr.f
--- original/ARPACK/SRC/ssaitr.f	1996-08-28 16:21:43.000000000 +0200
+++ ARPACK/SRC/ssaitr.f	2005-06-30 23:20:39.001778940 +0200
@@ -262,7 +262,7 @@
 c     %----------------------%
 c
       external   saxpy, scopy, sscal, sgemv, sgetv0, svout, smout,
-     &           slascl, ivout, second
+     &           slascl, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -300,7 +300,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = msaitr
 c 
 c        %------------------------------%
@@ -420,7 +420,7 @@
 c              %------------------------------------------------%
 c
                info = j - 1
-               call second (t1)
+               call arscnd (t1)
                tsaitr = tsaitr + (t1 - t0)
                ido = 99
                go to 9000
@@ -460,7 +460,7 @@
 c
          step3 = .true.
          nopx  = nopx + 1
-         call second (t2)
+         call arscnd (t2)
          call scopy (n, v(1,j), 1, workd(ivj), 1)
          ipntr(1) = ivj
          ipntr(2) = irj
@@ -479,7 +479,7 @@
 c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j}.   |
 c        %-----------------------------------%
 c
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
 c 
          step3 = .false.
@@ -500,7 +500,7 @@
 c        %-------------------------------------------%
 c
          if (mode .eq. 2) go to 65
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             step4 = .true.
@@ -524,7 +524,7 @@
 c        %-----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if 
 c
@@ -592,12 +592,12 @@
          else
             h(j,1) = rnorm
          end if
-         call second (t4)
+         call arscnd (t4)
 c 
          orth1 = .true.
          iter  = 0
 c 
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call scopy (n, resid, 1, workd(irj), 1)
@@ -621,7 +621,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -695,7 +695,7 @@
          h(j,2) = h(j,2) + workd(irj + j - 1)
 c 
          orth2 = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call scopy (n, resid, 1, workd(irj), 1)
@@ -719,7 +719,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c
@@ -791,7 +791,7 @@
          rstart = .false.
          orth2  = .false.
 c 
-         call second (t5)
+         call arscnd (t5)
          titref = titref + (t5 - t4)
 c 
 c        %----------------------------------------------------------%
@@ -815,7 +815,7 @@
 c
          j = j + 1
          if (j .gt. k+np) then
-            call second (t1)
+            call arscnd (t1)
             tsaitr = tsaitr + (t1 - t0)
             ido = 99
 c
diff -urN original/ARPACK/SRC/ssapps.f ARPACK/SRC/ssapps.f
--- original/ARPACK/SRC/ssapps.f	1998-05-20 16:58:58.000000000 +0200
+++ ARPACK/SRC/ssapps.f	2005-06-30 23:20:39.004778590 +0200
@@ -91,7 +91,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers. 
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     svout   ARPACK utility routine that prints vectors.
 c     slamch  LAPACK routine that determines machine constants.
 c     slartg  LAPACK Givens rotation construction routine.
@@ -176,7 +176,7 @@
 c     %----------------------%
 c
       external   saxpy, scopy, sscal, slacpy, slartg, slaset, svout, 
-     &           ivout, second, sgemv
+     &           ivout, arscnd, sgemv
 c
 c     %--------------------%
 c     | External Functions |
@@ -213,7 +213,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = msapps
 c 
       kplusp = kev + np 
@@ -503,7 +503,7 @@
          end if
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tsapps = tsapps + (t1 - t0)
 c 
  9000 continue 
diff -urN original/ARPACK/SRC/ssaup2.f ARPACK/SRC/ssaup2.f
--- original/ARPACK/SRC/ssaup2.f	1998-05-20 16:58:59.000000000 +0200
+++ ARPACK/SRC/ssaup2.f	2005-06-30 23:20:39.010777890 +0200
@@ -148,7 +148,7 @@
 c     ssgets  ARPACK reorder Ritz values and error bounds routine.
 c     ssortr  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     svout   ARPACK utility routine that prints vectors.
 c     slamch  LAPACK routine that determines machine constants.
 c     scopy   Level 1 BLAS that copies one vector to another.
@@ -235,7 +235,7 @@
 c     %----------------------%
 c
       external   scopy, sgetv0, ssaitr, sscal, ssconv, sseigt, ssgets, 
-     &           ssapps, ssortr, svout, ivout, second, sswap
+     &           ssapps, ssortr, svout, ivout, arscnd, sswap
 c
 c     %--------------------%
 c     | External Functions |
@@ -262,7 +262,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = msaup2
 c
 c        %---------------------------------%
@@ -770,7 +770,7 @@
 c        %---------------------------------------------%
 c
          cnorm = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call scopy (n, resid, 1, workd(n+1), 1)
@@ -795,7 +795,7 @@
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -837,7 +837,7 @@
 c     | Error exit |
 c     %------------%
 c
-      call second (t1)
+      call arscnd (t1)
       tsaup2 = t1 - t0
 c 
  9000 continue
diff -urN original/ARPACK/SRC/ssaupd.f ARPACK/SRC/ssaupd.f
--- original/ARPACK/SRC/ssaupd.f	2001-04-11 00:52:40.000000000 +0200
+++ ARPACK/SRC/ssaupd.f	2005-06-30 23:20:39.015777306 +0200
@@ -380,7 +380,7 @@
 c     sstats  ARPACK routine that initialize timing and other statistics
 c             variables.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     svout   ARPACK utility routine that prints vectors.
 c     slamch  LAPACK routine that determines machine constants.
 c
@@ -456,7 +456,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   ssaup2,  svout, ivout, second, sstats
+      external   ssaup2,  svout, ivout, arscnd, sstats
 c
 c     %--------------------%
 c     | External Functions |
@@ -478,7 +478,7 @@
 c        %-------------------------------%
 c
          call sstats
-         call second (t0)
+         call arscnd (t0)
          msglvl = msaupd
 c
          ierr   = 0
@@ -638,7 +638,7 @@
      &               '_saupd: corresponding error bounds')
       end if 
 c
-      call second (t1)
+      call arscnd (t1)
       tsaupd = t1 - t0
 c 
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/ssconv.f ARPACK/SRC/ssconv.f
--- original/ARPACK/SRC/ssconv.f	1996-04-20 18:27:27.000000000 +0200
+++ ARPACK/SRC/ssconv.f	2005-06-30 23:20:39.016777190 +0200
@@ -34,7 +34,7 @@
 c\BeginLib
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     slamch  LAPACK routine that determines machine constants. 
 c
 c\Author
@@ -106,7 +106,7 @@
 c     | Executable Statements |
 c     %-----------------------%
 c
-      call second (t0)
+      call arscnd (t0)
 c
       eps23 = slamch('Epsilon-Machine') 
       eps23 = eps23**(2.0E+0 / 3.0E+0)
@@ -126,7 +126,7 @@
 c
    10 continue
 c 
-      call second (t1)
+      call arscnd (t1)
       tsconv = tsconv + (t1 - t0)
 c 
       return
diff -urN original/ARPACK/SRC/sseigt.f ARPACK/SRC/sseigt.f
--- original/ARPACK/SRC/sseigt.f	1996-08-27 07:29:04.000000000 +0200
+++ ARPACK/SRC/sseigt.f	2005-06-30 23:20:39.017777073 +0200
@@ -59,7 +59,7 @@
 c     sstqrb  ARPACK routine that computes the eigenvalues and the
 c             last components of the eigenvectors of a symmetric
 c             and tridiagonal matrix.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     svout   ARPACK utility routine that prints vectors.
 c     scopy   Level 1 BLAS that copies one vector to another.
 c
@@ -127,7 +127,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   scopy, sstqrb, svout, second
+      external   scopy, sstqrb, svout, arscnd
 c
 c     %-----------------------%
 c     | Executable Statements |
@@ -138,7 +138,7 @@
 c     | & message level for debugging |
 c     %-------------------------------% 
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mseigt
 c
       if (msglvl .gt. 0) then
@@ -168,7 +168,7 @@
          bounds(k) = rnorm*abs(bounds(k))
    30 continue
 c 
-      call second (t1)
+      call arscnd (t1)
       tseigt = tseigt + (t1 - t0)
 c
  9000 continue
diff -urN original/ARPACK/SRC/ssgets.f ARPACK/SRC/ssgets.f
--- original/ARPACK/SRC/ssgets.f	1996-04-20 18:27:30.000000000 +0200
+++ ARPACK/SRC/ssgets.f	2005-06-30 23:20:39.026776022 +0200
@@ -65,7 +65,7 @@
 c\Routines called:
 c     ssortr  ARPACK utility sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     svout   ARPACK utility routine that prints vectors.
 c     scopy   Level 1 BLAS that copies one vector to another.
 c     sswap   Level 1 BLAS that swaps the contents of two vectors.
@@ -131,7 +131,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   sswap, scopy, ssortr, second
+      external   sswap, scopy, ssortr, arscnd
 c
 c     %---------------------%
 c     | Intrinsic Functions |
@@ -148,7 +148,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = msgets
 c 
       if (which .eq. 'BE') then
@@ -198,7 +198,7 @@
          call scopy (np, ritz, 1, shifts, 1)
       end if
 c 
-      call second (t1)
+      call arscnd (t1)
       tsgets = tsgets + (t1 - t0)
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/zgetv0.f ARPACK/SRC/zgetv0.f
--- original/ARPACK/SRC/zgetv0.f	1999-04-08 00:46:45.000000000 +0200
+++ ARPACK/SRC/zgetv0.f	2005-06-30 23:20:39.037774738 +0200
@@ -89,7 +89,7 @@
 c     pp 357-385.
 c
 c\Routines called:
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     zvout   ARPACK utility routine that prints vectors.
 c     zlarnv  LAPACK routine for generating a random vector. 
 c     zgemv   Level 2 BLAS routine for matrix vector multiplication.
@@ -168,7 +168,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   zcopy, zgemv, zlarnv, zvout, second
+      external   zcopy, zgemv, zlarnv, zvout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -211,7 +211,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mgetv0
 c 
          ierr   = 0
@@ -238,7 +238,7 @@
 c        | the generalized problem when B is possibly (singular).   |
 c        %----------------------------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nopx = nopx + 1
             ipntr(1) = 1
@@ -261,7 +261,7 @@
 c
       if (orth)  go to 40
 c 
-      call second (t3)
+      call arscnd (t3)
       tmvopx = tmvopx + (t3 - t2)
 c 
 c     %------------------------------------------------------%
@@ -269,7 +269,7 @@
 c     | Compute B-norm of starting vector.                   |
 c     %------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       first = .TRUE.
       if (bmat .eq. 'G') then
          nbx = nbx + 1
@@ -285,7 +285,7 @@
    20 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -328,7 +328,7 @@
 c     | Compute the B-norm of the orthogonalized starting vector |
 c     %----------------------------------------------------------%
 c
-      call second (t2)
+      call arscnd (t2)
       if (bmat .eq. 'G') then
          nbx = nbx + 1
          call zcopy (n, resid, 1, workd(n+1), 1)
@@ -343,7 +343,7 @@
    40 continue
 c
       if (bmat .eq. 'G') then
-         call second (t3)
+         call arscnd (t3)
          tmvbx = tmvbx + (t3 - t2)
       end if
 c 
@@ -401,7 +401,7 @@
       end if
       ido = 99
 c 
-      call second (t1)
+      call arscnd (t1)
       tgetv0 = tgetv0 + (t1 - t0)
 c 
  9000 continue
diff -urN original/ARPACK/SRC/znaitr.f ARPACK/SRC/znaitr.f
--- original/ARPACK/SRC/znaitr.f	1996-08-27 18:12:57.000000000 +0200
+++ ARPACK/SRC/znaitr.f	2005-06-30 23:20:39.042774155 +0200
@@ -131,7 +131,7 @@
 c\Routines called:
 c     zgetv0  ARPACK routine to generate the initial vector.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     zmout   ARPACK utility routine that prints matrices
 c     zvout   ARPACK utility routine that prints vectors.
 c     zlanhs  LAPACK routine that computes various norms of a matrix.
@@ -273,7 +273,7 @@
 c     %----------------------%
 c
       external   zaxpy, zcopy, zscal, zdscal, zgemv, zgetv0, 
-     &           dlabad, zvout, zmout, ivout, second
+     &           dlabad, zvout, zmout, ivout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -326,7 +326,7 @@
 c        | & message level for debugging |
 c        %-------------------------------%
 c
-         call second (t0)
+         call arscnd (t0)
          msglvl = mcaitr
 c 
 c        %------------------------------%
@@ -437,7 +437,7 @@
 c              %------------------------------------------------%
 c
                info = j - 1
-               call second (t1)
+               call arscnd (t1)
                tcaitr = tcaitr + (t1 - t0)
                ido = 99
                go to 9000
@@ -477,7 +477,7 @@
 c
          step3 = .true.
          nopx  = nopx + 1
-         call second (t2)
+         call arscnd (t2)
          call zcopy (n, v(1,j), 1, workd(ivj), 1)
          ipntr(1) = ivj
          ipntr(2) = irj
@@ -497,7 +497,7 @@
 c        | if step3 = .true.                |
 c        %----------------------------------%
 c
-         call second (t3)
+         call arscnd (t3)
          tmvopx = tmvopx + (t3 - t2)
  
          step3 = .false.
@@ -513,7 +513,7 @@
 c        |          factorization to length j.   |
 c        %---------------------------------------%
 c
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             step4 = .true.
@@ -538,7 +538,7 @@
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -583,11 +583,11 @@
 c
          if (j .gt. 1) h(j,j-1) = dcmplx(betaj, rzero)
 c
-         call second (t4)
+         call arscnd (t4)
 c 
          orth1 = .true.
 c 
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call zcopy (n, resid, 1, workd(irj), 1)
@@ -611,7 +611,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -689,7 +689,7 @@
          call zaxpy (j, one, workd(irj), 1, h(1,j), 1)
 c 
          orth2 = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call zcopy (n, resid, 1, workd(irj), 1)
@@ -713,7 +713,7 @@
 c        %---------------------------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if 
 c
@@ -791,7 +791,7 @@
          rstart = .false.
          orth2  = .false.
 c 
-         call second (t5)
+         call arscnd (t5)
          titref = titref + (t5 - t4)
 c 
 c        %------------------------------------%
@@ -800,7 +800,7 @@
 c
          j = j + 1
          if (j .gt. k+np) then
-            call second (t1)
+            call arscnd (t1)
             tcaitr = tcaitr + (t1 - t0)
             ido = 99
             do 110 i = max(1,k), k+np-1
diff -urN original/ARPACK/SRC/znapps.f ARPACK/SRC/znapps.f
--- original/ARPACK/SRC/znapps.f	1998-05-20 16:45:03.000000000 +0200
+++ ARPACK/SRC/znapps.f	2005-06-30 23:20:39.045773805 +0200
@@ -92,7 +92,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     zmout   ARPACK utility routine that prints matrices
 c     zvout   ARPACK utility routine that prints vectors.
 c     zlacpy  LAPACK matrix copy routine.
@@ -184,7 +184,7 @@
 c     %----------------------%
 c
       external   zaxpy, zcopy, zgemv, zscal, zlacpy, zlartg, 
-     &           zvout, zlaset, dlabad, zmout, second, ivout
+     &           zvout, zlaset, dlabad, zmout, arscnd, ivout
 c
 c     %--------------------%
 c     | External Functions |
@@ -240,7 +240,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mcapps
 c 
       kplusp = kev + np 
@@ -495,7 +495,7 @@
       end if
 c
  9000 continue
-      call second (t1)
+      call arscnd (t1)
       tcapps = tcapps + (t1 - t0)
 c 
       return
--- ARPACK/SRC/znaup2.f.orig	1996-09-24 22:40:01.000000000 +0200
+++ ARPACK/SRC/znaup2.f	2022-12-15 21:20:55.564520897 +0100
@@ -135,7 +135,7 @@ c     zneigh  ARPACK compute Ritz values
 c     zngets  ARPACK reorder Ritz values and error bounds routine.
 c     zsortc  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     zmout   ARPACK utility routine that prints matrices
 c     zvout   ARPACK utility routine that prints vectors.
 c     dvout   ARPACK utility routine that prints vectors.
@@ -239,7 +239,7 @@ c     | External Subroutines |
 c     %----------------------%
 c
       external   zcopy, zgetv0, znaitr, zneigh, zngets, znapps,
-     &           zsortc, zswap, zmout, zvout, ivout, second
+     &           zsortc, zswap, zmout, zvout, ivout, arscnd
 c
 c     %--------------------%
 c     | External functions |
@@ -263,7 +263,7 @@ c     %-----------------------%
 c
       if (ido .eq. 0) then
 c 
-         call second (t0)
+         call arscnd (t0)
 c 
          msglvl = mcaup2
 c 
@@ -723,7 +723,7 @@ c        | the first step of the next ca
 c        %---------------------------------------------%
 c
          cnorm = .true.
-         call second (t2)
+         call arscnd (t2)
          if (bmat .eq. 'G') then
             nbx = nbx + 1
             call zcopy (n, resid, 1, workd(n+1), 1)
@@ -748,7 +748,7 @@ c        | WORKD(1:N) := B*RESID
 c        %----------------------------------%
 c
          if (bmat .eq. 'G') then
-            call second (t3)
+            call arscnd (t3)
             tmvbx = tmvbx + (t3 - t2)
          end if
 c 
@@ -787,7 +787,7 @@ c     %------------%
 c     | Error Exit |
 c     %------------%
 c
-      call second (t1)
+      call arscnd (t1)
       tcaup2 = t1 - t0
 c     
  9000 continue
diff -urN original/ARPACK/SRC/znaupd.f ARPACK/SRC/znaupd.f
--- original/ARPACK/SRC/znaupd.f	2002-08-15 07:50:57.000000000 +0200
+++ ARPACK/SRC/znaupd.f	2005-06-30 23:20:39.055772637 +0200
@@ -356,7 +356,7 @@
 c     zstatn  ARPACK routine that initializes the timing variables.
 c     ivout   ARPACK utility routine that prints integers.
 c     zvout   ARPACK utility routine that prints vectors.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     dlamch  LAPACK routine that determines machine constants.
 c
 c\Author
@@ -429,7 +429,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   znaup2, zvout, ivout, second, zstatn
+      external   znaup2, zvout, ivout, arscnd, zstatn
 c
 c     %--------------------%
 c     | External Functions |
@@ -451,7 +451,7 @@
 c        %-------------------------------%
 c
          call zstatn
-         call second (t0)
+         call arscnd (t0)
          msglvl = mcaupd
 c
 c        %----------------%
@@ -611,7 +611,7 @@
      &               '_naupd: Associated Ritz estimates')
       end if
 c
-      call second (t1)
+      call arscnd (t1)
       tcaupd = t1 - t0
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/SRC/zneigh.f ARPACK/SRC/zneigh.f
--- original/ARPACK/SRC/zneigh.f	1996-04-20 19:16:01.000000000 +0200
+++ ARPACK/SRC/zneigh.f	2005-06-30 23:20:39.056772521 +0200
@@ -64,7 +64,7 @@
 c
 c\Routines called:
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     zmout   ARPACK utility routine that prints matrices
 c     zvout   ARPACK utility routine that prints vectors.
 c     dvout   ARPACK utility routine that prints vectors.
@@ -152,7 +152,7 @@
 c     %----------------------%
 c
       external   zlacpy, zlahqr, ztrevc, zcopy, 
-     &           zdscal, zmout, zvout, second
+     &           zdscal, zmout, zvout, arscnd
 c
 c     %--------------------%
 c     | External Functions |
@@ -171,7 +171,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c
-      call second (t0)
+      call arscnd (t0)
       msglvl = mceigh
 c 
       if (msglvl .gt. 2) then
@@ -244,7 +244,7 @@
      &              '_neigh: Ritz estimates for the eigenvalues of H')
       end if
 c
-      call second(t1)
+      call arscnd(t1)
       tceigh = tceigh + (t1 - t0)
 c
  9000 continue
diff -urN original/ARPACK/SRC/zngets.f ARPACK/SRC/zngets.f
--- original/ARPACK/SRC/zngets.f	1996-04-20 19:16:02.000000000 +0200
+++ ARPACK/SRC/zngets.f	2005-06-30 23:20:39.063771704 +0200
@@ -63,7 +63,7 @@
 c\Routines called:
 c     zsortc  ARPACK sorting routine.
 c     ivout   ARPACK utility routine that prints integers.
-c     second  ARPACK utility routine for timing.
+c     arscnd  ARPACK utility routine for timing.
 c     zvout   ARPACK utility routine that prints vectors.
 c
 c\Author
@@ -126,7 +126,7 @@
 c     | External Subroutines |
 c     %----------------------%
 c
-      external   zvout,  zsortc, second
+      external   zvout,  zsortc, arscnd
 c
 c     %-----------------------%
 c     | Executable Statements |
@@ -137,7 +137,7 @@
 c     | & message level for debugging |
 c     %-------------------------------%
 c 
-      call second (t0)
+      call arscnd (t0)
       msglvl = mcgets
 c 
       call zsortc (which, .true., kev+np, ritz, bounds)
@@ -157,7 +157,7 @@
 c
       end if
 c     
-      call second (t1)
+      call arscnd (t1)
       tcgets = tcgets + (t1 - t0)
 c
       if (msglvl .gt. 0) then
diff -urN original/ARPACK/UTIL/Makefile ARPACK/UTIL/Makefile
--- original/ARPACK/UTIL/Makefile	1996-09-09 16:30:48.000000000 +0200
+++ ARPACK/UTIL/Makefile	2005-06-30 23:23:46.746224656 +0200
@@ -31,7 +31,7 @@
 #  ../ARmake.inc and is created at the next higher directory level.
 #
 
-OBJS  = icnteq.o icopy.o iset.o iswap.o ivout.o second.o
+OBJS  = icnteq.o icopy.o iset.o iswap.o ivout.o arscnd.o
 
 SOBJ  = svout.o  smout.o
 
diff -urN original/ARPACK/UTIL/arscnd.f ARPACK/UTIL/arscnd.f
--- original/ARPACK/UTIL/arscnd.f	1970-01-01 01:00:00.000000000 +0100
+++ ARPACK/UTIL/arscnd.f	2005-06-30 23:22:55.982983811 +0200
@@ -0,0 +1,36 @@
+      SUBROUTINE ARSCND( T )
+*
+      REAL       T
+*
+*  -- LAPACK auxiliary routine (preliminary version) --
+*     Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
+*     Courant Institute, Argonne National Lab, and Rice University
+*     July 26, 1991
+*
+*  Purpose
+*  =======
+*
+*  ARSCND returns the user time for a process in seconds.
+*  This version gets the time from the system function ETIME.
+*
+*     .. Local Scalars ..
+      REAL               T1
+*     ..
+*     .. Local Arrays ..
+      REAL               TARRAY( 2 )
+*     ..
+*     .. External Functions ..
+      REAL               ETIME
+      EXTERNAL           ETIME
+*     ..
+*     .. Executable Statements ..
+*
+
+      T1 = ETIME( TARRAY )
+      T  = TARRAY( 1 )
+
+      RETURN
+*
+*     End of ARSCND
+*
+      END
diff -urN original/ARPACK/UTIL/second.f ARPACK/UTIL/second.f
--- original/ARPACK/UTIL/second.f	1995-10-02 19:56:44.000000000 +0100
+++ ARPACK/UTIL/second.f	1970-01-01 01:00:00.000000000 +0100
@@ -1,36 +0,0 @@
-      SUBROUTINE SECOND( T )
-*
-      REAL       T
-*
-*  -- LAPACK auxiliary routine (preliminary version) --
-*     Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
-*     Courant Institute, Argonne National Lab, and Rice University
-*     July 26, 1991
-*
-*  Purpose
-*  =======
-*
-*  SECOND returns the user time for a process in seconds.
-*  This version gets the time from the system function ETIME.
-*
-*     .. Local Scalars ..
-      REAL               T1
-*     ..
-*     .. Local Arrays ..
-      REAL               TARRAY( 2 )
-*     ..
-*     .. External Functions ..
-      REAL               ETIME
-      EXTERNAL           ETIME
-*     ..
-*     .. Executable Statements ..
-*
-
-      T1 = ETIME( TARRAY )
-      T  = TARRAY( 1 )
-
-      RETURN
-*
-*     End of SECOND
-*
-      END